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I. Introduction 

As early as the mid— nineteen twenties, when scientists were just 
beginning to accept the concept of the extraterrestrial origin of the 
radiation which discharged electroscopes carried aloft in manned balloons, 
and even when knowledge of the cosmos was still quite vague and unsettled, 
James Jeans (1925) speculated that the newly discovered cosmic radiation 
filled the entire universe. After two decades and a half, Richtmeyer and 
Teller (1949) proposed the other extreme kind of model wherein cosmic rays 
are supposed to be confined to the solar neighborhood; this model was 
masterfully developed later by Alfven (1949, 1950) when little was yet 
known of the nature of this radiation. Between then and now, there has 
been a spate of astounding and far reaching observational information on 
cosmic rays and astronon^. This in turn has stimulated the proposal and 
development of a variety of models and theories for the confinement and 
propagation of cosmic rays. 

For any model and theory of propagation of cosmic rays to be acceptable, 
it should be able, among other things, to explain satisfactorily the rele- 
vant observational data on cosmic rays. These include; (i) the chemical 
and isotopic composition of cosmic rays as a function of energy; (ii) the 
flux and energy spectrum of the individual nucleonic components; (iii) the 
flux and energy spectrum of the electronic component; (iv) the cosmic ray 
prehistory; and (v) the degree of isotropy in their arrival directions as 
a function of energy. The model should also be able to lead naturally and 
quantitatively to the observed intensities and spectra of electromagnetic 
radiations arising from interactions of cosmic rays with particles and 
fields in the space traversed by them such as pionic gamma rays from 
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nuclear Interactions of cosmic rays with matter and synchrotron radiation 
emitted by electrons spiralling along the feeble magnetic fields. As at 
present the model which has been able to bring to pass the greatest measure 
of success in this respect is the galactic confinement model. Proposals 
have also been made that cosmic rays may be confined to a cluster of 
galaxies of which our own, the Milky Way, is a member. 

1.1 Models of Cosmic Ray Confinement 

Of the various possibilities listed above, the solar model can be 
ruled out as untenable because of the many insuperable obstacles it en- 
counters. Among others they include: (i) the anti correlation between the 

cosmic ray intensity observed near the Earth and the 11-year solar activity 
cycle; (ii) the disagreement between the elemental abundance of cosmic 
rays and the solar composition; (iii) the significant difference between 
the steep spectral shape of particles emitted at times of solar flares and 
the relatively flat energy spectrum associated with cosmic rays at quiet 
times; (iv) the existence in cosmic rays of a continuation of the power 
law energy spectrum from a few GeV up to 10^^ eV and extending even many 
orders beyond; (v) the difficulty of explaining simultaneously the few 
g.cm”^ of matter that cosmic rays are known to have traversed before reaching 
us and the spectral shape of electrons; and (vi) the fact that we know from 
studies of the diffxise galactic radio noise, that the intensity of energe- 
tic cosmic ray electrons in interstellar space is of the same magnitude as 
is sampled near the Earth. 

Under extragalactic confinement models there are three important varia- 
tions: (a) the xjniversal model (Jeans, 1925; Lamaitre, 1931; Burbidge and 
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Hoyle, 1964); (b) the super cluster model (Burbidge, 1962); and (c) the 
local cluster model (Sciama, 1962) . By their very nature these models are 
associated with considerable amount of speculation arising from the severe 
lack of information on objects and space beyond our galaxy. On the other 
hand these models derive their strength from recent observations on the 
frequent occurrence of violent extragalactic events releasing colossal 
amounts of energy in energetic particles and hard quanta. The most impor- 
tant arguments against extragalactic models are the following: (i) The 

flux of hard quanta expected to result from the inverse Compton scattering 
of cosmic ray electrons residing within the confinement region (having an 
intensity same as in near- inter stellar space) , with photons associated with 
visible light and the universal black body radiation at 2.7®K, can be cal- 
culated as a function of the dimension of the confinement region. This 
can then be compared with the total observed flux of the isotropic compo- 
nent of gamma rays of energy in excess of 100 MeV (Fichtel et al. , 1973; 
Kraushaar et al, , 1972) to set an upper limit for the radiating distance 
which is found to be only 'VlO^^ cm (Daniel and Stephens, 1970); this dimen- 
sion may be compared with the diameter of our own galaxy which is 'V-IO^^ cm, 
the dimension of '\^10^^ cm for the local cluster of galaxies and of 'vlO^^ 
cm for the superclus ter. (il) In a similar fashion if one attributes the 
entire flux of the observed isotropic X-rays (Metzger et al, , 1964) to 
the inverse Compton scattering of electrons filling all space, with pho- 
tons of the universal black body radiation (Felten and Morrison, 1966), 
one notes that the flux of electrons needed is a factor of 40 or perhaps 
even 10^ smaller (Longair, 1970) than what is known to exist in interstel- 
lar space, (iii) A study of the relative Intensities of Be and B nuclei 
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(Brown et al. , 1973a; O'Dell et al, , 1973; Webber et al. , 1973a) and that 
of the spectral shape of cosmic ray electrons (Daniel and Stephens, 1967) 
strongly suggest that the residence time of cosmic rays in the storage 
space is leas than 10^ years whereas typical time scales involved in ex- 
tragalactic models is about two orders of magnitude larger, (iv) From a 
detailed exasiination of the energetics involved in extragalactic models, 
Ginzburg and Syrovatskii (1964, 1966, 1968) have come to the conclusion 
that the energy density of cosmic rays in metagalactic space is likely to 
be exceedingly small compared to that in our galaxy. In spite of all this, 
recently there have been attempts for the revival of the extragalactic 
model for cosmic rays though not on an universal scale (Setti and Woltjer, 
1971; Brecher and Burbidge, 1972). 

It is thus seen from the above that one is able to provide decisive 
arguments to reject the solar model of confinement of cosmic rays. On the 
other hand, though the situation regarding extragalactic models, particularly 
the cluster models, is not that decisive, the many individual arguments 
add up and pose forcible enough objections to consider them, if at all, as 
poor substitutes. This leaves only the Galactic Model for serious consid- 
eration for the propagation and confinement of cosmic rays, a model that 
has been able to explain eminently well all the observations so far avail- 
able, and one which has found acceptance from the widest circle of cosmic 
ray physicists and astrophysicists; it is in brief the best working hypo- 
thesis we have at present. 

1.2 The Galactic Model for Cosmic Ray Confinement: 

The most direct evidence for cosmic rays pervading the entire galactic 
space stems from the observation of the diffuse galactic radio noise which 
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is now universally accepted as due to synchrotron radiation emitted by 
cosmic ray electrons spiralling along the weak interstellar magnetic field 
lines; if there are electrons, there is no reason to doubt that energetic 
nuclei also exist therein. We also have equally strong evidence to demon- 
strate that estimates on the energetics and number densities of potential 
cosmic ray sources such as super novae and pulsars in the Galaxy are ade- 
quate to explain the cosmic ray intensities, as also its other features 
observed near the Earth. All these suggest a simple model in which cosmic 
rays are generated in galactic objects and injected into interstellar space; 
this radiation then propagates in the Galaxy until individual particles 
are removed from the beam through catastrophic processes or their leakage 
out of the Galaxy. While traversing galactic space, these particles under- 
go a variety of modifications, further to those which would depend crucially 
on the properties of interstellar medium such as matter, magnetic field 
and radiation field existing therein. Thus a careful and detailed study 
of cosmic rays reaching the vicinity of the Earth would, in principle, 
be capable of revealing some of them. 

The early observations on the diffuse galactic radio noise seemed to 
suggest that the emitting region had a quasi-spherical shape consisting of 
an intensely radio emitting Disk encompassed by a weakly radiating Halo 
(Shklovsky, 1952; Baldwin, 1955b). It then seemed very attractive to con- 
sider this entire volume as the dwelling region of cosmic rays; such a 
model was capable of explaining all the observations available then, in- 
cluding the high degree of isotropy of the radiation and the few g.cm"^ 
of matter traversed. However, it can easily be shown that the latter obser- 
vation needs a cosmic ray residence time of about 10® years in the halo 
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model. Furthermore detailed study on the background galactic radio emission 
suggests that the radio spurs, which were once thought to be large scale 
features extending right into the halo region could only be local phenomena 
(Baldwin, 1967). If this were so then the genuine galactic radio emitting 
region will also be considerably flattened and approach closer to the size 
of the radio disk. Furthermore, recent observations on cosmic rays strongly 
suggest that the residence time of cosmic rays is likely to be less than 
10^ years (Brown et al. , 1973a; Webber et al. , 1973a; O'Dell et al. , 1973; 
Daniel and Stephens, 1967). These observations would therefore indicate 
that even in the galactic model the effective cosmic ray storage region 
is the galactic disk. Possible intermediate regions between the galactic 
disk and galactic halo have also been referred to in the past (Daniel and 
Stephens, 1970). Thus in summary the current status in the field seems to 
be in favour of the galactic disk model though one could still argue about 
the effective extend of the thickness of the Disk for cosmic ray containment. 
1,3 Scope of the Present Review ; 

It is evident that if one is to investigate and deduce useful informa- 
tion from a model in which cosmic rays are generated, propagated and stored 
in our galaxy, one would have to make the best use of all available informa- 
tion on the Galaxy, its dimensions and its constituents. In particular, 
one would need as much data as possible on the interstellar medium and its 
physical state. We have, therefore, judiciously collected from the relevant 
literature all such information of direct concern to us here and summarized 
them in Chapter 2. Since in this model, the confinement region is adequately 
circumscribed, one is able to work out a theoretical formulation of the 
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propagation of cosmic rays in the Galaxy; this is described in Chapter 3. 

The present theoretical treatment is however limited to a situation in 
which the cosmic rays are in a state of equilibrium and are distributed 
homogeneously in the residence volume. While such a simplified treatment 
has been forced upon us because of the otherwise intrinsic complexities 
and difficulties encountered, the assumptions regarding the equilibrium 
nature of the cosmic rays and its homogeneous distribution in the Galaxy 
seem to be well founded on the basis of our existing knowledge. In the 
fourth Chapter, we examine and interpret the observational data on cosmic 
rays in relation to their propagation in the Galaxy. 

The fact that cosmic rays traverse a region of cosmic space pervaded 
by a tenuous gas , weak magnetic fields and radiation fields , implies that 
cosmic rays which are literally Immersed in them, should interact with 
them in a wide variety of ways. Many of these lead to the production of 
radiations ranging the entire electromagnetic spectrum - from radio waves 
to gamma rays; these are described in Chapter 5. Finally, recent studies 
in astrophysics have unmistakably highlighted the possible roles that 
cosmic rays are likely to play in some aspects of galactic dynamics such 
as the hydrostatic equilibrium of interstellar gas, gravitational insta- 
bility and formation of clouds, heating of Interstellar gas etc. These are 
briefly covered in Chapter 6. Though we still have a long way to go 
towards a fuller understanding of the propagation of cosmic rays in the 
Galaxy and its implications, we can already see that much progress has been 
achieved particularly during the last decade or so. 

Before proceeding further, we would like to emphasize the propagation 
of cosmic rays in the Galaxy is a subject in which theoretical 
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formulations including model building on the one hand, and experimental 
observations on the other, have to advance in parallel, continuously 
checking the findings of one with that of the other, and incoporating 
suitable modifications in the model and the theory to build-up an in- 
creasingly self consistent picture. Real progress is often made in being 
able to rule out or discredit models whose predictions are in contradiction 
with observations rather than in being able to propose new models which 
are consistent with existing observations. Furthermore, in this subject, 
wherein one finds room for highly speculative ideas , there is immense 
scope for widely differing viewpoints. It is therefore too much to expect 
in one review article extensive coverage to be accorded to all view points 
so far proposed; a certain amount of selectivity on the part of the authors 
seems inevitable to make the presentation coherent. The reader* s indul- 
gence and understanding are therefore solicited for any apparent biases 
and omissions that might have resulted in an attempt to achieve this. 
Finally, as supplementary reading we will like to refer the reader to a 
number of other review articles, and references therein, on topics of 
allied and overlapping Interest (Parker, 1969; Meyer, 1969; Daniel and 
Stephens, 1970; Shapiro and Silberberg, 1970; Simpson, 1971; Wentzel, 1974). 
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2. GENERAL DESCRIPTION OF THE GALAXY 

The Galaxy is an assembly of about 10^^ stars together with much gas 

and dust. Its spheroidal central region is enveloped by a flat disk of 

stars in which the spiral arms are embedded. The Galaxy is in a state of rotation 

main body 

about an axis through its center perpendicular to the galactic plane. Around this/ 
is a slowly rotating near-spherical halo containing a low density of high 
velocity stars and globular clusters. The Sun is situated close to the 
equatorial plane of the Galaxy and about two thirds of the distance from 
the centre to the edge of the Disk, In this Chapter* we shall summarize 
the general properties of the Galaxy relevant to our present purpose as 
inferred from observations to-date, 

2,1. Dimensions and Mass Density 

The Galaxy has no sharp physical boundaries , and its dimensions are 
not known with high precision. The distance Rq between the Sun and the 
Centre is taken to be 10 kpc (Arp, 1965; Schmidt, 1965), though this value 
could be as small as 8.5 kpc (Toomre, 1972). However, in this Chapter all 

parameters which depend upon Rq are given on the basis of Rq = 10 kpc. 

From the neutral hydrogen distribution, one can deduce the diameter of the 
Galaxy covering the major features (but excluding the extended Outer Arm), , 

as about 30 kpc. The thickness of the Disk varies from region to region, 
the optical disk being flatter than the radio disk. A study of the dis- 
tribution of the stellar population in a direction perpendicular to the 
galactic plane shows that the effective thickness defined as the distance 
between half density points is different for different classes of stars; 
for Population I Cepheids it is about 180 pc (Walraven et al, 1958) while 

for F8-G8 dwarfs it is about 300 pc (Elvius, 1965).” Also the thickness 
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increases considerably to more than about 1 kpc towards the Centre (Oort, 
1938), The non-thermal radio emission in the Galaxy is symmetric about 
the Centre within a radius of 8.7 kpc having a thickness of 750 ± 100 pc 
(Baldwin, 1967). This would mean that the Sun is outside this symmetric 
disk of emission. Beyond the 8.7 kpc region there is an outer disk of 
uniform emission which is relatively thinner with a probable thickness of 
about 500 pc. 

According to Schmidt (1965) the density distribution in the Galaxy 
can be well approximated by a spheroid having an axial ratio 1:20 with 
the following features; (i) equidensity surfaces in the form of similar 
spheroids, (ii) a mass point at the Centre and (iii) a shell outside. The 
mass distribution could then be represented as: 

p = 3.930 R~^ - 0.02489R for R < 9.71 kpc 


2.1 


p = 1449.2 R-** 


for R > 9.71 kpc 


where p is expressed in solar mass M per pc^ and R is the distance from 
the Centre in kpc; the central mass point has a value of 7 x 10^ that 
of the spheroid is 8.2 x 10^® and of the shell 9.3 x 10^° From 

this the total mass of the Galaxy is 1.8 x 10^^ M , half of which is within 

O 

a spheroidal surface drawn through the Sun; the total mass in the form of 

gas is only about 5 to 7 x 10^ M . Observationally , the hypothetical cen- 

0 

tral mass point consists of about 10^ within a radius of 100 pc, 10® 

within 10^ M within 1 pc (Oort, 1973). The density near the Sun is 0.145 

O 

M pc“® i.e. , 10”2 3 g cm”®. Further one may also add that the escape velo- 

0 

city near the Sun is nearly 380 km s”^;the age of the Galaxy is about 1.2 
X 10^® yrs; and the period of rotation in the neighborhood of the Sun is 


2,5 X 108 yrs. 
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2 , 2 S tructure 

Much of the information on the structure of the Galaxy comes from the 
study of the 21 cm emission profiles of the Galaxy. In Figure 2,1 is shown 
the neutral hydrogen distribution in the galactic plane. This figure is 
derived from Kerr and Westerhout (1965) and Kerr (1967) , but the scale of 
the radius has to be enlarged by a factor of 10/8.2 to agree with the new 
value of = kpc. In this figure the region marked 1 in the centre is 
the nuclear disk of radius about 900 pc with an inner and outer section, 
the inner one being a rapidly rotating disk and the outer one, probably, 
with some radial motion as well. The region marked 2 constitutes the 
"3-kpc Arm" at about 4 kpc from the centre; it has an outward motion to- 
gether with rolling. The total neutral hydrogen mass in the Nucleus is 

6 7 

estimated to be about 3 x 10 I^while in the 3 kpc Arm it is 3 x 10 

(Rougoor and Oort, 1960), At about 5 kpc from the centre is the "Scutum- 
Norma Arm" (region marked 3 in Figure 2.1); the "Sagittarius Arm" (region 
No. 4 is at a distance of 7.5 kpc. The Sun is situated in the inner edge 
of the "Orion Arm" or the "Carina Cygnus Arm" (region No. 5) which is about 
10 kpc from the Centre. At about 3 kpc beyond the Sun, we have the promi- 
nent "Perseus Arm" (region No, 6). Outside the Perseus Arm one can see 
the "Intermediate Arm" (region No. 7) beyond which is located the faint 
"Outer Arm", (region No. 8). The spiral feature of the Galaxy can in 
general be represented by a logarithmic multi-arm spiral with a pitch 
angle of about 85® (Pskovskii, 1965). It is difficult to decide from 
Figure 2.1 whether a single spiral arm continues over very large distances 
in the Galaxy or whether the observations can be interpreted as a collec- 
tion of short spiral arms. The first alternative will suggest that the 




Fig, 2,1 : Distribution of neutral hydrogen in the Galaxy according to Kerr and Westerhout (1965) 

in which the "3 kpc Arm” and the nucleus have been inserted from Kerr (1967). The scale of radius 
corresponds to Ro » 8,2 kpc. The inset near the bottom indicates the value of ny in atoms cm”^, which 
does not apply to the central region. 
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spiral arms exist over many rotation periods of the Galaxy while if the 
latter were the case, it would indicate that they form and dissolve even 
in a time duration as short as the galactic rotation period. 

Many spur features have been noticed from studies of the galactic 
non-themal radiation extending from the equator to high latitudes. Of 
the three which have been well recognized, the "North Polar Spur" is the 
most intense; the other two are the Cetus Arc and the Loop III. The radio 
brightness of these features falls off with increasing latitude; at the 
same time not much is known about their brightness at low latitudes, where 
they merge into the intense radiation from the Disk. Furthermore, although 
at low latitudes they appear to lie roughly perpendicular to the galactic 
equator, there is no evidence to reveal whether they cross the galactic 
plane. There are many ways of understanding the origin of these spurs: 

(a) old supernova remnants of large size relatively close to the Sun (Han- 
bury Brown et al, , 1960) ; (b) cosmic ray pressure in the Disk can, according 
to Parker (1965) , lead to instabilities in the galactic magnetic fields 
whereby loops of field may be blown out to large values of z to make them 
appear like long fingers roughly perpendicular to the galactic plane; (c) 
effect of helical tubes of magnetic field wrapped round the local spiral 
arm (Rougoor* 1966); and (d) large scale features in the galactic Halo, 
Nevertheless, there are as yet no decisive arguments to establish any one 
of the above possibilities nor to rule out others (Baldwin, 1967). 

Apart from the spurs, which are mostly high latitude features of 
radio continuum in the Galaxy, one can also identify many structures from 
the neutral hydrogen distribution in the Disk. The most prominent of 
these seen from Figure 2,1 are: (1) the extension of .the Orion Arm. towards 
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the Anti Centre; (ii) the outer strong feature in the longitude range 200° 
to 240° as the extension of the Persus Arm; and (iii) the connecting links 
between the Nucleus and the 3 kpc Ann. Many inter arm links, though not 
very prominent, can also be recognized providing support for the assumption 
that for many purposes, e.g. the cosmic ray propagation, the Galaxy is a 
simple uniform disk with many large and small scale structural features. 

It is also pertinent to point out that the hydrogen layer beyond R = 9 kpc 
shows a systematic distortion with a positive deviation in the region of 
galactic longitudes 20° < 1 < 140° and a negative deviation at 200° < 1 < 340° 
the deviation is nearly maximum in the direction perpendicular to the Sun- 
Centre direction (Kerr and Westerhout, 1965). It is suggestive that this 
distortion is a general feature of the outer region of the Galaxy. 

The association with the Galaxy of globular clusters well beyond 
z 'v2,5 kpc makes one surmise that the Galaxy extends also in the z-direction 
beyond the apparent boundary of the Disk mentioned in Section 2.1. More- 
over, there is evidence for the existence of neutral hydrogen at distances 
of 1 kpc above the Perseus Arm with an estimated density of nearly 3 x 10”^ 
hydrogen atoms cm“^, decreasing to about 10“^ hydrogen atoms cm“^ at 1.5 
kpc (Kepner , 1970). Also, from a study of the negative velocity clouds 
at high latitudes, one infers that there should exist considerable amotint 
of gas outside the Disk so as to offer an efficient mechanism to retard the 
inf ailing of the intergalactic gas (Oort, 1969). All these indicate the 
existence of a physical halo around the Disk of the Galaxy in which most 
of the gas is likely to be highly ionized. This deduction is further 
strengthened by the fact that it is difficult to imagine a sharp boundary 
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between the optical disk and the near-intergalactic space which differs 
markedly in the physical conditions, namely matter density and magnetic 
field (Ginzburg, 1967). The shape of the Halo may be nonhomogeneous 
spheroid, with no sharp boundary but perhaps extending up to about 10 kpc 
from the plane of the Disk. The non- thermal radiation at high latitudes 
also suggests the existence of a radio halo with a radius probably of 
about 18 kpc (Baldwin, 1967). However, the absolute brightness from the 
Halo is in question; this aspect will be further examined in Section 5,3, 
Notwithstanding all these, a detailed study of 18 spiral galaxies shows no 
evidence of any associated large radio halo postulated for our own galaxy 
(Mills, 1967). 

2.3 Interstellar Gas 

The existence of interstellar gas was first suggested by Hartmann (190A) 
following his observation of stationary absorption lines of ionized calcium 
in a spectroscopic binary. Today, the basic information on interstellar 
matter is derived from the intensity of neutral hydrogen emission line at 
21 cm. If as is generally believed atomic hydrogen is the major constituent 
of interstellar gas, one can reliably assign values for the density of gas 
in different regions of the Galaxy from the 21 cm observations. Following 
such a consideration it is seen that the distribution of the mean density 
in the galactic plane has a broad maximum of about 0.6 hydrogen atoms cm“^ 
between R = 5 and 12 kpc falling off both inwards and outwards (Kerr and 
Westerhout, 1965). Nevertheless the density of neutral hydrogen in the 
central region increases from about 0.2 atoms cm“^ at R = 600 pc to about 
5 atoms cm"^ at R - 100 pc (Oort, 1973). The density of gas may further 
increase towards the inner parts, where the observed density of ionized 
hydrogen itself is about 15 atoms cm“ (Lequeux, 1967). The ionized 
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hydrogen constitutes only about 3% of the total gas in interstellar space. 

The mean density of gas in the Halo, which is fully ionized is about (5 - 10) 

X 10”^ atoms cm"*^. The thickness of the layer of hydrogen gas obtained 
from a series of tangential measurements gives a value of about 280 pc be- 
tween half density points. However, around 4 kpc from the centre the 
thickness is only about 150 pc, which drops to a value of about 100 pc 
close to the Centre. This demonstrates that the total mass in the form 
of gas in the central region of the Galaxy is negligible and that the cen- 
tral mass density resides essentially in the stars. At large radii, the 
thickness of the gas layer increases considerably; it may be as high as 
several hundreds of parsecs beyond 15 kpc. The distribution of gas per- 
pendicular to the galactic plane is symmetric and is approximately of 
Gaussian shape, except at large distances where it is wider. The mean 
thickness of 280 pc between half density points corresponds in the z-component 
to a dispersion of = 120 pc. Mention may also be made here that the 
distribution of stars in the z-plane of the Galaxy in the neighborhood of 
the Sun is also approximately a Gaussian; and the Sun is situated at about 
10 pc north of the equatorial plane of symmetry (Elvius, 1965). 

Although the large scale physical properties of interstellar gas are 
known fairly well, the detailed structure could be very complex. The in- 
terstellar medium is generally known to be composed of cool, dense clouds 
embedded in a hot, rarefied intercloud medium. The clouds vary widely in 
size from about 0.1 pc to 100 pc (Van de Hulst, 1958) ; their density and 
distribution in the Galaxy is not uniform. Combining the observational 
data as summarized by Van Woerden (1967) , and the theoretical study on the 
pressure equilibrium between the clouds and the intercloud medium (Field 
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et al., 1969), one can describe the physical properties of an average 
cloud as follows: the diameter of the cloud is about 10 pc; the density 

of neutral hydrogen (20-30) atoms cm'^; the free electron density 

n^ 2 X 10’^ el. cm"^ (Stephens, 1971); and the temperature is about 100°K. 

The fraction of volume occupied by these clouds in the equatorial plane of 
the Galaxy is about (2-4)%. The distribution of clouds in the z plane of 
the Galaxy is similar to that of neutral hydrogen. The total mass in the 
form of clouds is about 70% of the interstellar gas. The in ter cloud medium 
which is considered to fill the Galaxy uniformly, has a density in the 
equatorial plane of n^ ^ 0.2 atoms cm~^ at a temperature of about (6 - 10) 

X 10^®K with a free electron density n^ 4 x lO'^el. cm”^, as inferred 
from the dispersion measures of pulsars (Prentice and Ter Haar, 1968). The 
distribution of this hot gas in the z plane is slightly wider than that of 
cold gas with a dispersion o = 160 pc. As can be seen from the free elec— 
tron densities in the clouds and intercloud medium, the fraction of ionized 
gas is about 3 x 10~^. 

Except for a few elements, the chemical composition of interstellar 
matter has not been determined directly from observations. The estimate 
made on the strength of interstellar absorption lines, though somewhat less 
certain, reveal that the composition of gas is about the same as in Popu- 
lation I stars (Spitzer, 1968a). The relative abundances of various ele- 
ments as summarized by Greenberg (1968) and Allen (1965) are given in 
Table 2.1. 

Recent observations by the Copernicus satellite indicate (Morton et al. , 1973 
Jenkins et al. , 1973) that the abundances of elements like C, 0, Si and Fe 
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Table 2.1 

Relative Abundances of Some of the Elements in Interstellar Mediima 
Relative Relative 


Element 

number 
of atoms 

Relative 

mass 

Element 

number 
of atoms 

Relative 
■ mass 

U 

1.0 



1.0 



Si 

2.5 

X 

10-5 

7.0 

X 10- 

He 

1.2 

X 

10-1 

4.8 

X 

10“ 1 

P 

3.2 

X 

10“ 7 

9.8 

X 10"5 

Li 

1.0 

X 

10" 5 

7.0 

X 

lO-O 

S 

2.8 

X 

10“ 5 

9.0 

X 10“ 

Be 

2.5 

X 

10-10 

2.3 

X 

lo-o 

Cl 

1.0 

X 

10" 5 

3.5 

X 10" 5 

B 

6.3 

X 

10-10 

6.8 

X 

10“5 

Ar 

7.6 

X 

10-6 

3.0 

X 10“"+ 

C 

2.5 

X 

10-*^ 

3.0 

X 

10-5 

K 

7.9 

X 

10-5 

3.1 

X 10“ 5 

N 

1.3 

X 

10“** 

1.7 

X 

10-5 

Ca 

1.6 

X 

10" 5 

6.4 

X 10-5 

0 

7.9 

X 

lO"'* 

1.3 

X 

10“2 

Ti 

5.0 

X 

10-5 

2.4 

X 10-5 

F 

2.5 

X 

10“7 

4.8 

X 

10-5 

Cr 

8.9 

X 

10-5 

4.7 

X 10-5 

Ne 

6.3 

X 

lO-** 

1.2 

X 

10“2 

Mu 

6.3 

X 

10" 5 

3.5 

X 10"5 

Na 

1.3 

X 

10-6 

2.9 

X 

10-5 

Fe 

2.5 

X 

10" 5 

1.4 

X 10“ 5 

Mg 

4.5 

X 

10-5 

1.1 

X 

10" 5 

Co 

5.4 

X 

10" 5 

3.2 

X 10"5 

A1 

1.8 

X 

10-5 

4.8 

X 

10" 5 

Ni 

1.1 

X 

10-5 

6.6 

X 10" 5 
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vary from region to region, and that some of these trace elements are 
under abundant compared to the universal abundances. This can be due to 
the locking up of these elements in dust grains (Spitzer, 1968b) ; such 
depletion of trace elements in clouds resulting from accretion on grains 
seems to gain support from the fact that the temperature of the clouds 
are higher than that predicted using normal abundances (Goldsmith et al. , 
1969) , Since energetic cosmic ray particles while traversing matter in 
space do not distinguish very much between individual atoms and small 
grains, pne can, for purposes of cosmic ray propagation, consider the 
cosmic abundances in Table 2.1 as representative of the interstellar gas. 
Apart from individual atoms, the interstellar gas also contains molecules. 
The abundance of molecular hydrogen is relatively low in common clouds and 
high in dense dark clouds; however, from observations one can only say 
that the total abundance in the Disk hardly exceeds that of atomic hydro- 
gen (Pikelner, 1967), During more recent years, molecules like OH, CH, 

CN and even more complex ones are being observed in interstellar space. 

As pointed out earlier, a small fraction of the gas also exists in the 
form of dust grains which are presumably formed in the clouds. It has 
been suggested that the interstellar extinction can be well understood if 
grains consist of graphite cores with ice mantles, while the intrinsic 
polarization can be ynderstpod as due to small graphite flakes (Greenberg, 
1967) . 

2,4 Magnetic fields 

The first observation on stellar polarization (Hiltner, 1949; Hall, 
1949) almost certainly required for its understanding the existence of 
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interstellar magnetic fields. At present there exist many different 
methods of determining the strength and orientation of the magnetic field 
in the Galaxy. They are; (a) Stellar polarization; (b) Faraday rota- 
tion; (c) Zeeman splitting; (d) radio polarization; and (e) radio bright- 
ness distribution. We shall now very briefly summarize the results from 
these studies. 

(a) Stellar polarization 

It Is generally accepted that stellar polarization is produced by 
the scattering of light by grains, which get aligned by the galactic mag- 
netic field. Early measurements by Hiltner (1949) and Hall (1949) showed 
that the magnetic field is predominantly parallel to the galactic plane 
and that the major component lies along the spiral arm, A helical field 
in the spiral arm has also been postulated (Hoyle and Ireland, 1961). 

From an extensive analysis of recent data Mathewson, (1968) has proposed 
a model in which the local field is wound in a spiral of pitch angle 7® 
arotind the surface of an elliptical cylinder, whose axis is in the galactic 
plane and at right angle to the direction of the centre; the spiral is 
sheared through 40® anti-clock wise as seen from the north galactic pole. 

(b) Faraday rotation 

A plane polarized radio wave of wavelength X, while traversing an 
ionized medium with a frozen- in magnetic field, undergoes a rotation of the 
plane of polarization. The rotation measure R = ijj /x^, where the angle of 
rotation is in radians and X is in meters, is numerically expressed as 
R * 8.1 X 10^ Bi 1 n ds, where Bii is the longitudinal component of the 
magnetic field in Gauss and ds is the element along the line of sight in 
pc. When the rotation measure obtained fot extragalactlc sources is plotted 
against the galactic latitude of the respective sources (Berge and S dels tad. 
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1967), it is seen that most of the rotation takes place in the Galaxy, 

From these observations, taking <n 10”^ towards the 

galactic pole, and substituting n^ = 3 x 10”^ el, cm~^ as the mean electron 
density for L = 300 p in the Disk across the galactic plane, one finds 
Bu 2yG as the mean longitudinal component of the local magnetic field 
in the Galaxy (Shakeshaft, 1968), Rotation measures from pulsars also 
indicate that Bji is about 3pG (Smith, 1968). These values imply that the 
strength of the galactic magnetic field would be ^3yG. In the light of 
rotation measures obtained from extragalactlc sources , a revised analysis 
of stellar polarization suggests (Mathewson and I^icholls, 1968), that the 
helical component of the magnetic field, needed to account for the stellar 
polarization, is superimposed on the large scale longitudinal field in the 
Galaxy, which is along the spiral arm, 

(c) Zeeman Splitting 

In the case of Zeeman splitting of excited lines , the separation be- 
tween the two components polarized circularly in opposite directions is 
given as Av = 2.8 x 10^ where is expressed in Gauss and v in Hz. 

Recent observations of Zeeman splitting of neutral atomic hydrogen line at 
21 cm, both in absorption against strong radio sources and in emission, 
towards different galactic directions, show a wide variation in the magni- 
tude of Bjj (Verschuur, 1969). It has been pointed out by Verschuur that 
these observations cau reveal information only about the fields in high 
density, low temperature H-1 clouds and not about interstellar fields in 
general; it is also suggested that the cloud field can be effectively de- 
tached from the more general field in interstellar space. Thus one finds 
that a detailed study of the Zeeman splitting would enable one to elicit 
information about magnetic fields in dense clouds existing in interstellar 


space. 
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(d) Radio Polarization 

Since the non-thermal coc^jonent of the radio continuum is due to 
synchrotron radiation emitted by relativistic electrons, one would expect 
that polarization measurements would shed light on the orientation of 
large scale magnetic fields in the Galaxy responsible for the emission. 
However, the Faraday rotation within ionized clouds tends to depolarize the 
radiation and hence one is unable to derive Information on the magnetic 
field without a knowledge of the exact distribution of ionized hydrogen 
along the line of sight. Observations in the local region of space, where 
the line of sight crosses not more than a few clouds of ionized hydrogen, 
suggest that the field is well aligned to the local arm (Lequeux, 1969). 

(e) Radio Brightness Distribution 

Information on the strength of large scale magnetic fields in the 
Galaxy comes from the unique relationship between the cosmic ray electron 
spectrum in Interstellar space and the spectral distribution of the radio 
continuum -along any given direction in the Galaxy. This method was first 
introduced by Biermann and Davies (1960) after which many such investigations 
have been made in recent years. The mean perpendicular component of the 
magnetic field in the Galaxy averaged over the line of sight in different 
directions in the galactic plane is 6pG (Section 5.3), This value com- 
bined with that derived in (b) above on the longitudinal component leads 
to a field strength of (6-8) pG for the large scale magnetic fields. The 
magnetic field in the Halo is probably about (2-3) pG. 

2.5 Radiation Fields 

Radiation fields existing in the Galaxy extend right from the radio 
to the gammaray domain. Part of this radiation is the result of diverse 
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processes occurring inside the Galaxy; the rest has its origin in the 
Universe at large. We shall briefly sunnnarize here our knowledge of the 
radiation energy density in the Galaxy in the different spectral regions, 
(a) Radio Region 

The electromagnetic radiation in this regime is essentially due to 
synchrotron radiation emitted by cosmic ray electrons in the Galaxy, At 
frequencies below a few Ifilz, the free-free absorption in interstellar space 
prevents this radiation from propagating through the Galaxy , while for the 
entire range in excess of this frequency, galactic space is transparent. 

If the differential spectral shape of the radio brightness in the sky is 
expressed as I a it is then found that a 0.0 at about 2 MHz slowly 

increasing to about 0.8 at frequencies above 2 x 10^ MHz; below 2 MHz, the 
spectrum has a positive slope due to interstellar absorption. The total 
energy density in this field is about 10 ^ eV cm“^ in the neighborhood 
of the Sun. 


(b) Microwave and Far Infra-red Regions 

This is the most important region of the radiation field existing in 
space. Considerable enhancement of the sky brightness in the microwave 
region far above that expected from a simple extrapolation of the galactic 
radio continuum was first discovered by Penzias and Wilson (1965) . This 
has been interpreted by Dicke et al. (1965) as due to the existence of the 
Universal Black Body Radiation, the inescapable remnant of the "big-bang" 
origin of the Universe (Gamow 1948) . If this were so, then this radiation 
would permeate the whole Universe; within the accuracy of measurements so 
far made on the isotropy of this radiation the universal nature of this 
radiation has been substantiated. The observed intensity in the frequency 
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region of 10^ - 10^ MHz, is consistent with the Rayleigh-Jeans part of 
the black body radiation curve at T = 2.68®^: (for details see Peebles, 

1971). If the black body nature of this radiation is correct, then the 
corresponding energy density would be 4aT**/c JxJ 4,73 x 10“^ T** = 0.25 eV. 
cm"^, where o is Stefan's constant; the mean energy of the photon <e> = 

2.7kT * 6.3 X lO”** eV, where k is Boltzmann's constant. 

(c) Infra-red Region 

Observations in the infra~red region are so far made using grotind 
based instruments; they are thus exposed to the thermal radiation-environment 
of the atmosphere. Because of this reason, determination of the extra- 
terrestrial continuum radiation becomes difficult and one is constrained 
to observations relating to sources only (Webbink and Jeffers, 1969). On 
the other hand, theoretical estimates of the infra-red radiation from inter- 
stellar grains have also been made in the past; recently Krishnaswamy (1970) 
has calculated the expected spectral distribution in the solar neighborhood 
from different grain models and from this it is found that the energy 
density in the region (3-15) x 10^ MHz is ^ 10” ^eV cm”^. The energy density 
decreases sharply at higher frequencies, 

(d) Visible Region 

The radiation field in this small band of electromagnetic spectrum 
comes from stellar emission. The interstellar radiation field calculated 
by summing up essentially the contributions from stars of various spectral 
types (Zlmmermann, 1964) is consistent with a grey body radiation at 
T 10** *K with a dilution factor of 10^**. The total energy density in 
this frequency band is about 0.5 eV cm”^, which includes a little of near 
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infra-red and near ultraviolet radiations; the mean photon energy is about 
2 eV. 

(e) Ultra Violet Region 

Our knowledge of the radiation field in this domain is again very 
poor, and measurements are yet to be made in the future. However, recent 
calculations of the radiation density in the frequency interval (1.2-3. 2) 

X 10^ MHz indicate that p . ^ x 10“^ eV cm'"^ (Habing, 1968). 

ph 

(f ) X and Qamma-ray Regions 

Reliable measurements are now available on the diffuse background of 
X-rays and gamma rays. These observations suggest that the radiation field 
above 1 keV (2.4 x 10^^ MHz) is of extragalactic origin. The total energy 
density in this do^nain is only about 4 x 10”^ eV cm""^; also it can be in- 
ferred from the spectral shape that the radiation density decreases with 
increasing energy. 

From the above summary, it becomes clear that the important radiation 
fields in the Galaxy are the Universal Black Body Radiation and the stellar 


optical emission. 



26 


3. C6sinic Ray Propagation In the Galaxy-Theoretical Aspects 

In what h^s been summarized In Chapter 1, we have identified the most 
probable region of confinement of cosmic rays to \>e the Galaxy; we then 
surveyed, in Chapter 2, those descriptions and features of the Galaxy of 
relevance to our considerations here. We now proceed to examine theories 
applicable to the problem of propagation of cosmic rays in galactic space. 
In attempting to. formulate a rigorous theory for this, one encounters 
diverse difficulties because of the many inherent uncertainties and com- 
plexities of the problem. Some of these difficulties can be lessened 
significantly if cosmic rays are assumed to be in a state of equilibrium, 
and are distributed homogeneously within the confinement volume. Indeed, 
it is found that there does exist a reasonable body of. evidence in favor of 
such an assumption though recently some workers have pointed out observa- 
tions which seem contrary to it. These aspects are first summarized in 
this Chapter. We then proceed with a brief description of models of diffu- 
sion of cosmic rays in interstellar space; their relative merits are also 
critically examined. Following this, a general equation describing the 
transformation in the chemical composition and energy spectrum of cosmic 
rays during their propagation is discussed. Finally, several theoretical 
models of propagation of nucleonic and electronic components of cosmic 
rays, often employed to understand the experimental results, are described. 
We will like to state here that in the theoretical formulation of the pro- 
pagation problem, we will not separately deal with propagation in the 
"source region" since the physical characteristics of this region, even if 
it turns out in the future to be important in the overall problem of cosmic 
ray propagation, is poorly understood at present. 
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3.1 Quasi-statlonary State and Spatial Homogenlty of Cosmic Rays 

In order to study the propagation of cosmic rays in the region of space 
where they reside, it is essential first to enquire whether cosmic rays ex- 
tend throughout this region in a quasi-stationary state or with notable 
fluctuations. Although the first alternative seems to be the conservative 
line of approach, it is widely adopted in all Interpretative studies since 
it is amenable to theoretical formulations and constructive advancement 
compared to the latter in which there are considerable difficulties and 
very many unknowns. We shall, therefore, first examine the evidence from 
cosmic ray prehistory to see whether one can justify the simple hypothesis 
that cosmic rays are in a quasi-stationary state in the Galaxy as a whole. 

Most of our information on cosmic ray prehistory is obtained from the 
analyses of stable and tinstable cosmogenic nuclei found in meteorites (Lai, 
1966). Such studies. Involving nuclei of half lives ranging from days to 
millions of years, reveal that the average cosmic ray flux, integrated over 
the mean lives of the respective radio isotopes, has remained unchanged 
(Geiss, 1963), Scheaffer et al, (1963) have shown that the mean Intensity 
over the last 500 yrs . has been the same as that averaged over 5 x 10^ yrs 
within ±10%, Furthermore, the gross cosmic ray intensity has also remained 
unchanged within a factor two over extended periods in the past up to 10^ 
yrs (Voshage and Hintenberger , 1963; Anders, 1965; Lipschutz et al, 1965). 

A similar analysis using moon rocks (Finkel et al, 1971) suggests that the 
time averaged flux of cosmic ray protons above a GeV has been the same 
during the past 5 million yrs as the contemporary one. Hence, one feels 
justified in taking the viewpoint that the cosmic ray intensity has not 
changed, over a few rotations of the Galaxy, by more than a factor of two. 
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It needs to be pointed out here that above 10^ yrs the lack of detailed 
information does not permit one to rule out possible large scale fluctua- 
tions to occur with periods which lie between the half-lives of the radio 
isotopes used. For instance, large scale fluctuations during the present 
epoch with characteristic periods of a few times 10*^ yrs will not contra- 
dict the observational deduction that the present cosmic ray intensity is 
nearly the same as that averaged over the last 10^ or 10^ yrs (Parker, 1968) 
However, the failure to detect any large scale variation so far would 
appear to favor cosmic rays to be in a state of quasi-equilibrium over long 
periods of time. 

As for spatial homogenity, a detailed examination of non- thermal radio 
continuum towards various galactic directions by Anand et al, (1968a) re- 
veals that the cosmic ray intensity in the Galaxy when averaged over the 
line of sight is the same in all regions within a factor of two. However, 
it is difficult to rule out, at this stage, possible gradients to exist 
either in localized regions of space, such as active sources and source 
regions, the spiral arms, or in the Halo perpendicular to the galactic 
plane. In fact, it has recently been suggested, on the basis of radio as- 
tronomical and cosmic gamma ray observations, that contemporary cosmic rays 
are not uniformly distributed in the Galaxy (Sign ami and Fichtel, 1974; 
Stacker et al, 1974). Nonetheless, the present status of these postulations 
are either tentative or speculative; also their quantitative effect on the 
theoretical treatment of the propagation problem is not clear. Finally, 
there are as yet no major and reliable attempts made to formulate theoreti- 
cally the propagation problem in a situation wherein cosmic rays are not 
.in a state of equilibrium in the Galaxy; it is, however, almost certain 
that it is likely to attract greater interest and attention of the re- 
searchers in the years to come. 
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In consequence of the above observations , we will confine ourselves , 
in this review, to the extensive work so far undertaken under the basic 
assumption that cosmic rays are in a state of quasi-equilibrium in the 
Gala^^, 

3, 2 Models of Diffusion of Cosmic Rays 

The motion of a cosmic ray particle in the Galaxy is influenced by 
the nature of the magnetic fields therein. In this section, we shall 
briefly describe some of the important models for the diffusion of cosmic 
rays in interstellar space. 

3.2,1 Motion of charged particles in a magnetic field and their difftision 
In its simplest form the motion of a charged particle in a magnetic 
field having large scale irregularities, is associated, with a drift velo- 
city across the magnetic field and is given as 

^d ^ 

where, 1 is the characteristic scale of field inhomogenity , r^ = E sin 0/3OOZB 
is the radius of curvature in cm of a relativistic particle of charge Z, E 
in eV is the energy of the particle moving at an angle 6 to the magnetic 
field of strength B expressed in Gauss, and v is the velocity of the parti- 
cle. Two points need to be noted here. One, for a typical situation in 
interstellar space, Equation 3.1 would require that for the escape from the Galaxy 
of particles of energy, say 1 GeV, the scale length has to be as small as 
10^^ cm. Two, this expression, to be applicable, requires that the motion 
of the particle is adiabatic all along the line of force and r^ << 1, a 
condition unlikely to be true in Interstellar situation. 
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The above difficulties can be circumvented by proposing a variation in 
the model wherein cosmic rays have to traverse scattering centers, such as 
magnetic clouds and shock waves, and that the random walk of particles re- 
sulting from multiple scattering in the medium leads to diffusion (Ginzburg 
and Syrovatskii, 1964; Parker, 1963). In case of magnetic clouds with size 
1 very much in excess of the radius of curvature of the particle r^^ inside 
the cloud (l.e. r^ « 1^) , one can write for the scattering mean free 
path 

X l.O/n.l^ 3.2 

s ^ c 

/ 

where n is the number density of clouds in space. It can then be seen that 
for typical interstellar conditions X is independent of energy at least 
up to 10^^ eV. Nevertheless, at much higher energies (i.e. r^^ > 1^), 

X depends strongly on the energy and can be written as 


^ n i 


3.3 


Since the diffusion coefficient D (1/3) X v, it is seen that when n ^5 x 

s ^ 

10“®^ cm”^ and 1 'v/ 20 pc, the value of D is x 10^^ cra^ . It can easily 

c 

be shown that such a large value of diffusion coefficient would allow cos- 
mic rays to escape from the Disk within a short period of '^lO** yrs unless 


one invokes strong reflection at the boundary for an effective confinement 
in the Disk; on the other hand, it would be quite consistent with a model 
\/jierein cosmic rays are confined to the Halo. It is important to point 
out that in these considerations one has neglected the existence of any 
ambient magnetic field in the intercloud region in which the particles are 
not likely to satisfy the relation r^ « 1 required of the treatment; 
rather, one has made use of the idea of single particle scattering. 
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Plentifiil realization has now emerged that, since the energy density 

* 

of cosmic rays is comparable to the magnetic field energy and the kinetic 
energy of the gas, the collective effect of the relativistic cosmic ray 
gas, which is essentially a collisionless hot plasma, assumes significant 
proportions. We shall now briefly examine the collective effects of cosmic 
rays on their diffusion. 

3,2.2 Self scattering of cosmic rays 

Earlier hopes that scattering of cosmic rays by small scale irregulari- 
ties in the general galactic magnetic field will be adequate for our under- 
standing of cosmic ray propagation and confinement in the Galaxy have not 
been well borne out, and it has become imperative to look for other possible 
scattering mechanisms as well. In this search a process which seems to 
offer considerable scope for the future is the collective effects of cosmic 
rays during their propagation in the Galaxy. For a recent survey on this 
subject we refer the reader to an interesting article by Wentzel (1974) . 

In interstellar space, the relativistic cosmic ray plasma, streams 
along the general large scale magnetic field lines; in this process they 
collectively excite a spectrum of hydromagnetic waves . These are Alfven 
waves which propagate in the ionized component of the interstellar gas along 
the streaming direction of the particles. Under the condition when the 
cosmic ray performs one gyration around the general magnetic field B while 
traversing one Alfven wavelength X, there is resonant interaction between 
the particles and these waves. The condition for this resonant self scat- 
tering (or cyclotron resonant scattering) can be written as X = r^ 
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where is the velocity of the particle along B, Also, the scattering 

mean free path X can be written as 
s 

^ 2 energy in the ambient magnetic field ^ ^ 

^s tT energy in the magnetic waves 

These relations highlight the following very interesting f eatures i (a) 
cosmic rays of different energies resonate with different waves; (b) since 
at each energy range cosmic rays make their own waves, off which they 
scatter, this self scattering process becomes increasingly ineffective at 
higher energies because of the falling power law nature of the primary 
spectrum; and (c) at a typical energy of about 1 GeV where the flux of 
particles, and hence the wave growth, is large, the scattering mean free 
path is exceedingly small and is only ''^10*"^ pc. In fact, it is the last 
feature which makes this process more attractive from the point of view of 
the diffusion of cosmic rays in the Galaxy. 

When there is a systematic general streaming or drift of cosmic rays 
relative to the backgroxmd interstellar plasma and if the drift velocity 
becomes comparable to the Alfven speed, the waves excited by cosmic rays 
scatter them and thereby inhibit the streaming (Wentzel, 1969; Kulsrud and 
Pearce, 1969). The damping of the Alfven waves in the Disk is due to fric- 
tion between the ionized hydrogen set in motion by the waves and the neutral 
hydrogen. From these ideas Wentzel (1974) has shown that the mean bulk 
velocity <v> is related to the Alven speed v^ = B/(4Trn^)^ and the number of 
cosmic rays n (>p) above a momentum p as follows! 

<v> = ^ )“■“ Nj^^g/N(>p) cm s"l 3.5 

Here T is the temperature of the medium and n^ and n^^ are the ionized and 
neutral densitites in the medittm, and 6 is the differential spectral index 
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of the cosmic ray spectnm. In the case of fully ionized medium, the 
damping of hydromagnetic waves is due to wave-wave interactions and as a 
result, the energy dependence of the bulk velocity (second term of Equation 3,5) 
becomes If one uses typical values for the parameters applicable 

to Interstellar space, one finds that particles of about 1 GeV are well 
scattered, those at about 10 GeV are weakly scattered and those above 100 
GeV are hardly scattered at all. 

Thus, in spite of the many attractive features of the mechanism of 
self-scattering, it is disappointing that in its present form it is lanable 
to explain the high degree of the observed isotropy of cosmic rays of 
energy 10^^ - 10^^ eV, 

3.2.3 Instability model 

In this model, wherein pioneering work has been carried out by Parker 
(1969), it is assumed that the cosmic ray gas diffuses freely along the gen- 
eral galactic field lines, which are parallel to the plane of the Galaxy, 
except in regions of small scale magnetic irregularities. In fact it has 
been shown recently (Earl, 1974) that if coherent propagation of cosmic 
rays is effective in interstellar space, they can be transported almost 
scatter free along the field lines. Diffusion perpendicular to the plane 
occurs either due to the ergodic nature of the field lines or the stochastic 
property of the field. The cosmic ray gas escapes from the Disk by in- 
flating the outer regions of the field. 

The presence of the supratherraal cosmic ray gas in the large scale 
magnetic field of the Galaxy does not affect the usual hydrodynamic waves 
even though the thermal plasma and the suprathermal gas are tied to the 
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magnetic lines of force; the physical reason for this is that the suprathermal 
gas is 80 mobile, except when self-scattering becomes important, that it can 
redistribute itself with ease along the lines of force, thereby avoiding 
compression and rarefaction associated with the hydromagnetic modes. Though 
the suprathermal wave propagation has also some effect on the thermal gas, 
the role of the thermal gas is to provide enough Inertia to make the field 
act as a rigid gxiide along which the suprathermal mode propagates. Thus, 
we find the intimate connection between the gas-field system and the propa- 
gation of cosmic rays. We shall not qualitatively examine the manner in 
which diffusion of cosmic rays occurs across the galactic plane. 

Firstly, on account of the ergodic nature of the galactic field along 
which the suprathermal gas streams, it gradually tends towards isotropy and 
uniform density distribution. The characteristic time in which the cosmic 
ray gas relaxes towards this uniform distribution is of the order of the 
dimension of the region of space divided by a few times the Alfven speed 
in the ionized component of the gas. The thermal gas, on the other hand, 
permits the field to expand and decrease its energy while streaming along 
the lines of force. The cosmic ray gas would then inflate the outer por- 
tions of the field where the cosmic ray pressure becomes larger than the 
magnetic pressure (P > B^/8tt) and therefore tends to escape. 

Secondly, the stochastic nature of the field which may prove even more 
effective, is a direct consequence of the random turbulence of the gas in 
which the field is embedded (Jokipii and Parker, 1968; 1969a). It has been 
shown (Jokipii and Parker, 1969b) that any given line of force random walks 
from the equatorial plane of the Galaxy to its surface, while traversing 
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distances of (5-10) x 10^ pc. Thus the cosmic rays are transported across 
the Disk by the very stochastic nature of the magnetic lines of force. If 
this is true, then the cosmic rays do not have to travel far to enter regions 
in which B^/8 ti < P; inflation would then set in regions at about every 10^ 
pc along each tube of magnetic flux. The instability resulting from the 
inflation is opposed by the gas perturbations and by the tension in the 
magnetic field, while the cosmic ray pressure increases the inflation. 

Such a situation may lead to a state of turbulence, and cosmic rays may 
escape out of the Galaxy either freely or probably carrying with them part 
of the magnetic lines of force. 

3, 3 The diffusion process and isotropy of cosmic rays 

It has been seen in Section 3.2 that in the various models described 
so far, cosmic rays stream along and across the interstellar magnetic field 
lines and diffuse in all directions. However, if the streaming of particles 
is not the same in all directions, it would result in an anisotropy. De- 
tection of such an anisotropy, in a region with dimensions of the order of 
the particle^s radius of curvature, would give information on the net 
streaming velocity of cosmic rays in any particular direction and thus lead 
to an insight into their diffusion and propagation. Thus far, observations 
have yielded only upper-limits to the magnitude of the anisotropy. We 
shall now try to imderstand the implications of these observational results 
on the basis of the diffusion of cosmic rays in interstellar space. 

The cosmic ray anisotropy is defined as 6 = (I - I , )/(I + . 

max min' ^ max ^ min' * 

where I is the Intensity in a given direction. From a study of the side- 
real variation of underground muon intensity, it has been shown by Elliot 
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et al, (1970) that at energies ^ 10^ GeV, the observed anisotropy is x 
10“*^. It has been pointed out by Parker (1969) that at these energies the 
irregularities in the interplanetary magnetic field itself could signifi- 
cantly isotropise cosmic rays because the gyroradius of a proton at 10^ GeV 
is only 0.4 AU in a typical magnetic field strength of 5 x lO”^ gauss. In- 
deed, it has been shown by trajectory calculations (Speller et al. 1972) 
that at rigidities less than a few hundred GV, the particle's arrival 
directions are well scattered in the solar system. Therefore, taking into 
account the response of their detector as a function of primary energy, 

Speller et al. (1972) have concluded that in the energy region 10^ - 10^ 

GeV, 6 < 1.5 X 10"*J. Since the anisotropy can be related to the streaming 
velocity of cosmic rays as v c6/(2+3), one finds that for the observed 
value of 6, the corresponding streaming velocity is ;^10 kms"^. Such a low 
anisotropy would require for particle escape from the Disk a lifetime of 
^3 X 10^ yrs in conflict with the estimates of '\'10^ yrs deduced for galactic 
disk models (Parker, 1969). 

The low streaming velocity of cosmic rays deduced from the low anisotropy 
observed by Elliot et al. (1970) have stimulated many workers to seek plau- 
sible explanations. It was first proposed by Jones (1970a) that the dis- 
crepancy in the cosmic ray life time for escape referred above could be due 
to fluctuations resulting from the statistical nature of the cosmic ray 
source distribution, which at the present happens by chance to give an ex- 
tremely low anisotropy at the Earth. Soon after, Ramaty et al. (1970a) pro- 
posed the statistical discrete- source model of local cosmic rays in which 
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they started with the acceptance that the general streaming velocity is 
Indeed as low as that suggested by the low observed anisotropy. These 
authors further assumed that cosmic ray sources are confined to the Disk 
and are supernova explosions, each treated as a random point source in 
space and time; they then characterized cosmic ray propagation in inter- 
stellar space by three dimensional isotropic diffusion with a constant mean 
free path X. From such an approach these authors noted that in order to 
account for the observed anisotropy, one needs X < 10“^ pc; furthermore to 
accommodate the observed ratio of the flux of light to medium nuclei at 
relativistic energies, one needs 10”^ pc A ^ 10“^ pc. Also, in this 
model the life time for escape from the Disk is ^10^ yrs and hence the life 
times of various nuclei are determined by their nuclear interactions and 
not by escape. These authors further pointed out that the discrepancy be- 
tween isotropy and the otherwise needed escape life time of 'vlO^ , can- 
not be accounted solely by fluctuations as suggested by Jones (1970a) unless 
the propagation is strictly one dimensional. Following this, Jones (1971) 
presented an argument in favor of an essentially one dimensional propagation 
of cosmic rays and showed that the distribution of streaming velocity is 
strongly peaked towards zero; thus the probability of seeing low observed 
anisotropy is more likely than any other value. Nevertheless, in this 
model one has to accept that the observed low anisotropy is only local to 
us and is due to an accident of our position in space and time, 

A serious difficulty in the statistical discrete-source model of 
Ramaty et al. (1970a) is the small scale size of ;^0.1 pc for the magnetic 
irregularities in interstellar space. Additionally, a proper understanding 
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of the observed abundance of positrons, deuterons and secondary light 
nuclei was also posing problems. In order to overcome these difficulties, 
Lingenfelter et al. (1971) modified their earlier model by considering the 
combined effect of one dimensional diffusion of cosmic rays along the magne- 
tic field lines and to the three dimensional random walk of these field 
lines, a treatment first discussed by Getmantsev (1963). Thus, the resul- 
tant motion of a particle is due to a compound diffusion with two mean free 
paths Xp and X^ to describe the scattering of particles along the field 
lines, and the random walk of the field lines, respectively. The quanti- 
tative treatment of the compound diffusion model of Lingenfelter et al, 
(1971) expressed in a simplified fashion by Allan (1972) will now be briefly 
described. As a result of the random walk of the field lines, the total 
length of the line of force along which the cosmic rays propagate from 

the source to the Earth can be written as L = d^/x where d is the dis- 

m m 

tance between the source and the Earth. Since particles have to traverse 


this distance by a process of one-dimensional diffusion, we can, by a 

similar argument, write for the total particle path length L = L ^/X = 

P m p 

d^/(X^X ). In consequence, the transmit time is proportional to d^ in this 
m p 

case instead of d^ for simple diffusion. Furthermore, the increase in 

transit time, compared with straight line motion, is by a factor d^/X ^X ; 

m p 

hence the outward streaming is reduced correspondingly and the mean aniso- 
tropy can be written now as 


6 

av 


2X /d3 
P 


3.6 


Setting Ip employing a value of 'v 10^ pc corresponding to 

a few g cm“^ of matter traversed by cosmic rays in interstellar space having 
a mean density of '^1 atom cm“^, one deduces a value of X 20 pc for 5 3 x 

10 “^ . 
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It is further argued by Allan (1972) that in the compound diffusion 
model each magnetic tube of force is required to have a minimum width much 
in excess of the distance coursed by the solar system with a velocity of 
*^20 km with respect to the galactic frame over the period of a few 
years when the measurements were in progress. Consequently, the observa- 
tional data corresponds to the microscopic anisotropy within a single tube 
of force 6 appropriate to the simple one dimensional diffusion, which is 

3 

written as 

6 ^ X /X = \ X /d^ 3.7 

s pm pm 

If this were so, then d 1 pc and X 10“^ pc leading us back to the same 
difficulty of too short a mean free path. 

Thus in summary, if one accepts seriously the low value of anisotropy 
experimentally determined so far by only one group (Elliot et al. 1970) , 
we are yet to propose an acceptable and internally consistent physical 
process which will explain it in a natural sequence. Therefore, it goes 
without saying that such a. consequential observation should be measured 
by other experimenters with greater accuracy and reliability in the years 
to come. And finally, one is left wondering that perhaps in all these 
attempts a fundamental concept in the propagation and confinement of cosmic 
rays in the Galaxy is still missing. 



40 


3.4 General Formulation of the Propagation of Cosmic Rays In Interstellar 
Space 

So far we endeavoured to understand the mode of propagation of. cosmic 
rays In interstellar space in relation to their dynamical properties. The 
obvioijs in^ression one gets from the existing theories is that there is 
perhaps no unique manner in which the diffusion of cosmic rays in the Galaxy 
can be treated quantitatively. Furthermore, during propagation in inter- 
stellar space, the composition and energy spectrum of cosmic rays are modi- 
fied due to nuclear interactions and other energy loss processes; incorpor- 
ating these in the diffiision theory is likely to introduce further difficul- 
ties, Nonetheless, it is essential to formulate quantitatively the propa- 
gation of cosmic rays in a simple form in order to make a beginning to 
understand and interpret observations. Modifications and improvements could 
then be incorporated, or new formulations attempted, depending on the degree 
of success of the theoretical predictions to explain observations. Towards 
this we shall first proceed to develop a general formulation for the propa- 
gation of cosmic rays, and then discuss various simplifying approximations 
which lead to different models of propagation; in this formulation the dy- 
namical aspects of cosmic ray diffusion will not be taken into account. 

The general equation of propagation, which allows for spatial diffusion, 
energy losses, acceleration and transformation of chemical composition in 
interstellar space can be written in the form (Ginzburg and Syrovatskii, 

1964) . 


V 

3t - 




3.8 
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th 

where is the density of particles of the i kind defined as 

Nj^(E,r,t) / (E.r.t.k) d 3.9 

Here, I^(E,_r,t,k) is the intensity in the direction _k, and d^ is the 
element of solid angle in k-spacej if the intensity is isotropic, then 
N = 4irl/v, We shall now examine Equation 3.8 in detail. The first term 
on the left hand side of Equation 3.8: This term represents the rate of change 

of density of cosmic rays per unit time due to processes described by the 
other terms in the equation. The second term on the left hand side of Equation 
3.8: This term is the rate of net flow of particles in the co-ordinate 
space, where (E,_r,t) is the net flow of particles defined by the mode 
of diffusion; it can be written in general as 

J = -D V N + v*N 3.10 

where ^ is the diffusion tensor and v' is the velocity of the scattering 
centers. Under isotropic diffusion approximation. 


J = -D V N 3.11 

where D is the diffusion coefficient. 

The third term on the left hand side of Equation 3.8: This includes the 

continuous energy loss by particles due to processes, such as ionization 
and synchrotron and inverse-Compton radiation and continuous energy gain 
from possible acceleration mechanisms as in the case of Fermi process. The 
net change in the energy per unit time can be written as 


f = 


3.12 


One can also include in principle an additional term - %9^/9E^(e.N^) 
which can take care of fluctuations in the continuous variation in energy. 
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In this case the term e^ =* ~ <AE^> is the mean square of the 

energy increment per unit time. However, we shall not introduce such a 
term in our considerations since, firstly, the amount of fluctuation expected 
in the energy loss by cosmic rays during their propagation in the Galaxy is 
expected to be small, and secondly. It is believed that the statistical 
acceleration mechanism does not play any important role in the propagation. 

The fourth term on the left hand side of Equation 3.8: This term corresponds 
to catastrophic energy losses, such as the emission of high energy brems- 
strahlung quanta by an electron and spallation of heavy nuclei resulting in 
the loss of particle identity; here p^ is the probability per unit time for 
such losses . 

The term on the right hand side of Equation 3.8: This represents the pro- 

duction rate of cosmic rays and can be written as 

Q^(E,r,t) = q^(E,r,t) + 1 \(E\r,t) dE’ 3.13 

k 

Here, the first term on the right represents the production rate from sources 
situated at a distance r^ at time t, while the second describes the secon- 
dary particle production as a result of catastrophic collisions. p^j^(E',E) 

th 

is the probability of production per unit time of the i type of particle 

th 

with energy E from catastrophic collisions of particles of the k type 
having energy E'. 

The general solution of Equation 3.8 can be obtained by making use of Greenes 
function which satisfies the following equation (Syrovatskil , 1959), 

-D7G + Is- (bG) -I- pG = «(r-r') 5(E-E’)6(t-t') 


3.1A 
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The solution of Equation 3.8 can then be written as 

N^(E,r,t) = /// dr’ / dE’ / dt’ Q^(E' ,r, t)G(E,r,t; E’,r*,t’) 3.15 

Here, the Green's fimction G is called the propagation function since Equation 

3.14 represents the propagation part of the general diffusion equation. 

For infinite space, the propagation function has the form 

G(E,r,t; E’ ,r’ ,t’) = — exp [-(r-r’ ) ^/4X-t ' /T]6 (t-t’ ,x’ ) 3.16 

b(E)(4irX)3/2 

where t' is the time taken for the particle of energy E’ at to reach £ 
with an energy E due to continuous energy loss processes and can be written 
as 

(E.E-) = t - f = ^ 3.17 


Also 


X 



D(E) 

b(E) 


dE 


3.18 


which becomes D t’ , when the diffusion is independent of particle energy; 

T = 1/p is the mean life time for catastrophic loss processes, where p is 
the probability of such processes per unit time. 

For quasi-steady state condition, where the explicit time derivatives 
in Equation 3.8 vanish, we get 

N,(E,r) = /// dr’ f ,^D -. exp[-(r-r’)2/4X-x’/Tl 

^ - " |b(E)|(4^X)3/2 

. - /// dr’ / — exp[-(r-r’)2/4DT’-r’/T] 3.19 

|b(E) 1 (4 ttDt') 


In the special case of an unbounded confinement space with a homogeneous 
distribution of sources therein, the above solution reduces to 
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\ 


N^(E) 


b(E) 


/g. exp 


1 .E dE” 
" T ^E’ b(E”) 


Q^(E»)dE* 


3.20 


The limits of the integral over E* is determined from the behavior of the 
energy loss parameter b(E) with energy. If T depends on the energy of the 
particle. Equation 3.20 can be written, without loss of generality, in the form 

Ni(E) - /g. exp I - dE"/{b(E") T(E")}] q^(E’) dE> 3.21 

This equation reveals that the density of particles does not depend upon 
the diffusion mechanism but only on the continuous and catastrophic energy 
loss processes. Such a situation is also applicable in space with finite 
boundary conditions wherein the diffusion of particles is very slow. 

The exact solution of Equation 3.8 in space with finite boundary and an 
arbitrary distribution of sources, is quite difficult to obtain even under 
steady state conditions. However, one can evaluate Equation 3.19 numerically for 
any arbitrary source distribution by in copo rating spatial dependence of the 
energy loss processes through Equation 3.17. In case of homogeneous distribution 
of sources both in space and time, it is customary to assume that the diffu- 
sion is very fast within the finite volume and that particles are lost by 


slow leakage from the boundary. In such a situation, the diffusion term in 
Equation 3.8 is replaced by a catastrophic loss term with an equivalent leakage 
life time t which is the residence time of cosmic rays in the finite volume 
of confinement; t can either be energy dependent or independent. Under this 

situation. Equation 3.21 would take the form 

1 E dE" 1 1 

" I (b(Eyi "^E’ ^ “"^E* b(E") ^ **' t (E") ^ 


This solution is widely used by researches in the field and it is called 
the steady state solution for cosmic ray propagation. 
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We considered so far the general formulation of the diffusion problem 
in the Galajqr. The next step is to apply it to a situation involving the 
various observations within the framework of models so far proposed. These 
are reviewed in what follows in Section 3,5, 

3.5 Quasi-Equilibrium Models 

Almost all interpretative studies on the chemical composition and energy 
spectrum of cosmic rays are based on propagation models under quasi-equilibrium 
conditions. For the sake of clarity in presentation, we shall separately 
consider the electron conqjonent because of the dissimilar nature of dominant 
energy loss processes associated with their propagation in space, 

3.5,1 Rucleonic Components 

In case of nucleonic components, at all energies except perhaps ^10^® eV, 
the only important continuous energy loss process is ionization; this loss 
becomes increasingly severe at sub-relativistic energies, where one may also 
have to take into account elastic collisions particularly with nuclei of low 
mass. The catastrophic losses are primarily due to inelastic collisions and 
leakage from the confinement space. In inelastic collisions, heavier nuclei 
suffer fragmentation and lose their identity, giving rise to nuclei of lower 
charges; the associated change in the velocity of the fragment is inappre- 
ciable except at sub-relativistic energies. The fragmentation parameters 
involved are, in general, energy dependent at low energies and seem to saturate 
at the relativistic domain. For protons, inelastic collisions at energies 
above the threshold for meson production is of relevance thoiigh its contri- 
bution in the general propagation can be disregarded except in special cir- 
cumstances where there is slow diffusion in a confinement region with high 
matter density, such as source regions. 
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Two major approaches have been adopted so far to study the propagation 
of cosmic ray nuclei in space. In the first, only interaction losses are 
taken into account by assum i ng that all particles arrive at a given point 
in space after the traversal of a fixed distance. An example of this is the 
"Slab Model" in which it is assumed that particles of a given energy tra- 
verse a fixed length of trajectory and hence the same amount of matter. In 
the second, leakage is also folded in by considering the confinement of cos- 
mic rays in a given region of space; the "Steady State Model" in which all 
possible vacuum trajectories are allowed with a finite probability belongs 
to this category. 


(a) The Slab Model 

The simplified diffusion equation which includes both ionization loss 
and fragmentation can be written as (Aizu et al. 1960; Ray, 1960; Fichtel 
and Reames , 1966) . 


Here, j defines the differential intensity of particles, E the energy per 
nucleon, w(E) = dE/dx is the rate of energy loss per nucleon per unit path 
length, x the path length in g cm”^ , the attenuation mean free path in 
g cm”^ and the fragmentation parameter defined earlier. f\_^ can be de- 
fined as X./(l-p..) where X. is the interaction mean free path and P.^ the 
fragmentation parameter relating to the probability of the particle retaining 
its identity after interaction; this relation also holds when i represents 
a group of nuclei but does not apply for protons in which case P^^ has to 
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be replaced by the inelasticity parameter. The solution of Equation 3.23 under 
such an approximation would be 


j^(E,x) = 


w(E^) 

w(E) 


j^(E’,o) exp [-/ 


dx’ ^x 

A,(E',x») i 


Vi (E * .y) 


w^(E) 


f" 


Z 

k>i 


jk(E",y) 
X J.(E") 



(E”) 


dy exp[-/* - 


dz 

^(E*,z) 


1 


3.24 


The first term here represents the single particle propagation, while the 
second corresponds to the propagation of the i type of nuclei produced 
during the traversal of nuclei heavier than the i^^ type, or simply, the 
secondary particle propagation; j^(E* ,0) is the intensity of particles at 
the source with energy E' such that ionization loss reduces it to E at the 
point of observation. Since X is a function of energy, its energy dependence 
is taken care of at each depth. In the second term, E* corresponds to the 
energy of the i*^ type of nuclei during production at y due to the corres- 
ponding k^ type of nuclei with the relevant energy E", The integration is 

carried out numerically from the heaviest particle onwards so that at each 

tli 

depth the flux of the k type would be available for the propagation of the 
lighter nuclei i < k. 

Even under the slab approximation, certain variations in the method of 
computation have been introduced leading to very different results. Two 
such variations may be mentioned here, (i) A distribution in the path 
length about the mean amount of matter traversed <x> by the particles, has 
been considered by Balasubrahmanyan et al. (1965). In this case the flux 
of particle can be written as 


j^(E) = /“ j^(E,x) p(x) dx 


3.25 
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where j^(E,x) is defined by Equation 3.23 and p(x)dx is the normalized distribution 

CO 

such that / p(x)dx = 1, From a consideration of a cosmic ray source situated 
o 

at a distance r^ from the observer, in a homogeneous medium with isotropic 
particle diffusion, it is found that p(x)dx is approximately Gaussian and 
the distribution of path length is of the form 

p(x)dx = i;8858<x 7 ^ V] 3.26 


This model necessarily predicts a spatial gradient of the intensity of cos- 
mic rays by virtue of its assumption that the source is at a finite distance 
from the observer and that diffusion takes place Isotropically in an infi- 
nite medium, (ii) In the second, the matter traversed by particles is 
assumed to be velocity dependent. This can be achieved by replacing the 
limit X in Equation 3.23 by x0, where 6 is the ratio of the velocity of the par- 
ticle to that of light (Durgaprasad, 1967) . This would imply that all par- 
ticles at the point of observation have spent the same amount of time in 
space, leading thereby to a concept that cosmic rays were produced by a sin- 
gle event at a time corresponding to t = x/pc, where p is the mean matter 
density in the medium. Such a hypothesis would result in temporal varia- 
tions and cannot be then brought under quasi-equilibrium models. 

(b) The Steady State Model 


The steady state model assumes space- time homogenity of cosmic ray 
sources within the confinement volume. Under this assumption, particles 


have finite probability of having finite path lengths before undergoing 
catastrophic loss either from collision or leakage from the boundary. As 
discussed earlier, Equation 3,21 represents the steady state solution for cosmic 
ray propagation. Transforming now the energy variable to matter traversed. 


we can rewrite the equation as 
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}] 3.27 


.th 


i ' ' 4tt o iwusj 1 - - ' ‘ ■ o 

where Q’CE'jX) is the x-lndependent production rate of the i'"“ particle de- 

_2 

fined by Equation 3.13, per gram per unit path length traversed at depth x g cm 
having energy E’ such that it reduces to E due to ionization at the point of 
observation; X is the mean matter traversed by cosmic rays corresponding to 
the mean residence time t. Once again, one can introduce here variations 
in the method of computation by making plausible assumptions regarding the 
energy dependence of X (Cowsik, et al. (1967; Fichtel and Reames , 1968; 

Ramaty and Lingenfelter, 1968; Shapiro, et al. 1970a). It is also pertinent 
to mention here that attempts have been made recently to analyze the observed 
data using the diffusion model with boundary conditions (Ramaty and Lingen- 
felter, 1971; Pacheco, 1971). 

3.5.2 Electron Component 

In case of the electron component, almost all energy loss processes are 
applicable. At energies less than about a few tens of MeV, ionization loss 
becomes increasingly important while at very high energies inverse Compton 
and synchrotron processes dominate. At intermediate energies escape loss 
is perhaps a more effective catastrophic energy loss than bremsstrahlung. 

Since the rate of energy loss due to inverse Compton scattering and synchro- 
tron radiation are proportional to the square of the electron energy, these 
processes lead to a heavy drain on the energy of energetic electrons. In 
consequence, observations on the shape of the electron energy spectrum could 
lead to information on the mode of propagation of cosmic rays and the re- 
gion of space where they, and hence all cosmic rays, are confined. 
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a) Leakage life-time model 

The simplest model, so far extensively used in all interpretative studies 
of the electron spectrum is to apply the leakage life time concept; accor- 
dingly Equation 3.8 can be written as 

II = Q - [-(a-«E + bE2)N] - [ i + i 1 N = 0 3.28 

Here, the positive numerical co-efficients a and b are related to the rate 
of energy loss through ionization and, synchrotron and inverse Compton pro- 
cesses respectively and a is related to the energy gain due to Fermi type 
of acceleration. The catastrophic losses are defined by their life times 
T and t, for bremsstrahlung and leakage respectively. The life time against 
bremsstrahlmg is energy independent and in the model of galactic confine- 
ment of cosmic rays T » t. The solution to Eqmtion. 3«28 can be written from 
Equation 3.21 as 

N(E)dE i L Q(E')dE' exp [- | 1 3.29 

(a-aE + bE^) ^ ^ ^ (a-aE" + bE"^) 

The upper limit to this integral is either E^, the energy at which the rate 

of energy loss equals the rate of energy gain, or <» depending upon whether 

E is less or greater than E respectively. It is obvious from the above solu- 

s 

tioh that for - 4ab - 0, N(E)dE at E = (a±/oi^-4^b) /2b . At these 

3 

energies, the rate of energy loss is equal to the energy gain and Equation 3,28 
has to be written as 

is. = Q - ( - + a-2b E)N = Q- i - + /a‘^-4"^") N = 0 3.30 

St T s T ^ 

and the corresponding solution at E = E^ is 
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N(E)dE 


Q(E) dE 


3.31 


— X /a^-4ab 

T ^ 

The singularity in this solution occurs when the physical conditions are 

such that at energy E^, the catastrophic loss of particles is same as the 

particle gain due to either energy loss or gain processes. Moreover, since 

the rate of energy loss is zero at E , one cannot expect steady state to 

s 

prevail at this energy and one can see from Equation 3.30 that 9N/3t = Q and 


hence N(E)dE ast-»-“atE = E . One can also reason from the above 

s 

discussion, without going through any mathematical treatment, that if the 
acceleration of cosmic rays is a continuous process like that of Fermi 
acceleration, then there will be an enhancement in the energy spectrum at 
the energy where the rate of energy loss is the same as the rate of energy 
gain. 

Having discussed the effect of acceleration in steady state propagation 
of cosmic ray electrons, we shall pass on to a particular case in which 

D 

a = 0 and the production spectrum is a power law of the type Q(E)dE = kE”^dE, 
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Under these conditions, one can see from Equation 3,28 that for very low 

energies, where ionization loss alone is effective, the equilibrium 
spectrum is of the form. 

N(E)dE = dE for E « a(3-l) 3.32a 

Thus for constant energy loss, the spectrum is flatter by one power. 

In the energy region, where the particle loss is dominated by the leakage 
from the confinement volume, the equilibrium spectrum retains its spec- 
tral shape, provided the leakage is independent of energy, and from Equation 
3.28 one obtains 

N(E)dE = kTE’^dE for a(3-l) « E « 3.32b 

At high energies, where the energy loss due to synchrotron and inverse 
Compton process dominates » the equilibriiim spectrim becomes 

^ 3.32c 

Here the spectrum is steeper by one power than the production spectrum. 

Thus one finds that the spectral shape changes at two characteristic 
energies and if such changes are detected, it would give information 
on the processes responsible for the same, 
b) Diffusion Model 

Attempts have been made to determine the solution of Equation 3,8 for an 
arbitrary source distribution but by neglecting the catastrophic loss 
term p = 0 and taking only the synchrotron and inverse Compton loss 
processes b(E) = -bE^ (Syr ovatskii , 1959 ; Jokipii and Meyer, 1968; Berkey 
and Shen, 1969; Webster, 1970; Jones, 1970; Bulanov et al, 1972). For 
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simple isotropic diffusion Equation 3.8 then reduces to 

II - 1^ (b£2) - DV^N = Q 3.33 

~6 

When the source function Q is expressed as Q = kE 6(tQ-t06(rQ“X*) * 
which implies a source at a distant position r^ injected an electron 

_Q 

spectrum of the type kE at a time t^, the non-stationary solution at 
a position r and at time t of Equation 3.33 can be written (Berkey and Shen, 
1969) as 

N(E,r,t) = exp [- ) 1 E~®[l-bE(t-t^)]®"^ 

o'^ 


for E < 


b(t-t^) 


for E> 


b(t-t^) 


Here, the first factor corresponds to the spacial dilution of particle from 
a point source, the second to the probability of the emitted particles to 


reach the point at r and the last factor gives the depreciation of the par- 


ticles in an energy interval due to energy loss processes. Obviously 
the spectrum would show a continuous steepening as a function of energy 
and drop to zero at an energy E ^ l/b(t-t^) beyond which the particles 
cannot survive due to the energy loss processes during a time (t-t^) even 
if they had infinite energy to start with. If the source emits 


particles continuously, the steady state solution can be obtained by 

integrating Equation 3,34 over t^ between the limits t - 1/bE and t. Even 

here, the equilibrium spectrum terminates at energies >l/bE. 

—6 

For a source function of the type Q(r\ E) = g(r^)E where g(r*) 
is arbitrary source distribution, the steady state solution of Equation 3,33 
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is, according to Webster (1970) 

N(E.r) - /// 4^D|v-f' T 

where f(a,B) = r(a-l)/»^ exp (-$)$^U(a-l, 3/2,3) and U is the confluent 
hypergeometric function of the second kind. In this equation f ^ 0, 
when E >> 4D/b|r-r*|^ and f 1 when E « 4D/b|r-r*|^. One can in 
principle evaluate numerically Equation 3.35 for any arbitrary source dis- 
tribution with boundary conditions. We shall now take a typical case 
of a flat disk (like the Galaxy) of radius R and height 2H, with uniform 
source distribution g(r) = const and determine the equilibrium spectral 


shape at the centre of the disk r = 0. 

For energies E < 4D/bR^, where the energy loss is so small that the 
particle can diffuse through the disk with very little or no modification 
to its energy, the function f in Equation 3.35 tends to one for all values of 
r* and hence the resultant spectral shape is the same as that of the 
source spectrum. As the energy increases , electrons from all parts of 
the disk cannot reach the observer at r = 0 and fewer sources contribute 
to the observed flux resulting in the steepening of the spectrum. Since 
the sources, which contribute to the flux at r = 0 are uniformly distri- 
buted in a thin disk, we have from Equation 3.35 

.2 


N(E 


. 0 ) 


f(3. ^ ) 2.rdr 


4irDr 


4D 


a E 


-(3+0.5) 


^ 4D r. 
for ■ ■ - n - < E < — ' 

hR^ ^ uti- 


3.36 


Here we find that the spectrum steepens by half power at an energy 
E ^ 4D/bR^, and this spectral shape continues up to an energy E ^ 4D/bH^ 
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beyond which the sources contributing to the flux are distributed within 
a spherical volume of radius H. As the energy increases even this 
spherical volume shrinks in size and from Equation 3.35 we obtain 


at ^ 

“ / hoi ^ ^ 


-( 6 + 1 . 0 ) 


for E > 


Thus at energies E > 4D/bH , the equilibrium spectral shape is steeper by 
one power compared to the production spectrum 

One can notice in the above treatment that for an axially symmetric 
source distribution of the type g(r)dr a rdr, the resultant spectral 
shape is steeper by half power, while for a spherically symmetric source 
distribution, g(r)<i£ « r^dr, it is steeper by one power. These clearly 
indicate that the spectral shape of the equilibrium spectrum is determined 
by the type of source distribution. Hence Webster (1970) showed that 
any arbitrary break of A 6 in the spectrum can be obtained by a source dis- 
tribution of the type g(r)dr a r^^^dr. 

We have so far assumed that diffusion is isotropic and in an infinite 
medium, but introduced different source distributions with boundary con- 
ditions. It is also possible to evaluate Equation 3.33 for a finite diffusive 
volume with a transmission co-efficient depending upon the physical con- 
ditions, which determine the boundary, and for anisotropic diffusion 
(Berkey and Shen, 1969; Bulanov et al. 1972). The uncertainties in the 


present experimental results do not warrant at this stage a detailed 
description of these complex models. However, it is easy to show (Berkey 
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and Shen, 1969) that for energies » 4D/bR^, the boundary at |r| = R 

plays no role because the particle can never reach the boundary without losing much 

of its energy; on the other hand at energies « 4D/bR^, the equilibrium 

spectrum depends upon the boundary conditions. It is important to point 

out ab this stage that none of the models predict an abrupt spectral 

Steepening, and the range of energy over which the spectral change occurs 

^ about a decade or more In energy depending upon the boundary conditions 

(Webster, 1970; Jones, 1970b). 
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4, Interpretation of the observational data 

Cosmic rays we sample in the neighborhood of the Earth reach 

us after having undergone a chain of transformations beginning with 

what occurs in their source regions. Often one refers to the composi- 
of 

tion/cosmic rays at the time of injection into interstellar space; on 
such occasions one has implicitly included modifications introduced 
during the process of acceleration and propagation, if any, within the 
source region. The interstellar Injection spectrum, thus defined, under- 
goes further changes while the particles propagate in galactic space and 
attains a state of quasi-equilibrium. On entering the solar system this 
equilibrium spectrum is further modulated at energies below a few GeV/n, 
when the particles have to work their way against the continuous outflow 
of solar wind from the Sun. The last phase of the journey of cosmic rays 
before being detected by instruments carried by balloons, rockets, and 
earth satellites, is their entry into the region of influence of the geo- 
magnetic field. However, there is no real difficulty in calculating 
quantitatively the influence of the geomagnetic field on cosmic rays 
entering the magnetosphere since our knowledge of the geomagnetic field 
and its effect on the motion of charged particles is now adequately under- 
stood. This knowledge then permits us to infer reliably the cosmic ray 
spectrum beyond the magnetosphere, i.e., in the interplanetary space, 
from observations made within; instruments on deep space probes, on the 
other hand, reveal the nature of the interplanetary cosmic rays directly. 
From the above brief accoimt, it is seen that the process of deducing 
information on the propagation of cosmic rays in the Galaxy from earth 
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based observations is a very complex one, subject to many uncertainties. 

In spite of this, considerable progress has been made in this field during 
recent years; and in what follows we shall try to deduce information on 
the propagation of cosmic rays by applying theoretical approaches de- 
scribed in the preceding section on the observational data. 

4. 1 Matter traversed by cosmic rays 

During their propagation in the Galaxy, cosmic rays interact with 
interstellar matter. When interactions are of electromagnetic type, the 
particles lose part of their energy by ionization. Nevertheless, there 
occur less frequent interactions which are catastrophic in nature; in 
these the following may take place: 

(i) heavy nuclei may break up into lighter ones; 

(ii) at energies greater than the threshold for meson production, all 
created unstable particles will decay in space and give rise to 
stable components like electrons and photons; 

(iii) the target and/or the colliding particle may be left in an excited 
or radioactive state, which will then decay to a stable one by the 
emission of particles and electromagnetic radiations; and 

(iv) in case of the electron component it can suffer bremsstrahlung 
radiation leading to gamna ray emission. 

Though in principle all these processes can lead to information on the 
amount of matter traversed by cosmic rays, thus far only a few could be 
profitably exploited; these will be reviewed now. 

4.1.1 The effect of ionization 

The rate of energy loss due to ionization of a particle, with charge 
Ze and velocity 8c, in space according to Ginzburg and Syrovatskii (1964) , 
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is 

^ = 1-524 X 10 ^**z£n ^ n^2 + ln(-^) -3^} MeV/sec 4,1 

dt 3 1-3^ 

where n is the number of hydrogen atoms cm"^. The above equation shows 
that the energy loss due to ionization decreases with the increasing 
velocity, reaching a near constant value of about 4Z^ MeV per g.cm ^ of 
hydrogen at relativistic energies. Since at subrelativistic energies the 
predominant energy drain for the nucleonic components is through ioniza- 
tion losses, one can in principle deduce their matter traversal by first 
starting with an injection energy spectrum and then comparing the spec- 
trum modified by ionization losses, for an assumed traversal of inter- 
stellar matter, with the observed spectrum. However, in practice, this 
is found to be a questionable procedure since, at non- relativistic ener- 
gies , the resultant spectral shape is sensitive to the propagation model 
employed, the form of the injection spectrum assmed, and the degree of 
solar modulation. It has been demonstrated quantitatively (Comstock, 

1969) that so far no single model is capable of reproducing the relative 
spectral shapes of the observed energy spectra of all components with a 

given spectral shape at injection and with the same parameters for solar 

> 

modulation. In the high energy region, GeV per nucleon, where the 
above difficulties don^t exist, the ionization loss is too small to have 
any observable effect on the energy spectrum. Notwithstanding all this, 
it has been recognized that a prior knowledge of the amount of matter 
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traversed from other considerations can be used to obtain information on 
the mode of propagation of cosmic rays at low energies. 

4.1.2 Fragmentation of nuclei 

It stands to reason to consider that the elements and isotopes, 
known to be nearly absent in astrophysical objects, are also absent in 
cosmic ray source material. Hence these nuclei, if observed in the 
cosmic rays reaching the solar system, can only stem from fragmentation 
of heavier nuclei in collisions with matter encountered enroute (Bradt 
and Peters, 1950). The most striking examples of elements and isotopes, 
which are believed to be nearly absent universally and which have been 
extensively used in the study of matter traversed by cosmic rays, are 
Li, Be, and B, and %e and respectively. With continuously improv- 
ing capability and sophistication of instrumental techniques, other rare 
elements and isotopes are also being increasingly studied during recent 
years. For the reason that in this method one cannot distinguish 

between matter traversed in the source region and in interstellar space, 
we shall assume in the following general treatment that the matter 
traversed by cosmic rays is in interstellar space. 

Before proceeding further, it is necessary to emphasize that the 
mean amount of matter traversed in this method also depends upon the 
propagation model employed because each model assumes different path 
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length distributions. This can be easily demonstrated by considering 
two extreme models, namely "the slab model" which assumes a 6-function 
for the distribution of path lengths, and "the steady state model" which 
assumes an exponential distribution. For simplicity we shall consider 
at relativistic energies the transformation of primary nuclei of the 

i^^ kind to secondary nuclei of the kind (Z^ > Z^) without taking into 
account the contribution from other nuclei either to k or i. Then according 
to slab model, one gets from Equation 3.23, 

j^(x) = j^(0) exp (-x/A^ 4.2 

and 

P., A. A. 

j|^(x) * j^(0) {exp(-xM^) - exp(-x/A^)} 4.3 

i k i 

where x is the constant amount of matter traversed by the primary cosmic 
ray (obviously the secondaries traverse less matter); the other quantities 
are defined in the same manner as in section 3.4. The ratio R(sl) of the 
fluxes of secondary to primary nuclei in slab model would be 


R(sl) 



A. A, 

1 k 

A^-A^ 


{exp [x(“ ~)] -1} 

i k 


4.4 


On the other hand in the case of steady state models, one obtains from 


Equation 3.26 


^i 


A^+X 



A, A^X^ 
k 1 

(Aj^+X) (A^+X) 



Qi 


and 
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A 4X X 
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ik 
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4.5 


4.6 


Here is the production rate per sr. sec. g.cm ^ and X is the mean matter 
traversed by the particles if they are not attenuated by nuclear interactions. 
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For a given observed value of R, the matter traversed as estimated by 
these two models can be related as 


A,(X+AJ 

JS i /JS L-\ 

Vi 


.(X+A^) 


4.7 


and for A, -A , the above equation reduces to x = XA, /(X+A.) , which shows 
k i K X 

that X is always less than X and the difference (X— x) increases with the 

decreasing value of A (i.e. , with increasing value of charge). When A^ 

and are very much larger than x. Equations 4,4 and 4.6 become identical and 


X = X = 


R(ob) 


4.8 


ik 


The above relations can be easily generalized to take into account all 
fragmentation leading to the i^^ and k nuclei and also for groups of 
nuclei . 

Keeping in mind the meaning of the mean amount of matter traversed 
by cosmic rays, let us interpret the observations. The observed ratio of 
the integral flux of light nuclei L (Li, Be, B) to that of medium nuclei M 
(C,N,0) at energies ^1.5 GeV/n is 0.25±0.02 (von Rosenvinge et al. , 1969c); 
the value extrapolated to the outside of the magnetosphere is 0.23±0.02 
(Shapiro and Silberberg, 1970). One can then use appropriate cross-sections 
for the production of L nuclei from nuclei of (Silberberg and Tsao, 

1973a, b) and by a method of successive comparison of the calculated ratio 
with the observed one, deduce the amount of mean matter traversed as about 
4 g.cm"^ for the slab model and as about 6 g.cm”^ for the steady state 
model. 

Likewise, X and X can also be deduced from the observed ratios of 
%/He and %e/(%e+%e) on the assumption that % and %e are absent in the 
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source. Note that the production of these isotopes are primarily from 
collisions of **He nuclei with interstellar gas and that the related attenu- 
ation mean free paths are much larger than the amount of matter traversed. 

For these reasons, an accurate determination of these ratios above a GeV/n 
would give us a more dependable value for the mean amount of matter trav- 
ersed by cosmic rays although the value thus derived will again be model 
dependent. Unfortunately, reliable finite values for these ratios are not 
available so far (Ganguli et al. , 1967; Damle, 1968; Apparao, 1973) and the 
observations lead only to upper limits of 'vlO g.cm ^ to the matter traversed. 
4.1.3 Charge composition of cosmic ray electrons 

Like in the case of universally rare nuclei, the positron component 
of the cosmic radiation is also expected to result from secondary processes 
occurring in interstellar space (also in the source regions) . Since the 
attenuation mean free path of electrons is very large (for example, the 
radiation length In hydrogen is 58 g.cm”^) compared to the mean matter 
traversed by cosmic rays, it can be seen that in the energy region, where 
continuous energy losses through ionization, inverse Compton scattering, 
and synchrotron radiation do not dominate, the electron component is ideally 
suited for the study of matter traversed by cosmic rays. The mean matter 
traversed, which in this case is model independent, can be written as 

X . A (E) dE f„t(E^.E2) j (E) j^(E)dE A.9 

El 

In this equation j^(E) is the differential flux of the total equilibrium 
electrons in interstellar space. One may recall here that the energy spec^ 
trum of equilibrium electrons is the same as the radio emitting electrons 
in the Galaxy (Section 5,3) and is not very different from that observed 
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in the neighborhood of the Earth at energies *^>2 GeV. The quantity (£^^,£ 2 ) 

is the observed fraction of positrons within the energy interval and £ 2 ^ 

Similarly, j (E) is the differential flux of equilibrixim secondary elec- 
s 

trons per g.cm"”^ of matter traversed by cosmic rays (Section 4,5) at energy 

E; f (E) is the corresponding fraction of positrons. 

3 

Equation 4.9 essentially relates the total flux of positrons in the 
interstellar space to its production per g.cm”^ of matter traversed by 
cosmic rays at energies of interest here. The main contribution of posi- 
trons to primary cosmic radiation comes from the decay 7T**‘-^y‘^->e‘*' . The 

charge composition of secondary electrons as a function of energy can be 
calculated from our knowledge of the characteristics of high energy inter- 
actions. In Figure -4.1 is shoim the charge ratio n^+/n^_ = f/(l-fy-aa a 
fimctlon of energy from two such recent calculations (Ramaty and Lingen- 
felter, 1966a; Perola et al. , 1967) above 50 MeV, where this ratio is 
insensitive to the propagation models. One notices that large differences 
exist between these calculations; these arise obviously from our incom- 
plete knowledge of the energy spectrom of and K-mesons in the labora- 

tory system over all angles of emission as a function of primary energy. 
However, using the existing information and the observed fraction f^^ 

(1.7 GeV, 4.2 GeV) = 0.046 ± 018 (Fanselow et al., 1969), one obtains 
values of 3.5 and 2.0 g.cm"^ for x corresponding to Curves I and II in 
Figure 4. r respectively. It needs to be pointed out here that, since the 
observed charge ratio for high energy atmospheric y-mesons is -1.2, a 
value close to the one according to Curve I in Fig. 4.1 (Wolfendale, 1969), 
the most probable value for the mean matter traversed by cosmic rays is 
likely to be about 3 g.'cm"^. It is hoped that more reliable determinations 
in the future of the positron fliox in the cosmic rays as a fimction of 
energy will lead to a better understanding of this subject. 




Fig. 4.2 ; The calculated ratio of positrons to electrons arising from nuclear collisions of 
cosmic rays with interstellar gas as a function of energy. 



66 


4*1.4 Gamma ray component of the galactic background radiation 

Of the many processes which lead to the emission of background gamma 
rays of galactic origin through interaction of cosmic rays with interstellar 
matter, there are two which, from a practical point, are capable of 
yielding information on the gas density in the Galaxy. These are (i) 
bremsstrahlung of cosmic ray electronsj and (ii) decay of neutral pions 
produced in nuclear interactions of cosmic ray nuclei with interstellar 
matter. It will be shown later (Section 5.1) that the contribution of 
bremsstrahlung radiation is about an order of magnitude smaller than that 
of plonic gamma rays and, from a comparison of the observed flux of gamma 
rays at energies ^100 MeV towards the galactic centre, one can set an 
upper limit of about 2 hydrogen atoms cm“3 as a mean matter density in the 
galactic disk. 

4.1.5 Summary 

At non-relativistic energies, where the ionization loss and solar 
modulation play. an important role in modifying the spectral shape of 
cosmic rays, it is difficult to obtain reliable information on the mean 
amount of matter traversed by cosmic rays from the observational data. 

As yet only two methods have been used with a good degree of success at 
relativistic energies. In the first the ratio of light to medium nuclei 
observed at the top of the atmosphere yields a value of about 4-6 g cm”^ 

The second method, wherein the flux of the positron component is employed, 
leads to a value close to 3 g cm” ^in this the mean amount of matter tra- 
versed is not sensitive to the propagation model used. From these obser- 
vations one may conclude that relativistic cosmic rays traverse a mean 


\ 

\ 
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amount of 3 to 6 g of matter in the Galaxy. It needs to be pointed 

out that we have so far defined the matter traversed in terms of the 
amount of hydrogen in g cm"^; when consideration is given to the composi- 
tion of the interstellar matter it is expected that it will effectively 
increase the amount of matter traversed as deduced above (Apparao and 
Ramadurai, 1967). Recent observations also suggest (Garcia-Munoz , 

1973) that the mean matter traversed by cosmic rays is probably an energy 
dependent parameter; this aspect will be discussed separately in section 
4.4. 

4.2 The chemical composition at the source 

The chemical composition of the nucleonic component of cosmic rays 
at the time of their injection into the interstellar space would give 
information on the elemental abundances (the source characteristics Including 
nucleosynthesis), acceleration processes and injection mechanisms operative 
in the source regions. It is possible to deduce the source composition 
from the observed elemental abundances by making use of the available 
fragmentation cross-sections of cosmic ray nuclei and plausible models of 
propagation in interstellar space. Many attempts have been made in the 
past to determine source composition, either for groups of nuclei or for 
prominent elements using the slab model for the propagation of cosmic rays 
(Hayakawa et al. , 1958; Aizu et al. , 1960; Badhwar et al, , 1962; Kristiansson, 
1966; Beck and Yiou, 1968; Waddington, 1969); more recently it has been 
carried out by using different propagation models (Shapiro et al, , 1970a 
1973; Ramaty and Ligenfelter, 1971; Webber et al. , 1972; Casse and Goret, 

1973; Cowslk and Wilson, 1973). In all these calculations, the fragmentation 
cross-sections used are mostly based on semi-empirical relations. We 
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shall first enumerate some of the difficulties associated with this kind of 
work so as to make the reader aware of the uncertainties involved , before 
we proceed to understand the recent studies on this aspect of cosmic ray 
propagation. 

The measurements of fragmentation cross-sections as a function of 
energy are still Incomplete, and theoretical estimates do not cover ade- 
quately the reactions of interest. Semi-empirical relations based pri- 
marily on the work of Rudstam (1955, 1966) and a few modifications 
introduced recently (Audouze et al. , 1967; Shapiro and Silberberg, 1970; 
Silberberg and Tsao, 1973n,b) are widely used for the purpose of computa- 
tions. These relations constructed by using measured cross-sections avail- 
able from time to time are being constantly revised. The typical standard 
deviation of the calculated cross-sections from a measured value (where the 
latter is available), according to Silberberg and Tsao (1973a ,b) , varies 
from about 30% for target elements of Z - 20 to about a factor of two or 
more for heavy elements. Further, most of these estimates refer to inter- 
actions with hydrogen only, though helium constitutes about 10% of the 
atoms of interstellar gas. New observations from accelerators show that 
in a carbon target, the cross-section for the production of isotopes like 
^Be and ^^Be with a-particle beams are a few times larger than the corre*- 
spending ones with proton beams (Fontes , et al. , 1971) thereby pointing to 
the need to measure in the laboratory the fragmentation cross-sections of 
cosmic ray nuclei with helium. 

The determination of the source composition of cosmic rays even at 
relativistic energies depends crucially on the type of propagation model 
employed. For example, let us consider the source abundance of iron and 



69 


oxygen nuclei. Using first the slab model, one can write for their source 
ratio according to Equation 4,2; 

R(Fe/0)g^ = exp (x/a^^ - x/aq) [ jp^(x) /jQ<x) ] 4.10 

In the case of the steady state model, one gets from Equation 4.5 
R(Fe/0)g^ = [(Ap^+X)/(Aq+X)1 (A^^/Aq) [ jp^(X) /J qCX) ] 4.11 

Here the quantity [j^^Cx) /j^Cx) ] = [jp^(X) /j(j(X) ] is the observed ratio of 
the flux of iron nuclei to that of oxygen nuclei andA andA^ are the 
attenuation mean free paths having values about 2.5 g cm~^ and 8 g cm"”^ 
respectively; x and X are the mean amount of matter traversed by cosmic 
rays as defined by the slab and steady state models respectively. One 
finds from these equations that for the same observed ratio, the slab 
model always predicts a higher ratio of iron to oxygen nuclei at source 
than the steady state model. Further, as x the ratio becomes infinite 

for the slab model while it reaches a constant value in the case of the 
steady state model. Taking now the ratio between Equations 4,10 and 4*11, one 
gets 

R(sl/st) = exp(x/Ap^-x/AQ)/[(Apg+X)Ap^/AQ(AQ+X)] 4.12 

which is independent of the observed ratio. Using x = 4 g cm”^ and X = 6 

thereby 

g cm”^ respectively for the slab and the steady state models ,/ insuring that the 
two models predict the same ratio of light to medium nuclei as observed 
(Section 4.1.2)pne gets R(sl/st) Z 2. From the above discussions, it 
appears that at present, the main source of error in the determination of 
source composition comes from the uncertainties in the propagation models 
and fragmentation cross-sections employed. 

In Table 4.1, we show the cosmic ray source abundances relative to 
carbon as calculated by Shapiro et al, , (1973), who made use of the 
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Table 4.1 

Elemental Abundance Relative to Carbon in the Source Region 


Element 

Abundance 

Element 

Abundance 

Z 

Abundance 

H 

5 X 105 

Si 

200130 

40-44 

6.213*3 X 10-3 

He 

26000 

S 

3016 

50-54 

1.010*3 " 

C* 

1000 

Ar 

715 

55-59 

8.212.9 X 10“*+ 

N 

110±20 

Ca 

2218 

75-79 

9.711.5 " 

0 

1070120 

Cr 

313 

80-84 

4. 1+1.0 ” 

Ne 

160+20 

Fe 

210130 

85-90 

2.911.2 " 

Na 

8+4 

Ni 

812 

90-94 

4.911.2 " 

Mg 

230+20 

Cu-Se 

0.6110,22 

Al 

414 X 10"5 

A1 

20±10 

Br-Y 

.011.004 




* Normalization 
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steady state model with X = 5 g These abundances are in agreement with 

those calculated independently by other researchers using the same model 
(Ramaty and Lingenfelter, 1971; Webber et al. , 1972; Casse and Goret, 

1973). The abtindances given in this Table 4.1 assume a power law spectrum 
in energy per nucleon, but if the spectrum obeys a power law in rigidity, 
then the value for hydrogen would be 2 x 10^ instead of 5 x lO'*. It has 
been observed that there is a striking similarity of the abundances of 
elements in solar cosmic rays to those in the photosphere of the Sun, 
thereby implying that, whatever it may be, the process of acceleration 
has no charge dependent preference (Biswas and Fichtel, 1965) . If this 
argument is extended to galactic cosmic rays, then their source composi- 
tion derived above can be related to the features of the chemical composi- 
tion of sources which generate cosmic rays. On this account, a comparison 
with the solar abundances can yield some interesting results on the origin 
of these elements. 

The ratio (CRS/SS) of the abundances of cosmic ray source elements 
to those of solar system relative to iron is plotted in Figure 4,2 r; as a 
function of atomic number (Shapiro et al. , 1973). In this plot, the solar 
abundances for neon and sulphur are taken from the measurements of solar 
flare particles (Berts ch> et al. , 1973; Mogro-Campero and Simpson, 19 72a, b; 
Crawford et al. , 1972; Teegarden et al. , 1973b). This figure reveals that 
while cosmic ray sources seem to be under-abundant in elements lighter 
than neon, there is a striking similarity with solar abundances above neon. 
This prompted Silberberg et al. (1973) and Kozlovsky and Ramaty (1973) to 
suggest that all elements in cosmic ray sources are produced through 
explosive nucleosynthesis (Arnett, 1969; Truran and Arnett, 1970; Arnett 
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and Clayton, 1970) and the enrichment of solar system abundances of 
elements lifter than neon can come from other processes such as hydro- 
static helium burning (Clayton, 1968). The under-abundance of hydrogen 
and helium in cosmic ray sources resembles that expected for materials 
ejected from supemovae (Truran and Cameron, 1971). 

Contrary to the above point of view, it has been argued by many 
that cosmic rays are not directly associated with violent events. It 
was suggested by Bradt and Peters (1950) that a selective electromagnetic 
acceleration process is necessary to account for the high abundance of 
heavy nuclei in cosmic rays; one can therefore examine whether this is 
related to ionization properties. In Figurfe/i.3 are pio-t-ted.^the^^raiio 
(CRS/SS) as a function of the first ionization potential of each element 
up to nickel; this plot exhibits a possible correlation that the rela- 
tive cosmic ray abundances decreases as the first ionization potential 
increases (Kristlansson , 1971, 1972; Havnes, 1971, 1973; Casse and Goret, 
1973). Havnes (19733 proposed that rotating magnetic stars could be the 
possible sources of cosmic rays, which can accelerate interstellar matter 
with a selective acceleration dependent upon the ionized state of the 
interstellar gas around the star. Cowsik (1971) suggested a preferen- 
tial injection mechanism for heavy elements into the accelerating regions 
d^ to the low equilibrium charge state of heavy elements compared to 
their rigidity (Korchak and Syrovatskii, 1958), the observed increase of 
the iron to helium ratio in large solar flares at small energies (Price 
et al. , 1971a) has also been attributed to this mechanism. Nevertheless, 
because of the large errors in the ratios (CRS/SS) , it is difficult at 
present to distinguish between explosive nucleosynthesis with further 
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acceleration and possible preferential acceleration mechanisms. 

The study of abundances of cosmic ray nuclei of charge >30 gives 
a clue to the process of nucleosynthesis. The abundance of ultra heavy 
nuclei at source given in Table 4,1 is based on the estimates of Fowler 
(1973) who made use of all the existing observations (Price et al. , 1971b; 
O’Sullivan et al. , 1971; Shirk et al. , 1973; Blanford et al. , 1973a,b; 
Binns et al,, 1973; Fowler et al, , 1973), These abundances seem to be 
in better agreement with r-process than with the total solar abundances 
relative to iron (Fowler et al. , 1970; Price et al. , 1971b); however, a 
more reliable distinction can be made only when individual abundances 
are measured, because in some of the charge groups used in Table 4.1 only 
a few elements , depending upon the process of nucleosynthesis, are ex- 
pected to dominate the observed flux (Seeger et al. , 1965). It is also 
found that the ratio Pb/U at the top of the atmosphere is about 1.3, 
indicating thereby that the UH nuclei are young. 

Before closing this section, let us briefly examine the possibility 
of obtaining some information on the propagation models from the source 
composition. As we have mentioned earlier that the slab model would 
predict a source abundance for iron a factor of two larger than the steady 
state model; at the same time, the slab model cannot explain the observed 
abundance of ultra heavy nuclei (Fowler et al. , 1970; Cowsik et al. , 

1970) . Ramaty and Lingenfelter (1971) have made an extensive study of 
the various propagation models such as simple 3-dimensional diffusion 
with boundary conditions, compound diffijsion and an infinite medium with 
exponential escape time. In doing this, these authors find that the 
observed abundances of fragmentation products from lithium to iron can 
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be understood provided the following conditions are satisfied in the 
respective models: (i) for simple diffusion with absorbing boundaries, 

the scattering mean free path Xg cannot be larger than 0.1 pc; (11) for 
simple diffusion with reflecting boundaries, the amount of fragmentation 
depends principally on the escape probability at the boundary and is 
nearly Independent of the scattering mean free path with the escape life 
time t^ 1,5 x 10^ yrs; (iii) for compound diffusion, the characteristic 
scale length in the 3-dimensional random field ~ ^ 30 (Section 

lU p 

3.3) with tg O' 2 X 10^ yrs; and (iv) for an exponential path length the 
value of tg 'v 3 X 10^ yrs corresponding to X - 5 g cm“^. Obviously, there 
is no clear justification to choose or to reject any one of these models 
on the basis of elemental abundances alone when there is a wide choice 
of physical parameters that can be varied. 

4.3 Modification of the energy spectrum at low energies 

It is found that the energy spectrum of cosmic ray nuclei observed 
in the terrestrial vicinity can be described by a simple power law at 
energies greater than a few GeV/n. Below this energy, the spectrum 
deviates from a simple power law such that the differential intensity 
first increases slower with decreasing energy, reaches a maximum at 
around a few hundred MeV/n, and then decreases at still lower energies; 
the spectral shape below a few GeV/n is not the same for all components. 
This change in the spectrum can result from the combined effect of solar 
modulation, energy loss in interstellar space during propagation, and 
perhaps to some extent from differences in the spectra at injection. In 
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this section we consider only the effect of cosmic ray propagation in 
interstellar space for different models using various injection spectra, 
and compare the results with observations. 

It can be shown from Jlqua.tion 3>ii9 tiiair. in- ener^. T^ (a) 

continuous energy loss processes are either negligible or proportional 
to energy and (b) catastrophic loss processes are independent of energy, 
the shape of the energy spectrum of cosmic rays, when expressed as a 
power law in kinetic energy per nucleon, does not change during the 
propagation. This condition is satisfied by the nucleonic components at 
relativistic energies where the ionization loss is negligible, and the 
nuclear interaction cross-sections and escape from the Galaxy are reason- 
ably independent of energy. At these energies experimental observation 
cannot distinguish between power law spectra either in rigidity j(R)«R“6 
or in kinetic energy per nucleon j(E)«E“^ or in total energy per nucleon 
j(V0®W”^. Since the source spectrinn at low energies is rather difficult 
to deduce directly from observations, one is tempted to assume that the 
spectral shape at high energies continues right down to the lowest ener- 
gies. In Figure 4 ar^^sbownf; -the . power spectra in kinetic energy-: per 
nucleon (Curve E) , rigidity (Curve R) , and total energy per nucleon 
(Curve W) plotted in the kinetic energy per nucleon scale for the same 
constant of proportionality. One can see from this figure that while 
these three forms of spectra become indistinguishable at high energies, 
they differ significantly at low energies. Currently, there is no 
apriori reason to prefer any one of these spectra over the others because 
each one has its own claim to be chosen: (i) the rigidity spectrum may 



GeV) 


- 78 - 



KINETIC ENERGY IN GeV 


Flfi, 4.4 ; The proton intensity as a function of kinetic energy, plotted 
in arbitrary units, for three different types of power law spectra. 
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result from acceleration and storage by the magnetic fields in the source 
region (Kaplon and Skadron, 1966); (ii) the shockwave acceleration in 
supemovae can lead to a power law in kinetic energy per nucleon (Col- 
gate, 1968); and (iii) the Fermi-type of acceleration (Fermi, 19A9) can 
lead to a power law in total energy per nucleon. We shall now attempt 
to bring out the merits and demerits of these three source spectra as 
we analyze the observed spectra of nuclei. 

The effect of ionization on the energy spectrvnn of cosmic rays in 
space can be well exemplified by starting with a power law spectrum in 
kinetic energy at the time of particle injection into interstellar space. 
The results are shown in Figure 4.5 for the spectral change of the differ- 
ent nuclear components by making use of the steady state model with X = 

6 g (solid curves) and the slab model with x = 4 g cm“^ (dotted 

curves). Curve I in this figure is the injection spectrum; Curves p, 0, 
and Fe are typically the equilibrium spectra of hydrogen, oxygen, and 
iron nuclei, which are normalized to the given injection spectrum around 
3 GeV/n; the Y-scale in this figure is expressed in arbitrary units. One 
can see from this figure that the effect of interstellar ionization losses 
is to reduce the flux at low energies; in the case of the steady state 
model, this implies that low energy particles come from nearby sources. 

It also becomes clear from this figure that the slab model suppresses 
the flux of particles of low energy more than the steady state model and 
this effect becomes more pronoimced as the particle charge increases. 
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KINETIC ENERGY IN GeV/n 

Fig. 4.5 : The effect of ionization during propagation is shown using 

a steady state model with X = 6 g cm"^ (solid curves) and slab model with 
X = 4 g crn’’^ (dotted curves) for a power law spectrum at Injection (Curve 
I); Curves P, 0 and Fe are for protons, oxygen and iron nuclei respectively. 
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In Figure (4.6 are shown, the equillbritCQ spectra o.f hydrogen and 
helium nuclei in Interstellar space for three different types of in- 
jection spectra. In getting these results, the steady state model 
has been used with X = 6 g cm"^. Curves E, R and W correspond to the 
power law Injection spectra in kinetic energy per nucleon, rigidity 
and total energy per nucleon respectively, which are normalized at 5 
GeV with the observed spectra (curves 0) during the period of minimum 
solar modulation.lt is apparent from this figure that if the equilibrium 
spectrum in kinetic energy exists below 100 MeV/n, in the same manner 
as indicated by this figure, the cosmic ray energy density in inter- 
stellar space would continue to increase as the energy decreases; 
however, this behavior is not exhibited by other spectra. The attempts 
made so far to deduce information on the spectral shape of cosmic rays 
at injection using different modulation theories yielded one conclusive 
result?iame]y f;he source spectrum is closer to a power law in total 
energy per nucleon (Meyer, 1971; Comstock et al. , 1972; Ramadurai and 
Biswas, 1972) or a rigidity spectrum which flattens below about 1,5 GV 
(Burger, 1971), However, none of these analyses could reproduce the 

entire spectral shape of the observed protons and helium nuclei right 

be 

from about 10 MeV/n to a few GeV/n; this can, perhaps, /due to our in- 
complete knowledge of the solar modulation and/or intrinsic 
nature of source spectra these nuclei have at the time of injection. 

nevertheless, solar modulation does not significantly affect the ratio 
fluxes of 

of /elements as a function of energy ; in what follows briefly 

review the work carried out so far in this direction. 
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^ « 6 . Equilibrium spectra of protons and bellum nuclei In Interstellar space ^ using steady 
state model with X = 6 g cnT^ for different injection spectra shown in Figure 4,5 (Curves E, R and W) , 
are compared with the observed spectra during the period of minimum solar modulation (Curve 0), 
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4.3.1 L/M Ratio 

Since the light nuclei L (Li, Be, B) in cosmic radiation are produced 
mainly by the fragmentation of medium nuclei M (C, N, 0), their ratio of 
Intensities as a function of energy has been extensively made use of to 
understand the propagation of cosmic rays in the Galaxy (Badhwar and Daniel, 
1963; Biswas, et al. , 1966; Cowsik, et al. , 1967; Fichtel and Reams, 1968; 
von Rosenvinge, et al. , 1969b; Gloeckler and Jokipii, 1970; Syrovatskii 
and Kuzhevsky, 1970; Garcia-Munoz and Simpson, 1970; Shapiro, et al. , 

1970b; Pacheco, 1971). In Figure 4,7 , t we shown. the' observed. ratio. :1L/M 
as a function of energy (Anand, et al. , 1969b; Freier, et al, , 1966; 

Hagge, et al, , 1968; Malmquist, 1967; Garcia-Mxmoz and Simpson, 1970; 

Reams and Fichtel, 1968; von Rosenvinge, et al. , 1969b; Webber and Ormes, 
1967; Mason, 1972) along with the theoretical estimates. In this figure. 
Curve A is the calculated ratio using steady state model with X = 6 
g cm”^ for a power law spectrum in rigidity at injection (Tandon, 1970) 
and Curve B is for a power law spectrum in total energy per nucleon 
(Mason, 1972). Notice that while Curve A is clearly in disagreement with 
the observations at low energies , Curve B too fails to give adequate agree- 
ment with the data points; it can also be shown that the calculated ratio 
using slab model as well as Gaussian distribution for path lengths would 
be close to Curve B (Ramadurai, 1970). It is thxis seen that a single 
power law in rigidity as a source spectrum is incompatible with the obser- 
vations and that either there could be an interference of a second component 
at energies below a few hundred MeV/n, or the adiabatic deceleration in 
the solar wind becomes important below about 200 MeV/n, 
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Fig, ^ , *7 ; The ratio of light to medium nuclei as a function of kinetic energy per nucleon. Curves 
A and B are the calculated ratios in interstellar space using steady state model for power law spectra in 
rigidity and total energy per nucleon respectively, at injection. 
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4.3.2 Ratio 

This ratio is insensitive to solar modulation since the charge to mass 
ratio of and **He is the same, and the observed ratio as a function of 
energy should reveal some aspects of interstellar propagation of cosmic 
rays. nuclei are primarily produced during collisions of p + p + 

^He, and a + ^He. Many theoretical estimates of the ratio ^H/*^He are now 
available as a function of energy (Meyer, et al. , 1968; Ramaty and Lingen- 
felter, 1969; Hsleh and Simpson, 1969; Comstock, et al. , 1972; Ramadurai 
and Biswas , 1972) . In these calculations one may have to take into 
accotint at low energies the effect of elastic collisions toward the re- 
distribution of energy of the colliding and target particles (Ramadurai 
and Biswas, 1972). In Figure 4,8, we display.. the recent measuremests; of ;;the 
ratio ^H/*^He (Hurford , et al. , 1973; Teegarden, et al., 1973b) , along with 
those summarized by Hsieh, et al. (1971); the year of measurement is also 
shown against the authors. Curves A and B in this figure are the inter- 
stellar ratios estimated by Meyer (1971) using steady state model with 
X - 6.3 g cm~^ for a source spectrum of the type (E + where = 

500 MeV and 0 respectively; Curves C and D are for steady state model with 
X = 7 g cm for a source spectra of the type and respectively 

(Biswas. and ^Ramadurai, .1971) . Prior to the 1972 data, it was concluded by 
Simpson (1971) that the observations imply a power law in total energy per 
nucleon at injection. However it is evident from Figure 4, 8; that, the ratios 
measured during 1972 at energies below 50 MeV/n are lower than the earlier 
ones by about an order of magnitude. 
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■£■^■5.*. ratio of deutrons to heliijm nuclei as a function of kinetic 

energy per nucleon; theoretical estimates of this ratio in interstellar space 
shown by solid curves are all on the basis of steady state model for different 
Injection spectra as described in the text. 
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Buggttttlng thereby that either the source spectrum is close to a 
power law in kinetic energy per nucleon without much adiabatic deceler- 
ation, or there could be a significant enhancement in the interstellar 
'♦He at low energy, which is noticeable at 1 AU even with adiabatic 
deceleration and without a corresponding increase in ( Stone, 1973), 

4.3.3 Ratio 

Since the mass to charge ratios of these two nuclei are different, 

also 

their solar modulation wlll/be different. This aspect was first em- 
ployed by Biswas et al (1967) to determine the modulation parameters. 

3 4 

In Figure .4.9 the observed ratio of ..He/. He is shown as compiled by 

Hsieh and Simpson (1970) along with the calculated ratio 1 AU, 
after correcting for solar modulation, including adiabatic deceleration 
as a function of energy. Curve A in this figure is the calculated 
ratio using a Gaussian path length distribution (Section 3.5.1) with 
<x> - 4 g cm“^ and a Fermi spectrima (Biswas and Ramadurai, 1973); 
curve B is that using steady state model with X = 4 g cm“^ and a power 
law in total energy per nucleon (Comstock et al, 1972), One can see 
from this figure that it is difficult to distinguish between the two 
models of propagation employed in the calculations because of the 
effect of adiabatic deceleration below about 200 MeV/n. 

4.3.4 He/M Ratio 

In Figure 4yl0, the observed ratio of the helium to medixjm nuclei 
as a function of energy is shown (Balasubrahmanyan, et al, 1966; Corn- 

stock et al. , 1969; Fan, et al. , 1968; Hagge, et al. , 1968; Mason, 
Garcia- 

1972; yMunoz and Simpson, 1970; Teegarden et al, , 1970; von Rosenvinge, 
et al. , 1969c), In this figure curve A is the estimated ratio using 
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Fig, 4,9 : The ratio of helium isotopes as a function of kinetic energy 

per nucleon. Curves A and B are the estimated ratios at 1 AU (after correcting 
for adiabatic deceleration in the solar system), using Gaussian distribution 
and steady state model respectively. 
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Fig. 4.10 : The ratio of helium to medium nuclei as a function of kinetic energy per nucleon. 
Curves A and B are the expected ratios in interstellar space on the basis of Gaussian distribution 
and steady state model respectively; Curve C is the expected ratio at 1 AU from Curve A. 
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a Gaussian distribution of path lengths with <x> = 4 g cm“^ for a 
power law in total energy per nucleon (von Rosenvinge, et al. , 1969c) 
and curve B is the expected ratio on the basis of steady state model 
with X - 7 g cm”^ (Tandon, 1970) ; curve C is the modulated spectrum 
from curve A (Mason, 1972) using adiabatic deceleration below 200 MeV/n. 
The trend in the data points in this figure below 100 MeV/n indicate an 
enhancement of medium nuclei with respect to helium. 

4.3.5 He/LH Ratio 

It has been shown earlier (Figure 4.5) that the effect of ionization 
increases with increasing charge and as a result, the ratio of helium 
to heavy nuclei should be sensitive to propagation models at low ener- 
gies. In Figure 4.11, we. have shown the observed ratio of helium nuclei 
to light heavy nuclei LH (10 1 Z - 15) as a function of energy (Reams 
and Fichtel, 1967; Fan, et al. , 1968; Comstock, et al. , 1969; von 
Rosenvinge et al. , 1969c), along with the theoretical estimates; here 
curves A and B are as described in Section 4,3.4. One notices 
from this figure that it would be rather difficult to explain these 
observations on the basis of propagation model with Gaussian distribu- 
tion for path lengths (curve A) even after correcting for adiabatic 
deceleration. 

4.3.6 He/VH Ratio 

Figure 4tl2 higihll^ts the ratio of helium nuclei to very heavy 
nuclei VH (20 i Z 5 28) as a function of energy,, where the data points 
are from the works of Reams and Fichtel (1968), Comstock, et al, , 

(1969) and von Rosenvinge et al. , (1969c); between 60 and 90 MeV/n, 
the data shown is deduced from the recent work of Price, et al. , (1973). 
In this figure curves A and B are the calculated ratio as described in 
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The ratio of helium to very heavy nuclei as a function of kinetic energy per nucleon; 
Curves A and B are the expected ratios similar to those in Figure 4.10. 



section 4.3.4 and it becomes quite obvious here that the Gaussian path 
length distribution (curve A) is incompatible with the observations. 
4.3.7 Summary 

From the above studies, the following conclusions can be drawn 
without much ambiquity 

(i) The source spectrum for cosmic rays cannot be a simple power law 
either in rigidity or in kinetic energy per nucleon, over the 
entire energy range so far covered. It is seen that, while a 
simple power law in total energy per nucleon is quite consistent 
with all the observations described so far, one may have to ex- 
press the source spectrum in terms of a rigidity spectrum with a 
desired flattening below about 1.5 GV to match the above require- 
ments, so that the interstellar electron spectrum can also be 
well explained by the same spectrum (Section 4.5). 

(ii) The steady state model consistently explains almost all observa- 
tions better than the other models of propagation. 

(iii) It becomes quite evident from Figures 4^7 to 4,13 that the obr 
servations below about 100 MeV/n cannot be understood purely by 
propagation effects in interstellar space; a similar inference is 
also obtained from a study of C/0 and N/0 ratios (von Rosenvinge, 
et al. , 1969b; Bhatia, et al. , 1970), These anomalies can be 
accounted if the adiabatic deceleration in the solar neighborhood 
amounts' to about = 200 MeV/n (Goldstein et al. , 1970 ; Gleeson and 
Ufarch, 1971), though one cannot rule out the possibility of the 
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low energy cosmic rays having a different origin (Biswas, et al. , 
1966; Comstock, 1969; Fichtel and Reaaes , . 19f8) . 

(iv) The observed turn up of the spectrum below a few tens of MeV/n 
(Mogro-Campero, et al. , 1953; Price, et al. , 1973; Hovestadt, 
et al., 1973; McDonald, et al. , 1974) and the anomalous charge 
composition Hovestadt, al, , 1973; McDonald, et al, , 1974) can- 

not be vinderstood at present by interstellar propagation and 
solar modulation described above, without invoking a new compo- 
nent having a different chemical composition. 
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4,4 Modification of energy spectrum at high energies 

In the proceeding sections, we Implicitly assumed for simplicity in 
all the calculations, in order to understand the observational data 
described so far, that (i) the energy spectra of all components of cosmic 
rays are similar at injection and hence an energy independent source com- 
position, (ii) the path- length distribution is identical for all components 
and, (ill) the matter density is constant over the entire storage volume, 
where the observed cosmic rays propagate. The preceedlng analyses did 

indicate that perhaps, not all these assumptions are valid, though at low 

definitive 

energies it is difficult to arrive at an jr / ' conclusion due to the effect 
of energy loss processes in interstellar space and solar modulation. At 
high energies, where the above difficulties do not exist, it is possible 
to obtain information on these aspects directly from observations. Indeed, 
the recent experimental investigations using different experimental tech- 
niques (Juiiusson et al. , 1972; Juliusson and Meyer, 1973; Smith et al., 
1973a; Balasubrahmanyan and Ormes, 1973; Ormes and Balasubrahmanyan, 

1973; Brown et al, ,1973b; Webber et al. , 1973b) clearly indicate a change 
in the charge composition of cosmic rays at high energies. In this sec- 
tion we briefly present these experimental data and critically examine all 
the theoretical explanations put forth so far to understand these obser- 
vations. 

4.4.1 Observed change in the abundances of secondary nuclei 

Any change in the abundance of secondary nuclei produced in inter- 
stellar space with energy can give some information on the possible de- 
pendence of propagation parameters on energy. In Figure 4, ,13 arp plotted 
the ratio of intensities of Li, Be, B and N to C and 0, at the top of 
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Fig. 4.13 : The ratio of secondary nuclei to C + 0 as a function of kinetic energy per nucleon; 

Curves A and B are two possible dependents of this ratio, one dn total energy and the other on kinetic 
energy respectively. 
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Fig, 4.14 ; The ratio of secondary nuclei to iron group as a function of kinetic energy per nucleon; 
curve shown here is similar to that of Curve A in Figure 4.13. 
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the atmosphere, as a function of kinetic energy per nucleon as summarized 
by Garcia-Munoz (1973) (Smith, et al., 1973b; Juliusson, 1974; Webber, et 
al. , 1973c) ; included also by crosses are the results from Golden, et al. 
(1973). For comparison, we have also shown in this figure two types of 
energy dependence, one varying as (Curve A) and the other varying as 

£-0,3 (Curve B) above 5 GeV/n, From this figure one can notice that, ex- 
cept for the results of Golden, et al. (1973) , the data points above 5 
GeV/n indicate an energy dependence as strong as that of Curve B. In 
Figure 4. 14., the observed ratios of sub-iron group to iron^group of nuclei 
are shown as a function of energy as given by Garcia-Munoz (1973) . In 
this figure, different authors (Smith, et al. , 1973b; Juliusson, 1974; 
Webber, et al. , 1973c; Balasubrahmanyan, et al. , 1973) have used different 
criteria in defining the charge groups and the observed scatter of points 
does not permit any definitive conclusions to be drawn as in the case from 
Figure 4., 13,^. for comparison we have also- drawn a aW Direct ^ 

measurements of the spectral indices of nuclei are given in Table 4.2 
from the work of Smith, et al. (1973a), Balasubrahmanyan and Ormes (1973), 
and Juliusson (1974), who have made use of power law spectra in rigidity, 
kinetic energy per nucleon, and total energy per nucleon, respectively. 

Again one notices from this table that the secondary nuclei are steeper 
than their primaries and their indices increase by about 0.2 between about 
a few GeV/n and 'vlOO GeV/n. 

It has been pointed out recently by Webber, et al. (1973c) that a 
stronger dependence of the type is operative right down to 100 

MeV/n, in contradiction to their previously reported results sximmarized 
in Figures. 4. 7 to 4*12. However, Price, et al. (1973) have observed nc such 



Table 4.2 


Observed Spectral Indices of Cnsmlc Ray Nuclei 


Elements 

Berkeley* 
R>8.5 GV/C 

Elements 

GSFC** 
E>3.3 GeV/n 

Elements 

Chicago! 

W>1,5 GeV/n W>20 GeV/n 

P 

2.63±.08 

P 

2.751.03 

Li 

2.951.12 

2.91.8 

He 

2.47±.03 

He 

2.771.05 

Be 

3.091.14 

3.21.7 

Li 

2.67±.13 

Li 

2.281,15 

B 

2.951.07 

2.81.4 

Be 

2,66±.12 

Be 

2 . 61.2 

C 

2.651.02 

2.91.2 

B 

2.76±.08 

B 

2.761.13 

N 

2.741.03 

3.11.3 

C 

2.541.04 

C 

2.521.06 

0 

2 . 531.02 

2 . 61.1 

N 

2.721.09 

N 

2.731.11 

Ne 

2.571.03 

2.91.2 

0 

2.521.05 

0 

2.571.06 

Mg 

2.561.03 

2.71.2 

Li,Be,B&N 

2.721.07 

B&N 

2.771.08 

Si 

2.501.03 

2.21.2 

C&O 

2.531.03 

C&O 

2.561.04 

S 

2.61.1 

— 

9iZll4 

2.521.07 

10<Z<14 

2.441.07 

Ar 

2.51.1 

— 



15<Z<23 

2.1 1.2 

Ca 

2.71.1 

2.51.3 





Cr 

2.61.1 

— 



Fe 

2.01.14 

Fe+tfo 

2.391.04 

2,21.2 

* for power 

law spectra in rigidity (Smith et 

al,, 1973a) 




** for power 

law spectra in 

kinetic energy per 

nucleon (Ormes 

and Balasubrahmanyan, 1973; 

Ba lasub r ahmany 


and Ormes, 1973) 


t for power law spectra in total energy per nucleon (Juliusson, 1974) 
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variation from about 2 GeV/n down to 30 MeV/n. Thus one can infer from 
Fi^.ree ,4%13 and 4*14r>anil:Table. 4.2r that .the j^urrent -observations s^gest 
that the ratio of secondary nuclei to primary nuclei varies either as 
VT”“ right from the smallest energies or as only above '''5 GeV/n, with 
a Z0.2±.l. This observed dependence would directly lead to the energy 
dependence of matter traversed by cosmic rays and, perhaps, the possible 
energy dependent confinement of cosmic rays either in the source or in 
the storage volume. 

4,4.2 Energy dependence in the abundance of primary nuclei 

The study of the change of abundance of primary nuclei as a function 
of energy can be used to examine whether its consequences are internally 
consistent with the propagation effects derived from the observed energy 
dependence of the ratio of secondary to primary nuclei; it can also give 
some additional information on source spectra. In Figure 4.15 are plotted 
the data on the ratio (C+0) / (Fe-group) as a function of energy from the 
observations summarized by Garcia-Munoz (1973) , after making necessary 
corrections to define the iron group as Z > 24. Likewise, in Figure. 4.16, 
the intensity ratios of LH nuclei (10 i Z i 14) to the Fe-group are 
plotted as a function of energy (Atallah, et al., 1973; Balasubrahmanyan, 
et al. , 1973; Juliusson, 1974); shown in these figures by shaded regions 
are the expected ratios in the source regions from Table 4.1. One may 
notice that there is an indication for the ratios to decrease with increas- 
ing energy and attain values close to the source ratios at high energies; 
this trend is also indicated in Table 4.2, where there is an observable 
change in the spectral indices with charge. 
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Fig, 4,15 : The ratio of C + 0 to iron nuclei as a function of kinetic energy per nucleon; Curve 

A is the expected ratio on the basis of steady state model using an energy dependence of the type shown 
in Figure 4.14, 
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Fig. 4.16: The ratio of light heavy nuclei to iron group as ^ a^ function of kinetic energy per 

nucleon; 'curve A is the expected ratio similar to that in Figure 4.15. 
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It Is, however, realized that the present status on observational 
data is to some extent unclear. For example, the observed ratio of 
(C+0)/(Fe-group) by Balasubrahmanyan, et al. (1973) are smaller than those 
obtained by other authors and are lower than the source ratio above 10 GeV/n. 
In contrast to this. Golden, et al. (1973) found no such variation with 
energy and the ratio (LH/Fe- group ) . Furthermore, it has been pointed out 
by Juliusson (197A) that the C/O ratio decreases at energies greater than 
20 GeV/n; the observed decrease is from about 1.1 around a few GeV to 
about 0.7 at '^lOO GeV/n, while the source ratio is about 0.9. Badhwar and 
Osborn (1973) found a similar trend while studying the ratio of Integral 
fluxes of C and 0 at energies >23 6eV/n. If these findings are confirmed 
by future experiments, it will imply that the source spectra are not the 
same for all nuclei! 

4.4.3 Theoretical models to explain the change in composition 

If one accepts for the moment the correctness of the observational 
findings for a possible energy dependence of the charge composition of 
cosmic rays, described above, they may be interpreted in many ways. There 
is no experimental evidence so far to suggest that the interaction mean 
free paths of nucleonic components increase with energy, in the energy 
range of a few GeV/n to 'vlOO GeV/n; in consequence, many of the explana- 
tions proposed so far have questioned the basic assumptions made in the 
beginning of Section 4.4 on propagation studies. These theories can be 
generally distinguished on the basis of either energy dependent galactic 
confinement of cosmic rays or energy dependent confinement in special 
regions within the galactic space, such as sources or clouds. In this 
section we briefly review these ideas and examine them critically. 
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(i) Energy dependent galactic confinement of cosmic rays 

noticed 

We have already / in section 3.2.2 that if self scattering of 
cosmic rays is the dominant mode by which cosmic rays are scattered in 
interstellar space, then the escape of cosmic rays from the Galaxy is 
energy dependent, resulting in a decrease of path length in the Disk 
with increasing energy. By taking into account the distribution of 
ionized and neutral matter perpendicular to the plane, Holmes (1974) 
showed that the energy dependence of path length is effective above a 
few GeV/n and can be approximated as above 30 GeV/n, where a<0.38. 

Many of the explanations so far put forward are based on this principle 
(Juliusson et al, 1972; Smith et al. 1973a; Audouze and Cesarsky, 1973; 
Meneguzzi, 1973a; Ramaty et al. 1973; Webber et al. 1973b) . For steady 
state models, the equilibrium intensity of i^^ kind of nuclei can be 
written from Eqn. 3.27 as 


+ kfl 


[l/A. + 1/X(E)] 


4.13 


and hence, the ratio of secondary nuclei to primary nuclei X(E) /(X(E)+A ) 

o r o 

while that of two different primaries ci(X(E)+A 2 ) /(X(E)+Ai) . From 

the observed change in composition, different workers obtained simple 


power law energy dependence for X with a value of a varied from about 
principle , 

0.2 to 1.0 ( in/ one can introduce higher order empirical relations 

like the one used by Cesarsky and Audouze, 1973). From Figures 4,13 and 

basis we have shown 

4.14, we find that the value of a is close to 0.2; and on this /in Figures 


4.15 and 4.16 the calculated ratios of primary nuclei (C+0) /(Fe-group) 
and (lOiZllA) /(Fe-group) respectively (Curve A) for a = 0.2. It can be 
seen from these figures that the observed variation of the ratios of 
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primaries is more than, that predicted from the variation of the ratios of 
secondary to primary nuclei. This has prompted Cesarsky and Audouze (1973) 
to postulate an energy dependent source composition and Ramaty, et al. (1973) 
to seek for an alternate source to account for the difference. The 
present experimental observations , with large systematic uncertainties, do 
not demand such hypothesizing at this stage. 

The energy dependent confinement of cosmic rays in the galaxy lead 
to the following difficulties (Rengarajan, et al, , 1973) , Firstly, the 
very low life time implied in this model at high energies would lead to 
large anisotropies in contradiction to the observed near isotropy (Section 
3.3), Secondly, since the energy dependent propagation requires that the 
observed spectrum is the steepened part of the equilibrium spectrum of 
nuclei, the spectrum at injection will have to be flatter than that ob- 
served; thereby raising the energy requirement from cosmic ray sources. 
Thirdly, this energy dependence cannot continue indefinitely but should 
cease at some energy corresponding to a lifetime of less than yrs . , 

the time taken for cosmic rays to stream along a field line at the speed 
of light to reach the surface of the Disk. Consequently, the cosmic ray 
spectrum should regain its original shape above unless there is a 
fortuitous coincidence of the right amount of steepening in the injection 
spectrum at E^.. In order to overcome these difficulties, it has been 
suggested that there could be energy dependent confinement of cosmic rays 
in some local regions in the Galaxy. 

(ii) Energy dependent confinement of cosmic rays in special regions 

If cosmic rays diffuse away from their sources in an energy dependent 
fashion, the matter traversed by them would also be so dependent provided 
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most of it is traversed within the source, where the density of matter 
is expected to be very large. By suggesting that cosmic rays above 3 
GeV/n have an energy dependent confinement of the type with Xg = 21/E g cm"^ 
in the source regions, Meneguzzi (1973b) was able to reproduce closely the 
observed spectral shape of nucleonic components; by virtue of this effi- 
cient confinement of low energy cosmic rays, these sources are expected to 
be sources of gamma rays through neutral pion decay. On the other hand, 
Cowsik and Wilson (1973) have speculated that the effective matter trav- 
ersed during diffusion in the source region varies as Xg I 1,7 exp (-E/ 

7.85) g cm"^ with a subsequent traversal of about 1.6 g cm“^ of interstellar 
matter during their confinement in the Galaxy. While these models over- 
come the difficulties of the earlier ones, they suffer from the following 

difficulties. First, these models assume exponential path length distrl- 

inside the source region 

butions with allowed zero path lengths/ as in the case of steady state models 
(which assumes uniform source distribution in an infinite diffusive medium) 
instead of using diffusion equations characterizing the diffusion of cosmic 
rays from a point source; in the latter case, there should be a minimum 
finite amount of matter traversed by cosmic rays inside the source leading 
to a large depletion of ultra-heavy elements. Secondly, the energy depen- 
dent diffusion in the source would lead to a steep steady state spectrum of 
cosmic rays inside the source region. In consequence, the resultant steep 
electron spectrum would give rise to synchrotron emission with a spectral 
index >0.8, which is in contradiction to the observed radio spectral in- 
dices of supernova remnants which is <0,5 below about 5GHz (Milne and Hill, 


1969) . 
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Rengarajan, et al. (1971) suggested that most of the matter traversed 
by cosmic rays Is in the dense clouds, where the fraction of time spent by 
them is about 5 times larger than the fractional volume occupied by clouds 
in the Galaxy. They also introduced (Rengarajan, et al. , 1973) an energy 
dependent streaming velocity for cosmic rays in the cloud and tried to 
explain the observed energy dependent charge composition. Similar to the 
source model, the dense clouds are expected to be strong sources of pionic 
gamma rays. The main difficulty of this model is that it is hard to 
achieve under steady-*-state conditions, a higher density of cosmic rays in- 
side the clouds than in in ter cloud medium, from where they propagate into 
the clouds. However, one cannot rule out the possibility that the gas 
density in interstellar medium fluctuate and, as a result, cosmic rays and 
magnetic fields are compressed and later released through ampipolar diffu- 
sion as required by this model. 
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4,5 Propagation of the electron component and its source characteristics 
Cosmic ray electrons in interstellar space consist of two components: 
(i) the directly accelerated electrons and (ii) the secondary electrons. 

The directly accelerated electrons are generally assumed to be negatively 
charged because of various reasons including the evidence for the negligible 
intensity of antiparticles in cosmic rays and the fact that the fraction 
of positrons observed at energies greater than a GeV seems to be consis- 
tent with that of secondary production. The secondary component would 
consist of positive and negative electrons and their relative abundances 
vary as a function of energy as seen below. In order to understand the 
propagation and source spectra of cosmic ray electrons, it is essential 
first to determine the characteristics of secondary electrons resulting 
from collisions of cosmic ray nuclei with interstellar gas, 

4,5.1 Secondary electrons in the Galaxy : 

The well known processes through which secondary electrons are pro- 
duced in the Galaxy are (i) decay of short lived particles like charged 
pions, (ii) decay of neutrons and radio nuclei through beta decay and 
(iii) knock-on of atomic electrons by fast moving particles. We shall 
briefly indicate the methods of calculating the production spectrum of 
secondary electrons through the above processes, 

(i) Decay of short lived particles : The most abundant of the created 

particles produced in interactions of high energy cosmic ray nuclei 
with matter, are pions. The charged pions decay to electrons through 
muons as 


IT ->■ M" + 


K 


e- + + Vg 


4.1A 
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Since the decay life times of these particles are very small compared to 
the time scales involved in the propagation of cosmic rays, one can assume 
the decays to be instantaneous in cosmic scale. The production spectrum 
of pions has been derived by many authors in the past as reviewed by 
Daniel and Stephens (1970) and recently by Chocate and Wayland (1972) 
and Dooher (1973). However, it needs to be mentioned here that the re- 
sults of these calculations are not in as good agreement with each other 
as one would like it to be mainly because the experimentally determined 
production spectrum of pions in the laboratory system, is not available 
at all energies of relevance here. 

We will now start with the production of pions (W)dW, which is 
defined as the number of pions of total energy W in the interval dW pro- 
duced per gram of hydrogen in interstellar space per second, to derive 
the production spectrum of electrons; here it denotes both posi- 

tive and negative pions unless stated otherwise. It needs to be pointed 
out that the contribution of electrons from the decay of other short 
lived particles like K-mesons can be taken care of by suitably adjusting 
the parameters in Qqy(W) . The production spectrum of muons can then be 
written as 

dW 

Q (W)dW = / Q (W’)dW’ 4.15 

W’ = W +E ^ 

Here e = m - m and il» (W) = W(l-f) , where f is the minimum fractional 

IT y y 

energy the muon gets from the pion and is given, by the relation 

f = ^ [(1 + m^ /m^ ) - (1 - m^ /m^ ) /I - m^ Ve 4.16 

*• y IT y IT IT 7T 

which is equal to 0.56 for >> m^. The production spectrum of elec- 


trons is then given as 



no 


Q^(W)dW = z"' Q^(W')dW* ip^(W’,W) dW 4.17 

where e = W for W > m and e = m when W < m ; (W* , W)dW is the proba- 

y y ye 

bility that, during the decay of muons of energy W* , the electron gets 
an energy W in dW and can be expressed as (Zatsepin and Kuzmin, 1962). 

i|;^(W , W)dW - ^r|)T [3(1-6^) - (3+b2)W/W^1 (WAJ^)2 dW/W^ 


for 0 < W/W 1 (l-e)/2 4.18a 

= ? <1 + 'f 

for (l-e)/2 1 W/W^ < (l+B)/2 4.18b 

where 8 is the velocity of muons relative to that of light. 

(li) Decay of neutrons and radio nuclei : Neutrons are produced mainly 

through charge exchange processes during proton collisions and the break 
up of helium nuclei, while radio nuclei result from the fragmentation 
and evaporation of incident and target nuclei respectively. The pro- 
duction spectrum of secondary neutrons and radio nuclei Q (E)dE per gram 


of interstellar material per second can be written as 

Qg(E)dE = 4 tt £ /dE^ ik ^^i^ 4.19 

i,k 


where E and E are the kinetic energies of the incident and secondary 
i 

nuclei respectively, is the flux of primary cosmic ray nuclei of the 
type i in interstellar space, n^ is the number of nuclei of the type k 
per gram of interstellar gas , a ..is the cross-section for the produc- 
tion of a particle of type s and ij; is the normalized energy distribution 
of the secondary particle during the collision. The cross section for 
neutron production and the mean energy E^ = / E)E dE as a function 

of E^ are summarized by Ramaty and Lingenf alter (1966a) ; the integral in 



Ill 


Equation 4,19 can then be replaced by substituting = fE, where E is the 
mean energy. Knowing the production spectrum of unstable nuclei , one 
can calculate the production spectrum of electrons using the relation 
(Ramaty, et al. * 1970b) 


1 1 /■ j.\ y ^ 


(Y)dY = z I dY / My . *) / + dT . : ... qs . (Tj 


4.20 


8 1 / Y ^- / y '^-1 

where y* , y» y* are the Lorentz factors of the unstable nuclei, electrons 
in the laboratory system and electrons in the rest system of the parent 
nuclei respectively, y^ is the maximum energy of the electrons in the 
rest system of the parent nuclei and ^^(y*) is the Fermi distribution 


function (Fermi, 1934): 


ip^(y*) = Y*(y*^ *“l)^ (y_ ~ y*) 

C lu 


4.21 


The constant g is the normalization factor for the Fermi distribution 


y 


m 


defined as g / ^|^(y*) = p, where p is the probability of electron decay 

1 

through the given mode. The limits in Equation 4.20 are given as 


= Y*y ± (y*^ - 1)^ (y^ - 


4.22 


The lowest value of y* is ' also determined by the minimum energy upto 
which the cosmic ray nucleon spectrum is ass\jmed to exist. Some of the 
radio nuclei like 0^^ and 0^^ decay through positron mode and 

hence positrons are also produced in this process. 

(iii) Knock’on process : In this process, energy is transferred to 

electrons of interstellar atoms in Coulomb collisions of cosmic ray 
nuclei. The cross section for knock-on production has been calculated 
by Bhabha (1938) and according to which, for spin 3^, we can write the 
differential probability, that a cosmic ray nucleus of charge and 
total energy per nucleon W* , collides with an atom of charge and 
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atomic mass A m giving rise to knock on electrons of total energy W, as 

ip 

(Daniel and Stephens, 1974) 

o 'i e 


W)dW = 


2itN Z. Z? m 


A. (W'2 - m^) 
1 P 


L (W “ m 




2 A.m W' 4- ra^ + Aj m^ - 'I 

i-s ^ ^ +— dW cm2 g-l 4.23 

2 aJ m (W - m ) 2 A? J 

1 e e' 1 

where. N is Avogadro’s number and r is the classical radius of the elec- 
o ® 

tron. The production spectrum of knock on electron can be written as 

Q (W)dW = 4ir /“ E Z <j> (W , W)dW j (W*)dW^ 4.24 

® W*=e i k 

th c 

where i (W’) is the flux of cosmic ray nuclei of k. type. One can, for 
k 

practical purposes , replace the summation sign from the above equation by 
a factor of 1.82 to take care of the composition of interstellar gas and 
cosmic ray nuclei, and defining all the quantities in terms of protons. 

We have shown in Figure 4.17, the differential production spectrum of 
secondary electrons per unit gram of hydrogen (after taking into account 
the composition of interstellar gas) , per second through various processes 
described above. The curves in this figure have been obtained as follows; 
(a) the knock on spectrum (Curve A) is the modified version of the spec- 
trum given by Abraham et al. (1966) ; (b) the electrons from neutron decay 
(Curve B) is taken from Ramaty and Lengenfelter (1966a) , (c) the pion decay 
spectra of both negative (Curve C) and positive (Curve D) electrons below 
a GeV have been evaluated from the work of Ramaty and Ligenfelter (1968) 
and Perola et al. (1967) and above 1 GeV from Stephens (1969) by assuming 
a charge ratio e+/e" 1.3 same as that of the observed muons at high 

energies ; and (d) the positron spectrum from radio nuclei is from the work 



ELECTRONS(g sec GeV) 
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of Ramaty, et al. (1970b) and Venua (1969)* One can see from this figure 
that at energies below 30 MeV, the knock on spectrum dominates over 
other e~ spectra, while at energies above 50 MeV, electrons from pion 
decay dominate. Though the e" spectrum from radio nuclei is not included 
here, it can be easily shown from arguments similar to those given in (ii) 
that their contribution would be much smaller than that from neutron de- 
cay at energies above a MeV, The positron spectrum is dominated by pion 
decay at all energies except below 2 MeV, where the decay of radio nuclei 
is important. The total secondary production spectrum is shown by curve 
T; it cannot be represented by a simple power law because of the many pro- 
cesses which contribute at different energies. However, the production 
spectrum below 10 MeV can be represented by a power law of the type 

Qg(E)dE = dE electrons g"^ sec“^ 4.25a 

where E is expressed in MeV; also the spectrum above a GeV can be repre- 
sented by another power law of the type 

Q (E)dE = 4,6 X 10“^ dE electrons g“^ sec”^ 4.25b 

0 

where E is expressed in GeV, 

When once the secondary electrons are produced, they lose energy 
through ionization, bremsstrahlung, synchrotron radiation and inverse 
Compton scattering during their propagation in interstellar space; some 
of them will also diffuse out of the confinement volume of space. We 
shall now calculate the equilibrium spectrum of electrons under the known 
astrophysical conditions in the Galaxy by incorporating in addition a 
continuous acceleration in Interstellar space. The equilibrium spectrum 
can be obtained by solving Equation 3^29; the values of the parameters a, b, 
a and t in this equation can be estimated in the following manner. 



TABLE 4.3 


Region 

J)isk 

It 

Halo 

II 


Parameters used in calculating the equilibrium spectrum of electrons 

Njj a b T a Curve in 

atoms cm~ ^ GeV sec~^ GeV“^ sec~^ sec sec~^ Fig« 4.18 

1.0 2.22 X lO’l^ 1.75 X 10“16 6 x 10^^ «l/x A 

II II 7 75 X 10-16 «i H B 

" « ” ” l/t^l-eTxlO^l*^ C 

1.0 X 10~2 2.22 X 10-18 7,5 x 10“17 6 x IQI^ «i/t D 

*» " ” ti 1 /t= 1. 67x10-18 E 



116 


(1) Since ionization is of importance only below about 20 MeV, 
a 2,22 X 10"^^njj GeV sec, where ng is the mean number of hydrogen 
atom cm”^ in Interstellar space, (ii) The value of b can be written 
as b ^ 10-^^ X (3.8 X 10“^ fif + p^) GeV“l sec“^, where is the mean 
perpendicular component of the magnetic field in pG and the total 
photon density in eV in the Galaxy, (iii) The residence time t is 
taken to be about 2 x 10^ yrs for the mean matter traversal of about 
3 g cm”^ of hydrogen by cosmic rays, (iv) The value of the accelera- 
tion parameter a is kept as a free parameter. Table 4,3 summarizes 
the parameters used in evaluating the equilibrium spectrum of secondary 
electrons. It needs to be mentioned here that in the case of the Disk 
the existence of the submillimeter radiation has not been taken into 
account in calculating the value of b for Curve A while for Curve B it 
has been included; for the Halo its existence has been neglected. 

In Figure 4,18 the equilibrium spectra of electrons are shown for 
various set of parameters summarized in Table 4.3. In the absence of 

acceleration, the spectral shape of the equilibrium electrons remains 
as the production spectrum (Figure 4.17) 
the same /in the energy region between about 20 MeV and a few GeV (Curves 

A, B and D) , Below 20 MeV, the spectrum is flatter due to ionization 

and 

loss (Equation 3. 32a) /above a few GeV, the spectrum steepens due to 
synchrotron and inverse Compton processes (Equation 3,32c); the energy at 
which this steepening occurs is inversely proportional to the product 
of the residence time and the parameter b. It can also be seen from 
Curve A that continuous acceleration is not at all effective for values 
of ot ^ »/4aF. For a > /4ab , it is only marginal if a << 1 /t even in the 
energy region a - - 4ab < 2bE < a + /ot*^ - 4 ab , where acceleration 
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Fig , 4.18 ; Equilibrium spectra of' secondary electrons in the 
Galaxy for various interstellar parameters summarized in Table A. 4; 
the observed intensities below a few tens of MeV are also shown. 
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dominates over continuous energy loss processes. However, the acceleration 
becomes important as the value of a approaches 1/t and when a = 1/t, the 
equilibrium condition is never reached at energy Eg = (a + /a^ - 4ab)/2b 
and the flux increases as the ratio of the age of the confinement volume 
of space to x. These calculations clearly demonstrate that if interstellar 
acceleration is effective, the flux of secondary electrons is very much 
enhanced in the region above a few hundred MeV. However, observations on 

f 

positrons do not indicate any such trend and hence it may be stated that 
continuous acceleration in interstellar space is either absent or ineffective, 
4.5.2 Directly accelerated electrons 

In order to deduce information on the directly accelerated electrons , 
it is essential to know first the equilibrium spectrum of the total cosmic 
ray electrons in interstellar space. Thereafter, by subtracting the contri- 
bution from secondary electrons, one can derive the intensity and spectral 
shape of directly accelerated electrons under equilibrium conditions 
(Daniel and Stephens, 1970a). The energy spectrum of the electron compo- 
nent in interstellar space can be obtained by different methods; (i) One 
can demodulate the observed electron spectrum by making use of the modula- 
tion parameters, deduced from an analysis of the temporal variation of the 
nucleonic and electronic components and a study of the propagation of nucle- 
onic component; (ii) It can be deduced from the galactic radio background; 
and (iii) It can also be inferred from a comparison of the observed posi- 
trons and the expected secondary positrons in interstellar space. 

In Figure 4.19, Curve A is the interstellar radio emitting electron spec- 
trum as deduced in Section 5,4.2 by using method (ii) and Curve B is the 
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Fig. 4.19 : The radio emitting electron spectrum in the Galaxy 
(Curve A) is compared with that deduced from the observed positrons 
(data points) ; Curves B and C are the equilibrium spectra of secon- 
dary and directly accelerated electrons respectively. 
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calculated equilibrium spectrum of secondary electrons in interstellar space. 
The points in Figure 4.19 are the calculated Intensities j^(E) of electrons 
in interstellar space using the observed fraction of positrons f^(E) as a 
function of energy (Beuermann, et al. , 1970; Fanselow, et al,, 1969) and 
the relation jg(E) = j^.(E) /f_j.(E) , where j^(E) is the estimated secondary 
positron flux in the Galaxy. It is seen from this figure that at energies 
> 200 MeV, the total electron spectrum deduced from method (ii) is in good 
agreement with that obtained by method (iii) . Hence, it is reasonable to 
extend Curve A below 200 MeV on the basis of the observed positron fraction; 
this is shown by the dotted curve in Figure 4.19. 

Having deduced the equilibrium spectrum of total electrons in inter- 
stellar space, the energy spectrum of directly accelerated electrons can be 
obtained by subtracting from the above the calculated spectrum of secondary 
electrons. This method, first employed by Daniel and Stephens (1970a), 
demonstrated that the directly accelerated electron spectrum flattens 
gradually below 2 GeV such that the flux probably becomes insignificant at 
very low energies; the revised spectrum as obtained from the present data 
is shown as Curve C in Figure 4.19. This spectrum is found to dominate over 
that of secondary electrons (Curve B) above 2 GeV by a factor of about 10, 
having a spectral index 3 - 2,6 as required by the radio spectra (section 
5.4) and is similar to those of nucleonic components. Below this energy 
the contribution of secondary electrons slowly increases and perhaps over- 
whelming below about a few tens of MeV, where the calculated secondary elec- 
trons of Figure 418 accounts fully for the observed flux; a similar conclu- 
sion was drawn recently by Cummings, et al. (1973a), 
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4,5.3 The origin of cosmic ray electrons 

The relatively low abundance of electrons compared to protons in 
cosmic rays demands an explanation as to why charge neutrality, which is 
universal in all natural phenomena, is not preserved in cosmic radiation. 
However, a satisfactory explanation to this most intriguing problem is still 
to be proposed in the future; therefore in what follows in this section, 
we shall briefly outline our present knowledge on the origin of cosmic ray 
electrons in a number of suitable energy intervals. 

(i) e^ below 1 MeV : The existence of cosmic ray electrons with energies 

E < 1 MeV in interplanetary space was first detected by Beedle , et al. , 
(1970) , who showed that the flux of electrons in the region of a few htm- 
dred keV seems to fall close to the extension of the observed spectrimi 
above a few MeV. Recently these measurements have been extended down to 
energies about a few tens of keV (Lin, et al. , 1972; Webber, et al., 1973d; 
Hurford, et al. , 1973a) and from the observed temporal variations, elec- 
trons above 100 keV are considered to be of galactic origin. One notices 
from the inset of Figure 4. 18 that at energies < 1 MeV, the flux values of 
Lin, et al. , are about an order of magnitude larger than those of Hurford, 
et al. ; this large steady flux of Lin, et al. , could be due, as pointed out 
by Hurford, et al. , to a substantial spacecraft background and long-term 
averaging. Since the electron spectrum of Hurford, et al. , and Webber, et 
al., are consistent with the extension of galactic knock-on spectrum (Figure 
4.18), one is inclined to consider these low energy electrons to be of 
secondary origin; the present situation therefore doesn't warrant any new 
mechanism to be invoked (Ramaty, et al. , 1972) to understand the observa- 


tions , 
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(ii) e~*~ around 1 MeV : In the energy region between 0.5 and 2 MeV, an 

anomalous finite flux of positrons was reported by Cline and Hones (1970) , 
while the observed positron intensity above 3 MeV is down by a factor of 
about 10^ (Cline and Porreca, 1970). In Fig. 4.20, we show the observed 
intensities of positrons from about 200 keV to 20 MeV (Cline and Hones, 
1970; Cline and Porreca, 1970; Beuermann, et al. , 1970; Hurford, et al. , 
1973b) along with theoretical estimates. Though recent measurements of 
Hurford, et al., provide only upper limits, the suggestion of a sharp de- 
crease of the intensity between 2 and 10 MeV indicate that perhaps these 
positrons arise from the decay of radio nuclei. Estimates have been made 
of the positron spectrum from radio nuclei, produced during the collision 
of cosmic ray nuclei with interstellar gas for various assxmied primary 
spectra (Verma, 1969; Ramaty, et al. , 1970b). Curve A in Figure. 4.20 Is a 
calculated spectrum by Ramaty, et al., for a power law spectrum of cosmic 
rays in kinetic energy to exist in interstellar space up to an energy = 

5 MeV/n. This explanation, however, requires a cosmic ray energy density 
of 'vlOO eV/cm^ in interstellar space, and a modulation factor vio to be 
operative for positrons only up to an energy just above 2 MeV. 

An alternate explanation was put forward by Burger, et al. (1970) and 
Colgate (1970), who suggested that these positrons could arise from the 
decay of ^^Co ^^Fe in the silicon burning shells of supemovae just after 
their ejection at relativistic energies. Curve B in Figure ^*20 is the calcu- 
lated equilibrium spectrum of positrons by Burger, et al. , using a power 
law spectrum in rigidity for the ^^Ni, which first decays to ^^Co, in the 
source region with R^, - 30 MV; Curve C is for = 100 MV when ionization 
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Fig, 4.20 ; The observed positron spectrum is compared with the 
theoretical spectra based on their production in Interstellar gas 
(Curve A) and in supernova envelopes (Curves B and C) , 
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loss is not taken into account. This explanation requires that most of 
the iron nuclei in cosmic rays result from the decay of ^%i, which are 
produced during efficient burning of silicon in the supernova shells 
(Bodansky, et al. , 1968), Subject to the reality of the above situation, 
this hypothesis is more attractive than the earlier one (Curve A) both 
from the point of view of better agreement with observations and more 
plausible assumptions regarding the primary spectrum. 

(iii) e~ between 1 and 10 MeV : The observed intensity and charge composi- 

tion of electrons in this energy region can be explained as due to the com- 
bination of interstellar knock-on and pion decay electrons (Figure 4.18). It 
needs to be pointed out at this stage that it is assumed here that at 
energies < 10 MeV the observed intensity truly represents the interstellar 
spectrum because of the absence of any detectable long-term variation of 
their intensity (Simnett and McDonald, 1969) . 

(iv) e~ between 10 and 100 MeV : In this energy region, the electrons 

undergo modulation and one is unable to fully account for the Intensity 
on the basis of the secondary hypothesis. Further, the observed charge 
composition in this energy region suggests that a comparable fraction of 
electrons in this energy range is directly accelerated (Figure 4,19). 

(v) e~ of energy > 100 MeV : Electrons in this energy domain. are .very 

important since directly accelerated electrons influence the spectral 
shape and account for about 90% of the total electrons above a few GeV. 

One would therefore like to examine whether the spectral shape of directly 
accelerated electrons in equilibrium is representative of the injection 
spectrum or is modified due to energy loss processes. Many arguments have 
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been put forward in the past, as summarized by Daniel and Stephens (1970), 
that the gradual steepening of the electron spectrum between a few hundred 
MeV and a few GeV is not due to energy loss processes but intrinsic to the 
injection spectrum. Further, one notices from Figure 4.19 that the spectral 
shape of directly accelerated electrons gradually changes by a power index 
I A3 1 ^1.5, which cannot be understood on the basis of energy loss processes 
in interstellar space. In this regard an accurate determination of the 
positron spectrum experimentally is very valuable because one can compare 
this with the production spectrum to look for any propagation effects. Not- 
withstanding this, the study of the electron spectrum above a few GeV plays 
an Important role in the understanding of the propagation of cosmic rays 
and we shall dlsc\iss this aspect in some detail. 

4,5.4 Interpretation of high energy electron spectrum 

The available data on high energy electron spectrum above a few GeV, 
published after 1970 are plotted in.: terms of energy flux units for con- 
venience in Figure 4. 21(Hovestadt, et al. , 1971; Marsden, et al. , 1971; 
Scheepmaker and Tanaka, 1971; Swanenburg, et al. , 1971; Zatsepin, 1971; 

Earl, et al., 1972; Anand, et al. , 1973a; Burger and Swanenburg, 1973; 

Muller and Meyer, 1973; Meegan and Earl, 1973; Silverberg, et al. , 1973; 
Ishii, et al. , 1973; Webber and Rockstroh, 1973): the earlier results have 

previously been summarized by Daniel and Stephens (1970) , Yash Pal (1969) 
and Meyer (1971) . One notices from this figure that the spread of the data 
points due to various experimenters is disproportionally large, compared to 
what is expected on the basis of statistical errors attributed to the indi- 
vidual values, thereby clearly demonstrating the existence of large 
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Fifi> 4,21 ; The observed energy spectrum of cosmic ray electrons of 
energy greater than a few GeV. 
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TABLE 4.4 

Observed Differential Spectral Index of High Energy Electrons 

Energy range Spectral Index 


Authors 

in GeV 

B 



Swanenberg et al. (1971) 

3-10 

2,6 

+ 

0.1 

Marsden et al, (1971) 

3-15 

2,6 

+ 

0.1 

Earl, et al. (1972) 

5-25 

3.3 

± 

0.1 

Webber and Rocks troh (1973) 

3-20 

3.0 

+ 

0.2 

Scheepmaker and Tanaka (1971) 

5-300 

2.7 

± 

0.1 

Zatsepin (1971) 

5-600 

2.7 

± 

0.1 

Muller and Meyer (1973) 

30-1000 

2.7 

± 

0.1 

Anand et al. (1973a) 

10-800 

2.7 

± 

0.1 

II II 

>150 

3.1 

+ 

0.3 

Ishli et al, (1973) 

>100 

3.2 

+ 

0.2 

Meegan and Earl (1973) 

10-100 

3.5 

+ 

0.1 

Silverberg et al, (1973) 

10-200 

3.2 

+ 

0.2 
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systematic errors between them. This is again brought out in Table 4.4, 
where we have summarized the spectral indices due to the same authors. It 
is, therefore, neither correct nor possible to look for any steepening in 
the spectral shape from a plot of "World data." However, one can still make 
the following statements about the electron spectrxim from Figure 4.21 and 
Table 4.4: (a) the cosmic ray electron spectrum extends at least up to 

about 700 GeV with non-trivial flux values; and (b) there is no indication 
of a detectable change in the spectral index up to 100 - 200 GeV. We shall, 
therefore, attempt now to bring out the possible implications of the above 

observations on the basis of the theoretical models described in Section 

3.5.2. 

source model : In the case of cosmic ray electrons originating 

from a single source, one can notice . ffoBl Equation 3.34 that electrons cannot 
reach the point of observation with energies E ~ l/b(t ~ t^) . Since we 
observe a non-trivial flux of cosmic ray electrons at least up to 700 GeV, 

(t - to) < 3 X 105 yrs for b V 5 X 10^ GeV^l sec~l (Table 4.3); a similar 
result has been arrived at recently from an analysis of the composition of 
UH nuclei (Kaiser, et al. , 1972). From the distribution of non-thermal 
radiation in the Galaxy, we infer (Section 5.4.4) that the intensity of 
cosmic rays in the neighborhood of the Sun is nearly the same as that else- 
where in the Galaxy and in order to achieve this, electrons should have 
enough time to diffuse to the point of observation from the source region. Thus, 
according to Equation 3.34, 4D(t - t^) > | r - r^|^ and taking Vela X to be 
the nearest plausible source of cosmic ray electrons (Shen, 1970), one can 
set a lower bound of 4 x lO^B cm^ sec"! for the diffusion coefficient D. 
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From the above it is quite clear that if a single source contributes to 
the observed cosmic ray electrons, then the mean free path for cosmic rays 
for isotropic diffusion is A > 1 pc which in turn leads to large anistropy. 

(ii) Flat disk model : Let us now consider the flat disk model in which 

sources are uniformly distributed and cosmic rays are in a state of equi- 
librium. In this case, one expects from Equations 3.36 and 3,37 two spectral 
steepenings, each with [ab| = 0.5 at energies and E 2 separated by a 
factor E 2 /E 2 R^/L^. It has already been mentioned in Section 4,5,3 that 

the observed electron spectrum shows a gradual steepening at about 2 GeV 
and in spite of -the fact that this spectral shape could be inherent to the 
directly accelerated electrons , one can still examine the consequences by 
ascribing this steepening to result from propagation effects under two 
possibilities. In the first, it may be said that this is the first spec- 
tral break corresponding to energy and from the observed lack of any 
detectable spectral change at least up to 'v.200 GeV, one could say that 
R^/H^~ 100. If one sets 2H 350 pc, the effective thickness of the gas 
disk, then the sources in the Disk which contribute to the observed radi- 
ation have to lie within a radius R - 1.75 kpc from the Sun. Further, from Equa- 
tion 3. 36 one gets a lower limit to the diffusion coefficient D ~ E-|^bR^/4 - 
2 X 10^^ cm^ sec“^ and a corresponding residence time t 'v H^/4D - 10^ yrs. 

Though these values seem attractive, there are a number of objections to be 
suitably explained. If, for example, the observed spectral index above a 
few GeV is close to 2.6' (Table 4,4), it means that the injection spectrum 
is significantly flatter than that of cosmic ray nuclei carrying with it 
other seriotis implications (Daniel and Stephens, 1970), On the other hand, 
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if It Is claimed that the spectral index is close to 3.1 (Table 4,4), then 
one runs into the following difficulties : (i) the radio spectral index 

should gradually steepen to a value of (3 - l)/2 ~ 1.05 beyond about a few 
hundred MHz; the available data suraniarized in Section 5,3,1 shows that the 
radio spectral index has a value of 0.8 from about 178 MHz to a few GHz, 

(ii) The positron spectrum should also have a spectral index of 3.1 beyond 
2 GeV; Indications from the presently available observations (Fanselow, et 
al, , 1969) are not consistent with this, though better measurements are 
needed for a final decision in the future. Finally, this explanation of a 
break at 'v2 GeV also requires a second break of |A3| = 0.5 to be established 
beyond a few hundred GeV. 

If we consider the other alternative that the apparent spectral change 
around 2 GeV is due to the second break at £ 2 * then from Equation 3.36 we have 
D 'v 2 X 10^^ cm sec“^ and x ■X' 10® yrs. These values are not compatible 
with. our knowledge of the propagation of cosmic rays in the Disk, Further, 
one requires that (i) the sources which contribute to the observed cosmic 
rays are distributed uniformly within a radius of about 175 pc, and (ii) the 
energy spectrum of electrons at injection should have a spectral index ~2.0, 
which is very different from that of cosmic ray nuclei, and (iil) the posi- 
tron spectral index has to be ~ 3.6 in contradiction to the observation. 

Thus, we find that this explanation is less plausible than the previous one. 
Perhaps before concluding, mention may also be made of a third but less 
attractive possibility that even the first break has not occurred at least 
up to 200 GeV. 
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(lii) Leakage lifetime model : On the basis of this model, the residence 

time of cosmic rays is about 10® yrs and - 10® yrs corresponding to the 
possible breaks by one power around 2 and 200 GeV, respectively. However, 
arguments mobilized in the previous paragraphs in this section would favor 
the second possibility. Furthermore, it is also possible to introduce an 
energy dependent leakage lifetime for electrons as in the case of nucleonic 
components (Silverberg and Ramaty, 1973). Accordingly, if one assumes a 
dependence of the type t a E®*® beyond about 5 GeV (Curve B in Figure 4.13), 
one gets for the lifetime of cosmic rays in the GeV region a value ~3 x 10® 
yrs on the basis that the electron spectrum doesn’t show a finite steepen- 
ing up to 200 GeV 5 some of the difficulties of this model have been summar- 
ized in Section 4.4.3. 

Thus , we find that there are many ways of interpreting the observed 
data on high energy electrons. It is evident that, though in principle 
some of the alternatives could be ruled out, it requires for this purpose 
more accurate and reliable observations on the spectra of both positrons 
and negatrons. In spite of the existing uncertainties, the overall evidence 
is in favor of a cosmic ray residence time 'vlO® yrs. It is also interesting 
to mention here that, since electron energy loss through synchrotron radi- 
ation is maximum when its pitch angle is perpendicular to the magnetic field 
lines, the high energy electrons in the local interstellar space would be 
anisotropic with bi-directional maxima defining the local direction of the 
large scale magnetic field (Earl and Lenchek, 1969). 
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4.6 Leakage Life Time of Cosmic Rays 

The leakage life time (also differently referred to as dwell time, 
residence time, confinement time and storage time) of cosmic rays is 
defined as the mean amotmt of time for which the cosmic rays reside 
within the confinement volume of space; it is characterised by the 
nature of the confinement volume, the distribution of sources within 
this volume of space and the diffusion mechanism. When diffusion is 
rapid and the storage volume has a reflecting boundary which inhibits 
free leakage, the mean age of cosmic rays sampled at any place within the 
confinement volume will be the same and will represent the residence time; 
also it will not depend upon the distribution of sources around the observer, 
except in case of a close proximity of a source or sources to the observer 
both in time and space. On the contrary, when diffusion is slow, the 
residence time is a meaningless concept because the mean life time of 
cosmic rays sampled by the observer is influenced more by the energy loss 
associated with the type of particle sampled and the distribution of 
sources around the point of observation, than by the nature of the confine- 
ment volume of space. In this section, we shall make use of the concept 
of leakage life time to understand the confinement voltme of space by 
attributing the mean age of cosmic rays sampled near the Earth to be 
identical with the residence time of cosmic rays. 

One powerful method of deducing information on the leakage life time is 
to take advantage of the sensitive dependence of the energy spectrum of 
cosmic ray electrons - among other factors - the residence time of cosmic 
rays. However, since this aspect is adequately covered in Sections 3.5.2 
(theoretical aspects) and 4,5.4 (interpretative aspects) we will not go 
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into it here. One other important aspect of leakage life time which 
we have avoided mentioning in this section for a similar reason is the 
following. The interesting observations recently made from a study of 
the energy spectra of "primary” and "secondary" cosmic ray nuclei in 
the energy range S-rlOO GeV/n, and the many theories proposed to under- 
stand these observations on the basis of an energy dependent confinement 
have been described in detail in Section 4.4. 

4,6.1 Mean Matter Traversed by Cosmic Rays and the Residence Time . 

It is a normal practice to relate, <x> the mean amount of matter 

traversed by cosmic rays, to the mean age of cosmic rays. Consequently, 

one can write x = <x>/pc, where p is the mean matter density in the 

confinement volume. It has been shown in Section 4.1, that presently 

—2 

available observations suggest that <x> 3-6 g. cm of hydrogen in the 

GeV region. Using this value for <x> , one gets x % (2-4) x 10^ yrs for 

—24 —3 

the Disk model where p = 1.67 x 10 g.cm of hydrogen while for the 
Halo model one gets x ^ (2-4) x 10® yrs, where p % 1.67 x 10 g.cm 
of hydrogen. These estimates of the leakage life time are justified 
under the assumption that the matter traversed by cosmic rays 1» In inter- 
stellar space and that matter is uniformly distributed over the entire 
volume of the storage space. However, it has been shown in Section 2,3 
that about 70% of the interstellar gas is in the form of dense clouds. 
Even so, if cosmic rays spend equal amounts of time in equal volumes of 
space, the above estimate of x will not be affected by the exact nature 
of the interstellar gas. Notwithstanding this, it has been found that 
magnetic fields in the clouds are high and are effectively detached from 
the intercloud medium. If these were so, one would expect cosmic rays 
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to be confined inside clouds for a longer time than in a corresponding 
volume of Intercloud medium; this situation would then tend to reduce 
the value of t. In a similar manner, if part of the matter is traversed 
In the source region, one would over-estimate the value of t by attri- 
buting the entire matter traversed to be in interstellar space. Thus, 
in either case, the values of t deduced above are likely to be upper 
limits. It needs to be emphasized at this stage that in this method 
T can be estimated only by assuming a model for the cosmic ray confinement. 
4.6.2 Leakage Life Time From Long Lived Radio Nuclei in Cosmic Rays 

A direct method of estimating the mean age of cosmic rays is to 
determine, in the primairy radiation, the abundance of a suitable long 
lived radio nucleus produced in secondary collisions. It has been 
suggested by Hayakawa et al (1958) and Peters (1963) 'that Be^^ , which 
decays to with a mean life of 1.6 x 10^ yrs , could be used as a 
clock to measure the mean age of cosmic rays. The first attempt towards this 
was made by Daniel and Durgaprasad (1966) and later by many others, by 
making use of the ratio Be/B experimentally determined. We briefly 
describe first the method of evaluating -the. age of cosmic rays using a 
radio active tracer and then summarize the present status in this subject. 
Since the estimate of leakage life time depends upon the propagation 
model employed in determining the flux of radio nuclei and their decay 
product, we shall derive some useful relations for the flux of nuclei 
tising slab model and steady state model. Let radio nuclei of the type 

i , where the subscript denotes the given isotope of element i, decay to 
m 

nuclei 1 with a mean decay time T at rest. We shall assume that the 
m 

attenuation mean free path due to interaction is the same for all 
isotopes of the element i. In order to take into account the decay 
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probability of the radio nuclei, one has to add another term in the 
general diffusion equation. In the case of slab model we can add a 
term j(E,x)/s(E) to Equation 3.23, where s(E) - T c3y* and the solution can 
be written as 


j (E,x) = !: r* (E'.x')|/ |w. (E)|} {j,(E",x')/X.(E")) 

- k>i o*' - ^ ^ 


m 


m 


m 


^(E”) exp[- / {l/A^(E*z) + l/s(E*))dz]dx» 4.26 

m X* 

Here, eg and y are the velocity and Lorentz factor of the particle respectively 
Similarly, the flux of 1^ nuclei which are produced in Interstellar 
space, directly from fragmentation of heavier elements and also through 
the decay of i , can be written as 


(E,x) = f {|w^ (E',x')|/|w^ (E)|) Z {j (E",x') P» (E")/X (E") 
mom m k>£ m 


+j^ (E",x*)/s(E")} exp {- / dz/A£<E»,z)}dx' 
m X 


4.27 


The first term in the above equation is same as the secondary production 

term in Equation 3.24 and (E",x*) in the second term is defined by 

m 

Equation 4.26. n the case of steady state model, the flux of nuclei 

(E) can be obtained from Equation 3,26 by Introducing an additional 
m 

attenuation due to decay. Thus 

ft ** X 

^i ^ ^"i (E*,x)|/|w^ (E)|}Q^ (E') exp {- / 


[l/A^(E',y) + 1/x + l/s(E')3-dy}dx 


4.28 
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Where Q. (E)dx is the production spectrum of the nucleus 1 per unit gram 
m 

of interstellar matter, which can be estimated from Equation 4,19. Simi- 
larly, (E) can also be written as 
m 

(E) = £3c {|w^ (E’,x)|/|w^ (E)|> (E') exp {- 

m 4n m mm ° 

[l/A^(E*,y) + 1/X] dy}dx 4.29 

where (E^)dx = 411 (E") ^dx + Q« (E^dx 4.30 

m p3c m s(E") m 

Here, E" is the energy of the particle when it gives rise to a particle t 

m 

of energy E* and (E^) is the production spectrum of t nucleus directly 

m ™ 

through interactions of cosmic ray nuclei. At relativistic energies when 

Af % A^ Equation 4,29 sirq>llfies to 

(E) = /(l/s(E) + 1/X + 1/M}/(A» + X) 4.31 

In principle, one can make use of the above equations to calculate the 
ratio of Be/B nuclei as a function of energy be replacing s(E) by TyX/t for 
various values of t and compare it directly with the observations. However, 
the present experimental uncertainties do not permit such a comparison to 
estimate the exact age. Hence it is a general practice to compare the 
experimental results to the calculated ratio by considering only two 
possibilities, namely, (1) complete survival of l^Be which gives T«ty 
and (ii) complete decay of l^Be which gives T>>ty" 

In Figure 4.22 are plotted the available experimental data on Be/B 
as a function of mean energy (Casse et al., 1971; Brown et al. , 1973a; 
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Fig, 4,22 ; The observed ratios of Be/B at the top of the atmosphere are plotted as a function 
of kinetic energy per nucleon; for comparison are shown Curves A and B, the expected ratios for com- 
plete survival and decay respectively of ^^Be using steady state model. 
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0*Dell et al., 1973; Smith et al. , 1973a; Webber et al. , 1973a; jullusson, 

1974) . Also shown in the figure are two curves A and B representing the 
calculated values of the ratio for complete survival and decay respectively 
of the ^°Be produced. On the basis of a general inspection of the data 
points and the associated errors in Figure 4.22, one is inclined to infer 
that no meaningful conclusion can be made about the life time of the 
radiation. Nevertheless, O'Dell et al (1971, 1973) have averaged the 
individual values obtained by some workers and claimed that there is 
substantive evidence to suggest a life time less than a few million years. 

But Raisbeck and Yiou (1973) have made a careful assessment of the various 
statistical and systematic errors associated with the experimental and 
estimated cross-sections used in these calculations, and have arrived at 
the conclusion that with the existing errors, no meaningful inference on 
the time of cosmic rays can be made from the study of the elemental 

abundance of Be and B . 

Other recent experimental data which lead: to ' small values for t. 
are: (a) Price et al (1971) from a study of the abundance ratio of nuclei 

with Z>85 to those with 70<z583 conclude that the life time of cosmic rays 
is <10^ yrs and possibly as low as 10^ yrs ; (b) since the nucleus ^3^11 
decays only by K-capture with a mean life on two million years, Reames (1971) 
has concluded from an analysis of the observed ratio of Cr/Mn as a function 
of energy, that the cosmic ray age is consistent with the confinement of 
cosmic rays in the Disk rather than in the Halo. 

While all attempts in the past to use ^^Be as a clock to measure the life 
time of cosmic rays have been confined to the study of the ratio Be/B, the 
ultimate aim will be to study the relative intensitives of the isotopes, 
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where the effect of the ^^Be decay will be maximal and will not be diluted 
as in the c^e of the ratio of Be to B nuclei. With the growing capability 
and sophistication of detector instrumentation, such intrinsic resolution 
needed for the efficient separation of individual isotopes of Be and B 
nuclei has just become a reality (Garcia-Munoz et al, 1973; Webber et al, 
1973a) . These first measurements have not yet been able to lead to defin*- 
itive result^s but have clearly demonstrated that in the near future we will 
have a reliable estimate of the cosmic ray life time, 

4, 7 Anti-Nuclei in Cosmic Raya 

It is well known fact that detection of antiparticles in cosmic radiation 
is perhaps the only direct proof to the existence of antimatter outside the 
solar system. After the early measurements of the east-west effect (Johnson, 
1938), which showed that cosmic rays composed of positively charged particles, 
the major attempts to look for anti-nuclei during late fifties gave an upper 
limit of 0.1% at energies less than a few hundred MeV/n (Aizu et al, 1961; 
Grigorov et al, 1961). Nonetheless, annihilation of antimatter with matter 
has often been suggested as an energy source in those astrophysical objects 
for which nucleosynthesis or gravitation may be inadequate. Such a process 
to occur on a large scale in the baryon-symmetric Universe has been invoked 
by Stacker et al. (1971) to explain the diffuse component of the background 
gamma ray spectrum. Further, Alfven (1965) proposed a model allowing the 
co-existence of matter and anti-matter within our own galaxy. As a result 
of all these, systematic search for anti-nuclei has been undertaken over a 
wide range of energies using magnetic spectrometers. In this section we 
briefly summarize the available experimental results and the theoretical 
estimates of the antiproton flux of secondary origin in the Galaxy, 
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4.7.1 Anti-Nuclei of Secondary Origin 

When high energy cosmic ray nuclei interact with interstellar gas, 
baryon pairs are also created along with other hadrons. Hence, even in 
the absence of anti-matter in the Universe, one would expect to detect 
anti-protons (p) in cosmic rays ; indeed the presence of positrons in cosmic 
radiation is known well for many years. Many theoretical estimates of the 
fraction of p in cosmic radiation have been made in the past, of which the 
early ones were only speculative in nature (Fradkin, 1956; Burb ridge and 
Hoyle, 1957; Hayakawa et al. 1958). With the availability of the experi- 
mental cross-sections for the production of p over a wide range of energies, 
more detailed calculations have been attempted to determine the equilibrium 
ratio of p/p as a function of energy using various propagation models for 
cosmic rays (Rosen, 1967; Shen and Berkey, 1968; Way land and Bowen, 1968; 
Suh, 1971; Chen, 1972; Gaisser and Maurer, 1973; Badhwar et al, 1973). 

The production spectrum Q_(E)dE per (cm^ sec) can be written from 

P 

Equation 4.19 as 

q_(E)dE = 4ln dE'j(E’) <T_(E') i{»_(E’,E)dE 4.32 

^ t 

P P P 

where E, is the threshold for the production of a baryon pair, which is 
about 6 GeV in p-p interactions and n is the number of interstellar gas 

_3 

atoms cm . From the available accelerator data, the cross-section o(E') 

for the production of p, Includlng the ^Itlplicity, can te approximated 

(Badhwar et al, 1973) as cj_(E) % 9.26x10 ^(E-C), where C = lO.SGeV is the 

P 

effective threshold energy; this relation seems to be valid up to about 
500 GeV. The energy distribution function is given by Shen and Berkey 



141 


(1968) as 

i|;_(E*,E)<iE » HE exp (-nE^/4<E>^/2<E>^ 4.33 

P 

where <E> is the mean energy of p, which is aE’ on the basis of scaling 
hypothesis and is according to statistical model. Knowing the 

production spectrum of p,. one can then calculate the equilibrium spectrum 
using any propagation model. This spectrum then can be compared with the 
equilibrium proton spectrum to obtain the ratio p/p as a function of 
energy . 

In Figure 4.23, we show the calculated equilibrium ratio p/p as a 
function of energy (Gaisser and Maurer, 1973; Badhwar et al, 1973) for a 
mean amount of 5 g, cm^ of matter traversed by cosmic rays in interstellar 
space. Since the effective threshold energy for production is about 
10 GeV/n, the flux of antiprotons of energy less than <E> decreased while 
proton flux continues to Increase, resulting in a rapid decrease of the 
ratio below about 2 GeV. Further, the tincertainty In this energy region 
is rather large and depends upon the exactness of the Equation 4.33 for 
values of E much smaller than <E>, the propagation model and the uncertainty 
in the equilibrium intensity of protons in interstellar space; however, the 
experimental upper limit of 3x10 ^ at energies <200 MeV (Apparao, 1967) is 
about a few order of magnitude larger than the calculated value. At higher 
energies, the calculation of Gaisser and Maurer shows that the ratio 
approaches an asymptotic value of about 4.6x10 . It has been pointed out 

by Badhwar et al that p/p ratio is dependent on whether the matter traversed 
by cosmic rays is during or after the acceleration of particles to cosmic 
ray energies. We also show in this figure the only one available upper 
limit above a GeV (Bogomolov et al, 1971) and one notes that this limit 
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needs to be improved by at least two orders of magnitude to make meaningful 
comparison with the theory, 

4.7.2 Anti-Nuclei of Primary Origin 

As is evident from the preceding section, the detection of antiprotons 
in cosmic radiation cannot unambiguously prove the existence of antimatter in 
cosmic space before one could separate those resulting from high energy 
interaction. Hence, most of the experiments carried out so far are aimed at 
detecting antinuclei of Z^2. In Table 4.5, we show the available results 
from such studies as suumarized by Garcia-Munoz (1973). One may conclude 
from this that the fraction of antimatter, if it exists in the Galaxy, is 
less than axlO^'* of the normal matter. Though it is interesting to set a 
more stringent upper limit to this ratio by a few orders of magnitude, the 
present result itself might perhaps indicate that in the baryon-symmetric 
universe, matter and antimatter might have separated from each other on an 
extragalactic scale. Would that then tempt us to look for anti-matter at 
energies >10^^ GeV, where one would expect antimatter to reach us from 
distant galaxies in the Universe? 
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Table 4.5 

Upper Limits to the Ratio of Anti-Nuclei to Nuclei In . Cosmic Radiation 


Nuclei 

Rigidity/ 

energy range 

Upper Limits 

Reference 

Z=2 

0.2-4. 3 

GeV/n 

1.4 

X 

10-^ 

Eavanson (1972) 

Z^2 

5-33 

GV 

2 

X 

10-“ 

Smoot et al (1973) 

Z^2 

33-100 

« 

1.5 

X 

10“^ 

Smoot et al (1973) 

Z^3 

4-125 

ft 

3 

X 

lO""^ 

Golden et al (1973b) 

Z^3 

4-10 


8 

X 

10"^ 

Golden et al (1973b) 

Z^3 

10-50 

II 

5 

X 

io“^ 

Golden et al (1973b) 

Z-3 

60-125 

M 

6 

X 

10“^ 

Golden et al (1973b) 

Z^2 

15-30 


9 

X 

10"^ 

Verma et al (1972) 

Z-2 

30-50 

M 

2.6 

X 

10-^ 

Venna et al (1972) 

Z^2 

50-100 

II 

7 

X 

io“^ 

Verma et al (1972) 

Z^6 

< 200 

MeV/n 

6 

X 

10“^ 

Ivanova et al (1968) 

Z^6 

5-9 

GeV/n 

7.5 

X 

CM 

t 

Greenhill et al (1971) 
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5. Cosmogenic Electtomagafetlc Radiations In the Galaxy^ 

5.1 Introduction 

The motion of cosmic rays within the confines of the Galaxy would 
necessarily imply their interaction with matter, magnetic fields and 
radiations existing in interstellar space. Some of these interactions 
lead to the emission of electromagnetic radiations which may fall anywhere 
In the spectrum ranging from the radio regime to the gamma ray regime. 
Though the rate at which cosmic ray particles lose energy in this manner 
Is exceedingly small, their integrated effect becomes of considerable 
significance under cosmic conditions of space and time. Furthermore, 
a distinguishing feature of these emitted electromagnetic radiations is 
that they originate from all regions of galactic space permeated by cosmic 
rays and hence lead to a diffuse glowing or sky background in the relevant 
frequency region. It Is also evident, that as viewed from the solar 
neighborhood, the intensity of these galactic background radiations as 
a function of celestial directions will exhibit a close correlation with 
galactic coordinates. One may also anticipate the possible existence 
of a background radiation of extragalactic origin. However, in contrast 
to the former, the characteristic feature of the extragalactic radiation 
will be its isotropic nature. Observationally , the two components will 
be superimposed on one another and it will be necessary to resolve the 
two contributions before they can be interpreted. 

Cosmic ray interactions in galactic space leading to backgroimd 
electromagnetic radiations are of three kinds: they are interactions with 

Interstellar matter, radiation fields and magnetic fields. The 
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intensity of such a radiation received by a terrestrial observer from any 
direction in the Galaxy, defined by the galactic longitude 1 and latitude b, 
can be written as 

L(l,b) 

I (E)dE = ^ e i//^(E,s)ds photons cm“^ sr“^ sec“^ 5.1 

where ij;^(E,s)dE is the emissivity at a distance s through any of the pro- 
cesses mentioned above, f(E,s) is the optical depth and L(l,b) is the total 
distance along the line of sight in the Galaxy. In the case of the galactic 
Disk, Equation 5.1 is valid for telescopes having a resolution 9^ such that 

tan 8 ^ z /L(l,b), where z is the thickness of the Disk. For larger 

m o * o 

values of 0 , (see Figure 5.1), one can rewrite Equation 5.1 by expressing 
m 

the intensity from the Disk as 

jp L(l,b) 

<|)(E)dE ^ ^ f cos b db f e ^ i|j^(E,s)ds photons. 

- eV2 s = 0 

m 

cm"’^ rad^^ sec”^ 5.2 

The value of L(l,b) can be obtained from Figure 5#1 as 
L(l,b) ^ L(l)/Cos b for L (1) tan b <z^/2 

= z /Sin b for L (1) tan b ^z^/2 5.3 

and L(l) = {R^ + Ro^ - 2Ro [Ro Sin^l - Cos 1 (R^ - Ro^ Sin^l) ] 5.4 

Using this general formulation, we shall estimate the intensity of various 
galactic background radiations arising through the respective physical pro- 
cesses. We will then compare them with observations and therefrom deduce 
information on the astrophysical conditions in the Galaxy. 
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5,2 Cosmic Ray Interactions with Interstellar Gas 

In principle there are many physical processes Involving cosmic ray 
interactions with Interstellar gas resulting In the emission of electro- 
magnetic radiations. To name a few: (1) Greenberg (1969) has discussed 

the radio emission lines produced by electron cascades following recombi- 
nation to highly excited states; (il) the characteristic x-ray lines from 
transitions following inner shell ionization of ambient nuclei by cosmic 
rays have also been studied (Gould and Burbidge, 1963; Hayakawa and 
Matsuoka, 1964; Lampton et al. , 1971; Verma, 1971); (ill) Silk and 
Steigman (1969) have estimated the doppler-broadened Lyman a-like radiations 
emitted by low energy cosmic ray nuclei while capturing electrons into 
excited states followed by the cascading to the ground state; and (iv) 
nuclear interactions can result in the excitation of the incident or 
target nucleus which while reverting to its ground state emits gamma rays. 

In all these processes , except perhaps (iv) , the contributions originate 
from interactions of very low energy cosmic rays, about which we know very 
little at present. 

Of the two processes which make significant and calculable contributions 
of gamma rays, the first relates to the decay of neutral pions created in 
nuclear interactions of cosmic rays of energy above a few hundred MeV/n 
with interstellar matter. The second arises through the bremsstrahlung of 
cosmic ray electrons in electromagnetic interactions with nuclei of inter- 
stellar atoms. These two will now be discussed in some detail. 
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5.2.1 Plonic gaiama rays 

When cosmic ray nuclei, having energies in excess of the threshold 
for pion production, interact with Interstellar matter, gamma rays result 
from the decay of neutral pions created in such collisions; a small part 
also arises from the decay of S hyperons and tio mesons . From a knowledge 
of the pion production spectrum, the interstellar cosmic ray spectrum and 
the matter density at s, one can write the eraissivlty function i{)(E,s)dE 
defined in Equation 5.1 associated with this process as 

CQ 

i|;^(E.s)dE - 87tiijj(s)dE J ===== J jp(Ep.s) 

£ TT El 

Here c = E + m^^/4E; Hj^Cs) is the neutral hydrogen density at s; jp(Ep,s) 
is the cosmic ray proton flux at s; and a^(E^,E^) is the cross -sect ion for 
the production of mesons defined in this case as 

"" I ^ ®^^k»^p^ Vl ,5.6 

1C ^ 

The summation is carried over all interactions of cosmic ray nuclei of the 
type k with interstellar nuclei of the type 1; f^ is the fraction of cosmic 
ray nuclei of the type k to that of protons, f^ is the relative abundance 
of interstellar nuclei with respect to hydrogen and g(E^,E^) is a function 
which takes care of the energy thresholds for different nucleus— nucleus 
collisions. It can be easily shown that the pionic gamma-ray spectrum 
should exhibit a maximum at an energy of about 67.5 MeV, corresponding to 
half the rest mass of the neutral pion. Extensive attempts made to 
evaluate Equation 5.5 by different authors (Stecher, 1970, 1971; Cavallo and 
Gould, 1971; Levy and Goldsmith, 1972) differ by a factor ^^2 in the integral 
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intensity of gamma rays above 100 MeV, and hence there could exist a signi- 
ficant uncertainty in the estimated differential spectrum. In Figure 5.2 
is shown the differential intensity of gamma-rays for an interstellar gas 

_3 

density of one hydrogen atom cm . In this figure, the spectrum below a GeV 
is taken from Stecker (1970) , which is then extended to higher energies with 
a spectral index of -2,6. The constant spectral index above a GeV is justified 
by the fact that the observed spectrum of gamma rays in the atmosphere has 
a spectral index of -2.6 (Stephens, 1970a; Anand et al, 1973b). It might 
also be mentioned that the East-West asymmetry of atmospheric gamma rays over 
low latitudes suggests the strong influence of the isobaric decay as the source 
for the gamma rays (Stephens, 1970b). 

In order to estimate the intensity of gamma rays towards any direction 
of the Galaxy, one has to substitute Equation 5.5 in Equation 5.1. Since 
absorption of gamma rays in interstellar space is negligible one can rewrite 
Equation 5.1 as 

i|; (E)dE 

\<E)dE = J n^(s) Pj, r(s) ds 5.7 

0 


where (E) is plotted in Figure 5.2; p t,(s) is the cosmic ray density as 
normalized to that observed in the neighborhood of the Earth. If one assumes 
a mean value for the cosmic ray density and the gas density, the above inte- 
gral reduces to (E)dE/4iT <n^ consequence if one observes a finite 

flux of cosmic gamma rays in the Galaxy, one can in principle estimate the 
quantity <n^ ^ L> along that direction. Since bremsstrahlung radiation 

leads to a similar deduction, we shall compare the calculated spectra with 
the observed data after dealing with the bremsstrahlung gamma rays as well. 



(E) IN PHOTONS/{cm'^ sec GeV) 
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5.2«2 Bremsstrahlung radiation 

When cosmic ray electrons Interact with Interstellar gas , they give 
rise to bremsstrahlung radiation whose emlssivity can be written as 

CO 

Ip (E,s) dE - 4tt njj(s) dE J j^(E\s) a^(E,E’)dE 5.8 

E 

where j (E) is the equilibrium spectrtim of galactic electrons as discussed 
in section 4.5; c^(E,E*) is the differential cross-section for the produc- 
tion of gamma rays of energy E from an electron of energy E* summed up over 
all nuclei of the ambient gas. At energies greater than a few GeV, where 
the interstellar electron spectrum is a simple power law, the emlssivity 
can be written as (Stephens, 1969) 

ip^(E,s) =» 2.8 X lO”^^ photons cm ^sec ^GeV ^ 5.9 

Comparing this with Figure 5.2, one finds that this contribution is an order 
of magnitude smaller than the pionic gamma rays. The bremsstrahlung spectrum 
below a GeV flattens since the equilibrium spectrum of cosmic ray electrons 
also flattens. 

In Figure 5.3 are shown the Integral spectra of gamma rays resulting 

from the interaction of cosmic rays with Interstellar matter in the direction 

of the galactic centre (Anand and Stephens, 1972), for a typical detector 

with 0 ^ 30®. The pionic (curve A) and bremsstrahlung (curve B) gamma ray 

m 

intensities have been calculated by assuming Pq effective 

thickness of z = 360 pc for the gas disk, which corresponds to a thickness 
o 

of 280 pc as defined by the half density points (Section 2.3); Curve 
C is the sum of Curves A and B. It is seen that the bremsstrahlung radi- 
ation contributes a small fraction of the total radiation above 100 



FLUX IN PHOTONS/(cm -sec rad) 


O SHARE ef al (1974) 

# KRAUSHAAR «t ol (1972) 
///zKNIFFEN et Ql (1973) 
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Fig, 5i3 : The observed integral intensities of gamma rays towards 

the galactic center are compared with the calculated spectra for a detector 
with = 30® ; Curves A and B respectively are the spectra of plonic and 
bremsstrahlung gamma rays, and Curve C is the sum of Curves A and B. 
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MeV but becomes significant at lower energies. It has been shown by Maraschi 
et al (1968) that the estimated bremsstrahlung radiation below 1 MeV is 
very much smaller than the observed radiation, even with extreme assumptions 
regarding the spectral shape of cosmic ray electrons at the corresponding 
energies. We have also shewn in this figure, the recently observed flux 
values towards the galactic centre (Kraushaar et al,1972; Kniffen et al, 

1973; Share et al, 1974), confirming the original findings of a line emission 
from the galactic plane by Clark et al, (1968); there also exist some upper 
limits in this energy region which are in agreement with the above results 
(Dahlbacka et al, 1973; Bennett et al, 1972) except that of Frye et al, 

(1971). Though one may have to correct for the different values of 6^ 
used in these experiments in order to compare with theoretical estimates 
(except for the detector of Kraushaar et al, for which 6^ = 30®), one can 
infer from this figure that the observed value is a factor of about two higher 
than the estimated one. It has been further shown by Anand and Stephens 
(1972) that the central region of the Galaxy does not contribute much to 
this total radiation. If the observed intensity of gamma rays has to be 
accounted for by interaction of cosmic rays with ambient matter, then it 
is necessary that the value njj(s)p^ j^(s)ds/L 2n^(0)p^^^ (0) in the 
direction of the galactic centre. 

Many suggestions have been put forward to understand this enhancement 
of the observed flux of gamma rays towards the galactic centre on the basis 
of cosmic ray interactions with the ambient gas (Stecker, 1969; Anand and 
Stephens, 1971; Ginzburg and Khazan, 1972; Strong et al, 1973; Black and 
Fazio, 1973, Stecker et al, 1974; Bignami and Fichtel, 1974). In effect 
all these suggestions indicate that there could be an increase of cosmic 
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ray density with or without a corresponding increase of matter density 
towards the inner region of the Galaxy because the observed flux value in 
the direction of the anticentre (Kraushaar et al, 1972) is in good agreement 
with the estimated value in hat direction. 


5.3 Cosmic Ray Interactions with Radiation Fields 

When cosmic ray electrons undergo invers e-Compton scattering with 
ambient photons in the Galaxy, they impart part of their energy to the soft 
photons thereby giving rise to galactic background radiation in the X- and 
Y^ray regions. The emisslvity function in this case can be written in the 



where n(e,s)de is the number of photons per unit volume of energy e at s; 
j(E',s) is the flux of cosmic ray electrons; and o(E,E',e) is the cross- 
section for the production of a photon of energy E by an electron of energy 


scattering an ambient photon of energy e and is, according to Ginzburg 
and Syrovatskli (1964), 


o(E,E’ ,e)dE 


4re^ EdE f 
4 



5.11 


Here and m are the classical radius, mass and Lorentz factor of the 

electron, respectively. The limits and in Equation 5.10 are defined 
by the energy band in which the ambient photon radiation exists and = 

mc^ (E/4e)^. 

Estimates of the inverse Compton spectrum from the Galaxy has been made 
by many authors in the past (eg. Felten and Morrison, 1966; Cowsik and 
Pal, 1969; Shen, 1969; O'Connell and Verma, 1969; Ipavich and Lenchek, 1970; 
Maraschi and Treves, 1970). In Figure 5.4, the Integral spectra of 
the Inverse Compton photons in the direction of the Centre, assuming n(e) 
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Fig. 5.4 : The calculated Integral spectra of gamma rays through 

Inverse Compton scattering of black body, starlight and submillimeter 
photons are shown by Curves A, B and C respectively, and compared with 
the observations. 
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and j(E’) are independent of s, and = 1, for a typical detector with 

e = 30® are given (Anand and Stephens, 1972); also shown here are the observed 
a 

flux values summarized in Figure 5.3. Curves A, B and C of Figure 5.4 are the 
inverse Compton spectra from the Universal Black Body photons at 2.7 K, the 
star light photons and the subrolllimeter radiations respectively; the energy 
densities for these ambient photons are taken to be 0.25 eV, 0.5 eV and 
6 eV respectively. Here one has to note that the possible existence of the 
submillimeter radiation first reported by Shivanandan et al, (1968) is not 
confirmed by the recent experiments (Williamson et al, 1973); furthermore, 
its existence is shown to be inconsistent with the gamma ray observations 
(Anand and Stephens, 1972; Cowsik, 1972), Thus, if one ignores Curve C in 
Figure 5.4, it is quite evident that the contribution of inverse-Compton 
process from the Disk, in general, is very much smaller than that due to cosmic 
ray interaction with matter (Section 5.2). It may be noted that in the above 
calculation, it is assumed that n(e) is independent of s, which may not be 
true for the optical photons in the Galaxy. Recently, Sullivan and Cowsik 
(1973) showed that if the optical luminosity is proportional to the mass 
distribution in the Galaxy, one could account for the enhanced flux of gamma 
rays towards the Centre. On the other hand, it is found that if one makes 
use of the observed infra-red radiation from the central region of the Galaxy, 
rhe inverse Compton process contributes only a small fraction of the observed 
flux (Anand and Stephens, 1972). 

5.4 Synchrotron Radiation 

Cosmic ray electrons in the Galaxy emit synchrotron radiation as they 
spiral along the weak interstellar magnetic field lines. This synchrotron 
emission is recognized to be the source of the galactic nonthermal background 
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radio noise, which carries with it the signature of the mean magnetic fields 
and the electron spectrum involved. Our problem therefore is to disentangle 
the information carried by the continuum radio emission by matching it suitably 
with the electron spectriim seen near the Earth and other astrophysical para- 
meters. Attempts to connect the interstellar electron spectrum and the cosmic 
radio continuum through galactic magnetic fields , have been made by many workers 
(Bierman and Davis, 1960; Sironi, 1965; Felten, 1966; Ramaty and Ligenfelter, 
1966a; Okuda and Tanaka, 1968; Anand et al, 1968a, b,c; Verma, 1968; Webber, 

1968; Alexander et al, 1970; Goldstein et al, 1970; Stephens, 1971; Burger, 

1971; Bulanov et al, 1972; Ciimmlngs, et al, 1973b). 

Since radio intensity is usually described in terms of the frequency of 

emission, one can rewrite Equation 5.1 as 
L(l,b) 

I(v) - ^ j e *^^Tp(v ,s)ds Watt, cm sr ^ Hz ^ 5.12 

0 

Here the emissivity function can be written as 


iKv,s) = 1.04 X lO”**^ / B (s) G(v,E’)jg(E’ ,s)dE* Watt. cm~^ Hz~^ 5.13 


Where B_l (s) is the perpendicular component of the magnetic field at s 
expressed in Gauss; j^(E' ,s) is the cosmic ray electron intensity where E' is 
expressed in GeV; and G(v,E*) is the power spectrum of the emitted radiation 
by an electron of energy E' (same as the function F(x) defined by Ginzburg 
and Syrovatskil, 1965). The optical depth f (v,s) can be written from 
Ginzburg and Syrovatskil (1965) as. 


f(v,s) = 


3.1 X 10 


16 


n? fj(s) 


I {17.7 + ln(T 3/2/v)} k.ds 

J i T 3/2 ^ 


5.14 


where n^ is the electron density at z = 0; is the temperature of the medium; 

f^(s) is the distribution of the medium along the line of sight; and is the 



- 159 - 


fraction of the volume occupied by the medium of the type i (summation over 
i is needed when the medium is a complex one) . We shall first summarize the 
observational data before we proceed to estimate the synchrotron spectrum 
using the above equation, 

5,4.1 Observational Data 

During the past decade or so, detailed and systematic radio surveys 
of the Galaxy have been carried out in a broad band of frequencies using 
wide angle and pencil beams. Making use of these surveys, radio spectra 
have been constructed for different galactic directions, of which the 
following are of interest to us: (i) the Anticentre, (ii) the Centre and 

(iii) the Halo. 

(i) The Anticentre 

The Anticentre has been defined here over a broad region corresponding 
to -C 140-190® and b 10®N-10®S. In order to deduce the radio spectrum, 
surveys with wide angle beams in the range 10—1407 MHz are used (Costain, 

1960; Pauliny-Toth and Shakeshaft, 1962; Turtle et al, 1962; Parthasarathy 
and Lerfald, 1965; Andrew, 1966; Purton, 1966; Bridle, 1967; Howell, 1970; 
Sironi, 1974; Webster, 1974). The radio spectrtjm thus derived for the Anti- 
centre in the range 10-600 MHz is shown in Figure 5.5 by Curve A. If the 

spectral form at any frequency is expressed as « v then it is evident 

from this figure that the spectral index ot varies from about 0,3 at low fre- 
quencies increasing gradually to about 0.8 at about a few 100 MHz and remains 
constant thereafter at least up to 1407 MHz (Webster, 1974). 

(li) The Centre 

In this case two particular directions are chosen towards the Centre, 

but clearly avoiding the Nucleus; they are a) 1 = 0° , b = 3.6®N and 3.6®S 

and b) 1 = 20° and 340® ; b = 0® . For the purpose of constructing the associated 
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Fig. 5,5 : Radio spectra in the direction of the Anti-centre (Curve 

A) and from the Matagalaxy (Curve M) , 
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Fig. 5.6 : Radio spectra in the general direction of the galactic 

centre but slightly displaced from the centre, in latitude (Curve Rj) 
and in longitude (Curve R^) ; the references to data points are given 
in the text. 
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radio spectrum, surveys made with pencil beams in the frequency range 80 - 
4080 MHz are used (Baldwin, 1955a; Denisse et al, 1955; Kraus and Ko, 1955; 
Piddlngton and Trent, 1956; Hill et al, 1958; Westerhout, 1958; Seeger et al, 
1960; Wilson and Bolton, 1960; Large et al, 1961; Mathewson et al, 1962; 
Braccesl and Vespiganl, 1964; Moron, 1965; Seeger et al, 1965; Koinesaroff, 
1966; Penzias and Wilson, 1966; Wielebinski et al, 1968). These flux values 
shown in Figure 5.6 are the mean values in the two directions b = 3.6®N and 
3.6°S in the case of and 1 » 20° and 340° for R^, except In surveys, in 

which data is available for only one. It is apparent from this figure that 
the radio spectra R^ and R 2 are similar in nature except that at frequencies 
greater than a GHz, there seems to be a suggestion of a flattening of the 
spectrum In the direction b * 0 which is presumably due to thermal radiations 
from the galactic plane. 

(ill) The Halo 

The halo region is defined over a wide region of the sky away from 
the plane of the Galaxy. In the case of the North Halo, the regions cen- 
tered around the halo minimum radiation (1 ^ 190°, b ^ 50°) have been 
chosen. Using the same surveys, as in the case of the Anticentre, the 
radio spectrum in this direction of the sky has been obtained at frequencies 
^ 10 MHz; at low frequencies, the works of Getmantsev et al. (1968) and 
Alexander et al (1969) have been used. These flux values plotted in Figure 
5.7 show that below a few MHz, the values decrease with decreasing frequency 
presumably as a result of absorption. A similar spectrum has been obtained 
towards the south galactic Pole (Yates and Wielebinski, 1966; Alexander et al, 
1969) and is shown in Figure 5.8. 




Fig. 5,7 ; Radio spectra in the direction of north halo minimum radiation. Curves H and D are 
respectively the spectra from the Halo and Disk after taking into account the self absorption, while 
De Is the emission spectrum in the Disk. The solid Curve A and the dash-dot curve are the expected 
spectra for two different models of interstellar gas. •• . , 
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It Is imperative that before attempting to interpret the background 
radio emission as of galactic origin, one should enquire whether there 
could exist a substantial component of metagalactic origin in the observed 
radio Intensities. Such attempts have been made (Shain, 1958; Turtle et 
al, 1962; Yates. and Wielebinski, 1966; Purton, 1966; Bridle, 1967) and 
the information available is also Included in Figure 5.5 as curve M. It 
can be seen that the metagalactic contribution at 20 MHz is about 20 ^ of the 
total radiation towards the Anticentre. Though there is some indication that 
this spectrum is perhaps steeper than the galactic radiation, an attempt to 
detect its contribution at 2 MHz revealed its near absence (Bridle, 1968). 
Since the latter observation indicates that the spectral shape of the meta- 
galactic component has probably become flat at low frequencies, one can justi- 
fiably Ignore its influence in the analysis of the galactic radiation (Daniel 
and Stephens, 1970; Stephens, 1971). 

Now we can proceed to Interpret the radio spectrum by combining it with 
our knowledge on cosmic ray electrons sampled near the Earth. 

5.4,2 The Energy Spectrum of Electrons in Interstellar Space 

Anand et al (1968a, b) first deduced a reliable and detailed energy 
spectrum of galactic electrons, as it exists in the near interstellar space 
by making use of the radio spectrum towards the Anticentre. In these calcu- 
lations it was assumed that (a) the emissivlty ip(v,s) is independent of s, 
and (b) the energy spectrum of electrons observed near the Earth above 5 GeV 
is well preserved in the near interstellar space; also single values for the 
mean magnetic field <Bjj> and the cosmic ray density ^ were assumed 

along the line of sight. Under these asstmptlons Equation 5.12 is evaluated 
through Equation 5.13 by taking an effective path length of 4 kpc towards 
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the Anticentre. It can be shewn that, since the second assumption requires 
an interstellar electron spectrum below 5 GeV, which preseimably should smoothly 
join the observed spectrum above 5 GeV, it is possible to assign a single 
value for <Bj_> and a unique smoothly varying electron spectrum to match 
the radio spectrum. The electron spectrum thus deduced for the Anticentre 
is shown as the radio emitting electron spectrum by curve A in Figure 4.19; 
the corresponding value of <Bj_> is 5yG. The uncertainty in the deduced 
electron spectrum increases with decreasing energy and is about a factor 2 
around 300 MeV increasing to very large value below 100 MeV (Cummings et al, 
1973b), due to the effect of interstellar absorption which dominates below 
a few MHz. Nevertheless, we feel that the smooth spectrum shown by Curve A 
in Figure 4,19 could be the most probable electron spectriim in Interstellar 
space because it is quite obvious from Figures 5.5, 5.7 and 5.8 that the 
radio spectral index changes gradually from 0.8 to 0.3 In the frequency 
range 200-10 MHz, where there is no absorption effect. 

Having derived the electron spectrum towards the anticentre region, it 
is natural to enquire whether this spectrum could also explain the obser- 
vations in the other galactic directions- In the direction of the Centre, 
it is found that the same electron spectrum could reproduce the observed 
radio spectra (curves and of Figure 5.6) with a single value of <B_jj> = 

7.2yG for the line of sight 1 = 13 kpc and L = 22 kpc. In the same 

Kl K2 

manner, the halo spectra can also be well explained above a few MHz with 
a magnetic field <Bjj> = 2yG for = 11 kpc towards the north halo minimum 
(Figure 5.7) and Ljj = 15 kpc towards the south Pole (Figure 5.8). From 
this, one may infer that the spectral shape of cosmic ray electrons remains 
the same in all regions of the Galaxy. 
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(il) Magnetic Fields in the Galaxy 

We have noted already that If * 1 applies to all regions in the 

Galaxy, then one can obtain a unique value for the mean perpendicular com- 
ponent of the magnetic field along any galactic direction from an analysis 
of the radio continuum in the Galaxy. Such values of <Bj> in the direction 
of the Anticentre and Centre in the Galaxy are about 5 and 7 pG respectively. 
This variation in the magnetic field strength can either be real or be a 
reflection of the effect due to p_ _ being i^X st all regions, combined with 
a possible variation of the magnetic field. If one assumes random field 
orientations, the above results indicate that the total magnetic field varies 
from about 6pG to about 9pG in the Disk. As for the Halo, the derived 
magnetic field has a value of about 2.5yG; the uncertainty in this value 
depends upon the value of p_ _ and the extent of the cosmic ray halo. It 
may be noted that the magnetic field strength deduced for the Disk here seems 
to be much larger than the value of about 3yG, commonly used by astrophysicists 
In order to resolve this difference, the following explanation has been 
advanced by Cowsik and Mltteldorf (1974) . The fluctuations in the magnetic 
field strength, whose magnitude is found to be of the order of the average 

field itself <AB> B (Jokipil and Lerche, 1969), can arise from the 

nns o 

compression and rarefaction of the ambient field by gas motions, which are 
coupled to the lines of force in interstellar space. This process can be 
accompanied by a correlated variation in the density of cosmic ray electrons 
due to betatron process. As a consequence, the synchrotron emissivity varies 
non-linearly (Equation 5.13) resulting in an enhanced mean emissivity depending 
upon the value of this could be as large as a factor 'v^lO for <AB>^^ 

B . 
o 
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5.4«4 Mean Density of Cosmic Rays In Different Regions of the Galaxy 

It has been demonstrated above that a consistent picture can be con- 
structed to explain satisfactorily the observed non-thermal radiation from 
different regions of the Galaxy, including the Halo* At the same time it 
is instructive to examine how much variations can be permitted in the value 

of p used. One way of doing this is to assume the same value of <Bj_> 

C * R* 

deduced for the anticentre direction to apply for the direction of the 
galactic centre and then to estimate the value of the latter 

direction. If one does this it leads to a mean cosmic ray density 
towards the Centre. Thus one finds that the cosmic ray density averaged over 
the direction as the Centre is unlikely to exceed by very much more than a 
factor of 2 in comparison with the outer regions. It has also been shown 
by Anand et al (1968b) that the cosmic ray density in the Halo cannot be 
lower than in the Disk by a factor as much as 2 or more. However, if cosmic 
rays are confined to the Disk and they slowly leak into the halo region, where 
the mean free path could be much larger than in the Disk, the density of 
cosmic rays in the Halo has to be very much smaller than in the Disk (Laster, 
1964) . This apparent contradiction arises from the fact that we have assumed 
a uniform distribution of cosmic rays. Since this is a simple minded picture 
one should keep an open mind for possible large scale variations of the cosmic 
ray intensity in localized regions as well as a possible general gradient 
from the Nucleus towards the Poles or to the periphery of the Disk. 

5.4.5 Some Inferences on the Properties of Interstellar Medium 

Since the optical depth depends upon the physical properties of the 
interstellar medium, it has been pointed out (Stephens, 1971) that the 
observed radio spectrum in the direction of the Halo, and the electron spec- 
trum derived earlier, can be used to distinguish between various models of 
interstellar medium. In Figure 5.7, the estimated emission spectrum from 
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the Disk towards the north halo minimum is shown by Curve D^, for a model 

in which cold Interstellar clouds are postulated to be immersed in the 

hot Inter cloud medium. Curve D in this figure is the expected spectrum 

when absorption is taken into account by setting in Equation 5.14 for the 

—3 

intercloud medium the values n^ * 0.04 cm , = 9000®K, kj^ = 1 and fj^(s) 

= exp (“S^ sin^b/2a^) , where = 160 pc; for the clouds one has n^= 0.02 
cm''^, = 90°K, k^ = 0.02 and f^Cs) = exp (rrs^sin2b/2a2^) , where 

= 120 pc; these parameters are similar to those derived by Field et al 
(1969) from other consideration. Curve H is the radio spectrum from the 
Halo after taking into account the absorption in the Halo by setting in 
Equation 5.14, n = 0.005 cm“^, T = 10^ ®K, f(s) *= k = 1. Curve A in the 

same figure is the sum of the radiations calculated for the Disk (Curve D) 
and for the Halo after taking into account the absorption of halo radiation 
in the Disk. The dash-dot curves are the calculated spectra for another 

_3 

set of parameters with = lOOO^K, = 50®K and n^ = 0.03 cm , the other 
parameters remaining the same, as deduced by Hjellming, Gordon and Gordon 
(1969). One finds that while curve A is in good agreement, the dash-dot 
curve is at variance with the observations. It is also found (Stephens, 
1971) that even in the absence of clouds, the fit obtained by curve A 
demands a temperature T^ > 4000° K, thereby ruling out the possibility that 
the temperature of the intercloud medium could be as low as 1000°K. A 
similar estimate of the radio spectrum towards the south Pole using the 
first set of parameters confirms the earlier deductions as one can see from 
Figure 5.8. 

5.4.6, The Galactic Halo 

The existence of a galactic halo of near spherical shape, glowing in 
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the radio region by synchrotron radiation emitted by relativistic electrons 
spiralling in weak magnetic fields existing therein, was first postulated 
by Shklovsky (1952). This suggestion was made credible by Baldwin (1955b), 
who on the basis of cosmic background radio observations, attributed a 
diameter of 20-30 kpc for this halo. The existence of the radio halo, 
which received almost universal support from observational radio astronomers, 
astrophysicists, and cosmic ray physicists for over a decade, and was some- 
times acclaimed as one of the most important discoveries of our times, is 
now being increasingly questioned during recent years. The arguments so 
far advanced against its existence have been drawn from recent observations 
on the galactic background radio emission which show structural features. 

If one considers that these structural features, as revealed by pencil beam 
surveys are localized phenomena within the Disk and that the general back- 
ground is indeed composed of the integrated emission from such features, 
then there can be no evidence for the existence of the radio halo. However, 
this extreme hypothesis can be understood only if the general magnetic field 

in the galaxy is exceedingly small («10 gauss) and that the emission takes 

-5 

place essentially in discrete regions where the magnetic field (>10 gauss) 
and the cosmic ray density are very high; this however seems very difficult 
to accept from the study of optical polarization (Mathewson and Nichols, 

1968) . Granting then that the radio emissivity can be considered to be nearly 
uniform in the Disk, which seems reasonable, it has been shown (Anand et al, 
1968c; Daniel and Stephens, 1970) that the present observational evidence 
cannot rule out the existence of a radio Halo from the following arguments; 

(i) It is evident from Figures 5.5 and 5.7 that the combined contribution 
from the Disk and the Metagalaxy towards the minimum halo direction may not 
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exceed 50% of the total radiation observed in that direction; the excess 
radiation should therefore have its origin from the Halo; (ii) It can 
also be shown that the ratios of the excess radiation towards different 
directions in the Halo to that towards the minimum radiation are consistent 
with a volume mucJi thicker than the Disk. Thus one is faced with the 
situation that the cosmic ray observations suggest either no Halo or a leaky 
halo, while the radio observation cannot at present rule out the existence 
of a radio halo. 

6. The Role of Cosmic Rays in Galactic Dynamics : 

We have so far considered in detail problems relating to the propagation 
of cosmic rays in galactic space and the many consequences arising from the 
particulate nature of this radiation. At the same time one recognizes that 
cosmic rays in bulk constitute a relativistic gas with its own contributory 
pressure and energy density in galactic space, comparable to those due to 
magnetic fields and thermal motion of gas in interstellar space. Though 
these facts were generally known for a long time, it is only during recent 
years their real significance has been increasingly realized resulting in a 
quest to investigate the consequences of the existence of the cosmic ray 
gas in galactic space. 

The role that cosmic rays play in galactic dynamics may be briefly 
stated as follows. The interstellar gas is confined to the Disk by gravi- 
tational potential resulting from the distribution of stars perpendicular 
to the plane of the Galaxy; in this the gas has its own random motions. At 
the same time, the magnetic field is tied to the gas because it is partly 
ionized; cosmic rays in turn are tied to the magnetic field lines. All 
these forms of energy exert pressure against the gravitational force. Hence, 
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cosmic rays play a,n important role in the hydrostatic equilibrium and 
gravitational stability of gas in the Disk. Further, the subrelativistic 
cosmic rays, during their propagation, lose their energy rapidly through 
ionization before they are thermalized; as a consequence they heat and 
ionize the interstellar gas. The stability and thermal phase of the gas 
then depend on the density of low energy cosmic rays existing in inter- 
stellar space. In the ensuing sections, we present some of the basic 
approaches that have been tried so far to demonstrate these effects, and 
the reader is cautioned that no attempt has been made to review this 
subject . 

6.1 Hydrostatic gquilibrjuni of the Gaseous Component of the Galaxy : 

The thickness of the Disk, composed mainly of stars, is maintained by 
the gravitational acqeleration g(z) , resulting from the distribution of 
stars and gas, against the random motion of stars and can be described by 
the hydrostatic equilibrium equation 

In P (z) =» 6.1 

S dZ 8 aZ 

where <v'^(z)> is the root-mean square of the random velocities of stars, 
s 

which is about 18 km. s^c ^ in ths neighborhood of the Sun (Wooley, 1965), 
pg(z) is the star density and <|i(z) is the gravitational potential. The 
two dimensions in the plane of the Disk, in which the centrifugal force of 
rotation of the Galaxy is balanced by the gravitational attraction of the 
Nucleus, are not directly relevant here and hence the hydrostatic equation 
relating to the plane may be written as (Parker, 1969) 

^ {p(z) + B^(z)/8tt + P(z)> - p(z)g(z) 


6.2 
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Here, p(z) « <v^( 2 )>p(z) Is the gas pressure, B(z) Is the magnetic field 
strength and P(z) is the cosmic ray pressure, which is one third of its 
energy density. 

In the above equations, it is implicitly assumed that the external 
pressures due to magnetic field and cosmic rays existing outside the Disk 
and the weight of halo gas are negligible. If we make a further assumption 
that the distribution of magnetic field strength and cosmic ray density 
in the z-plane is similar to that of gas namely, B^(z)/p(z) = b2(o)/p(o) 
and P(z)/p(z) = P(o)/p(o), then Equation 6.2 reduces to 

p(z)} = g(z) 6.3 

where (o) /{ Sirp (o) } + P(o)/p(o). This equation can be 

integrated to get the distribution of gas as 

1 z 

p(z) = p(o) exp { Q 2 S g(z)dz} 6.4 

^ o 

The acceleration parameter in this equation can be determined as a function 

of z from an analysis of star counts and radial velocities of K~giants using 

Equation 6.1; according to Oort (1960) the values of |g(z)| at z = 0, 100, 

200, 300, 400 and 500 pc are respectively 0.0, 2.5, 4.3, 5.4 6.2 and 6.8 x 
-9 -1 

10 cm. sec . In Figure 6.1 is shoxm the distribution of gas p(z)/p(o) 
as a function of z. The dashed curve is the experimentally determined 
distribution by Schmidt (1957) corrected for the galactic scale Ro » 10 kpc 
(Section 2.1). The solid line is the calculated curve using Equation 6.4 
(Kellman, 1972a) for Q = 9.84 km, sec . One can see the good agreement 
between the calculated and observed distributions within about 250 pc and one 
obtains 

p(o)<v^(z)>+b2(o)/8tt+P(o) = 9.7 X lO^^p(o) 


6.5 
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Fig, 6.1 : The observed distribution of gas perpendicular to the 

galactic plane (dotted curve) is compared with the calculated one (solid 
curve) . 
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From the above equation it Is possible to derive the value for any one 
parameter, If one knows the other two. However, in practice one notes that 
at present, all the three quantities are known with comparable degrees of 
uncertainty with the result that there is no special merit in deducing the 
value of one by using the other two. Nonetheless, to get a feel for the 
magnitude of these quantities one can assume that they contribute equally 
and thereby deduce the following values: <v^>^ 5.7 km.s”^, P 5 ,2 dynes, 

cm and B ^ 3.6yG for p(o) =» 1 hydrogen atom, cm . These numbers can 
be now compared with those available from observations. 

From a study of the velocity profiles of 21 cm hydrogen line, it is 
found that the velocity dispersion for clouds (Westerhout, 1957) as well 
as for intercloud medium (Helles, 1967) is about 6 km.s this is in good 
agreement with the value deduced above. In case of cosmic ray pressure, 
though the equipartition value is consistent with that observed in the 
neighborhood of the Earth, one suspects that the true pressure in inter- 
stellar space is likely to be much larger. Similarly, the mean magnetic 
field obtained from the radio brightness distribution (Section 5.3.3) is 
also larger than the value deduced above. This could imply that the internal 
pressure might be numerically larger than the gravitational force to achieve 
hydrostatic equilibrium. On the other hand in our assumptions, we have 
neglected the external pressures, such as those due to the presence of gas, 
cosmic rays and magnetic field in the halo region, which could in principle 
compensate this apparent higher internal pressure. One also notices from 
Figure 6.1 that the theoretical curve deviates from the observed one at 
large values of z. This again is due to the assumption that the distribution 
of p, P and B are identical, which is not true because we know that the radio 
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disk is much thicker than the gas disk (Section 2,1). Hence Equation 6.2 
needs to be solved more rigorously and the final fit to the observed data 
would give information on the possible distribution of cosmic rays and 
magnetic fields away from the galactic plane. 

6.2 Stability of Self Gravitating Gas and the Formation of Clouds 

The condition for the stability of self-gravitating gas can be under- 
stood from the following concept (Jean, 1928). Let us consider a uniform 
distribution of gas, in which a small condensation at a point P is 
accompanied by a compensating rarefaction at a nearby point Q, This 
disturbance results in a gain in the total internal energy, which then 
travels about with the velocity of sound until it is dissipated by viscosity 
in the form of heat energy. Such a process gives rise on one hand, to an 
increase in the thermodynamical energy of the gas but on the other, it causes 
a decrease in the gravitational energy. If P and Q are sufficiently far 
apart, the decrease in the gravitational energy becomes numerically greater 
than the increase in the thermodynamical energy and as a result the medium 
acquires kinetic energy at no expense; this will continue to grow indefinitely. 
This form of instability resulting through displacements in which condensation 
and rarefaction occur in pairs at sufficiently distant points in the medium 
is called the "Jeans gravitational Instability". 

For a uniform gas with a uniform gravitational potential throughout, 

Jeans (1928) has showed that the hydrodynamical equation relating to the 
gas motions described above can be written as 

= 4tt Gp s + C2 v^s 6.6 

where s is the condensation parameter 6p/p, C is the speed of sound in the 
medium and G is the gravitational constant. Considering now a pure wave 
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motion along the X-axis, with s being proportional to Cos (2nx/X) , the 
above equation becomes 

= [4irGp - c2]s 6.7 

where X is the wavelength. The solution to Equation 6.7, in which s is 
proportional to exp (+ Iqt) , represents the wave motion, where q^ = (2 tt/X)^C^ 
- AirpG, with a velocity of propagation 

V = [C^ - (X/2tt)^ 4xGp]^ 6.8 

One can see from this equation that, when the gravitational term is omitted, 
the wave motion which travels with a uniform velocity C is independent of X. 
The restoration of the gravitational term decreases the velocity of wave 
propagation, and since this term is proportional to X^, it is ineffective 
for short wave lengths. As we slowly increase the value of X, Equation 
6.8 becomes zero and then becomes imaginary. For such values of X there 
can not be any normal propagation of waves and the initial condensation 
and rarefaction will increase exponentially leading to unstable motion. 

The value of X for marginal instability, i.e. when the wave motion ceases 
to exist (according to Equation 6.8) is = c/iT/Gp. The radius of gas Rj 
under marginal instability is X/2 which is about 0.4 kpc for <v^>^ Rj 5 km. 
sec” and p 1 hydrogen atom, cm . Thus, under simple Jeans* instability 
the interstellar gas will start forming clouds in scale sizes larger than 
0.4 kpc in radius. The time scale for the onset of instability to the 
formation of clouds is about q ^ 'v 10® yrs. 

The effect of the inclusion of magnetic fields and cosmic rays in the 
self gravitating gas is to increase the internal energy of the system. From 
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an analysis of hydrodynamic equations relating to this gas-field system, it 
has been shewn by Kellraan (1972b) that the radius of marginally unstable 
gas increases by a factor (1 + a + 8)^, where a and 3 are the ratios of 
magnetic pressure to gas pressure, and cosmic ray pressure to gas pressure 
respectively. When these parameters are comparable, the instability starts 
at wavelengths corresponding to R == Rj. It has been further shown by 
Ledoux (1951) that for plane parallel, non-rotating gas layer such as the 
interstellar gas, the radius of marginally unstable gas at the plane of 
symmetry is larger by an additional factor Thus, for interstellar gas 

the value of R R,» 

The above treatment shows that the self gravitating gas tends to be 
stable against the simple gravitational instability under the influence of 
magnetic field and cosmic rays which add to the internal pressure. This 
situation is true for the interstellar gas only at the plane of symmetry 
at z = 0, where the magnetic field is parallel to the plane and is uniform. 
However, it was pointed out earlier (Section 3.2.3) that the instability 
in the gas field system inflates the magnetic field lines and the gas tends 
to slide along the field lines under the gravitational force. In Figure 
6,2 is shown a schematic diagram of an idealized vertical section through 
the Disk, taken along the lines of force (Parker 1968b) . When the lines of 
force are perturbed as shown by the sinusoidal line, the gas tends to slide 
down into the troughs along the field. This burdens the low regions causing 
them to sink further and unburdens the high places permitting them to expand 
upward. Similarly, the cosmic ray pressure is lower than the ambient pressure 
at high places contributing to the inflation and expansion of the field. 

The following simple treatment by Parker (1969) Illustrates this process 
clearly , 
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the gas field system perpendicular 
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Let the magnetic field ^ be along the y-direction and the gravitational 
acceleration be in the negative z-direction. Consider two elements of mass 
m^ and suspended on a line of force. The weight m^ distorts the field 
and changes the field from (0, B, 0) to ® position 

of m^; the slope of the lines of force through m^ is approximately 
As a result m 2 is pushed towards the negative direction of y along the 
magnetic field with a force In case of the magnetic field being 

pervaded by a tenuous background plasma with negligible pressure and weight, 
the distorsion of field, at a remote position (x, y, z) due to mass m^ at 
the origin, exert a force on given by 

F(r) — rajm2g2(z)y[l + (y^ + -ZxV/Cx^ + 6.9 

where r^ = x^ + y^ + z^. This can be compared with the force due to direct 
acceleration between and m^, namely ^" 0 ( 2 ^) ” Gm^ra^y/t^. The first term 
in Equation 6.9, which is zero at the plane of symmetry where g(o) = 0, 
increases numerically with |z| and exceeds very much larger than the simple 
self gravitating force. Hence, the total effect on mass m^ is proportional 

to G “ g^/B^ instead of G. Parker has further shown that the time scale 

h 7 

between the onset of instability and the formation of cloi^s is (A/g) '^10' 

yrs, where A is the scale height of the gas in the z-direction; this value 
is very much smaller than for the simple Jeans* instability. 

One can see from Equation 6.9 that, while the second and third terras 
average to zero over any plane y = constant, they are significant for large 
values of y. In the limit of large y. Equation 6.9 reduces to 

F(r) = -mjm 2 g ^(z)I(z^-x^)/(z^-Hx^)^] /B^ 6.10 

This force is attractive when z>x. The gas, which is displaced in z below 
or above a cloud is powerfully attracted towards the cloud while those dis- 
placed in X are repelled. In this picture, the initial instabilities may 
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grow over scales of the order of 100 pc along the lines of force and 10 
to 100 pc across the lines of force. Further, the state of motion of gas 
clumps would continue until the gas density becomes sufficiently large for 
self gravitation to take over the repulsive force described above and the 
internal forces resulting from thermodynamical, magnetic and cosmic ray 
pressures. 

6.3 Thermal Equilibrium of Interstellar Gas ; 

The observed physical properties of interstellar gas have been inter- 
preted (Section 2.3) as due to the gas being composed of cool, dense clouds 
embedded in a hot, rarefied intercloud medium. In this context, Hayakawa 
et al (1961) were the first to suggest that low energy cosmic rays are 
quite effective in heating the interstellar medium. The heating can proceed 
through primary ionization of gas atoms by the direct interaction of low 
energy cosmic rays and through secondary ionization by electrons liberated 
in the primary ionization process. The cooling of gas, which is very 
essential to maintain thermal equilibrium, is mainly through the de-excitation 
of atoms, both ionized and neutral, and free— bound and free— free emission of 
free electrons. The evaluation of the equilibrium phase of the gas is a 
complex mathematical exercise (Field, 1970; Dalgarno and McCray, 1972) and 
hence in this section we only Indicate the results from one such calculation, 
and summarize the various processes through which cosmic rays could contri- 
bute to the heating of interstellar gas. 

It is quite obvious that for a given rate of heating, the equilibrium 
free electron density and the temperature of the medium depend upon the 
density of the gas. In Figure 6.3, we show the result from a typical calcu- 
lation (Field et al, 1969), in which the equilibrium temperature is plotted 
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Fig * 6 > 3 ; The thermal state of the Interstellar gas as a function 
of gas density. 
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against the number density of interstellar gas atoms. One can see from 

this figure that the thermal phases, which are not very sensitive to the 

change in the number density, around 10**^**K and below 10^®K are stable. 

—1 —3 

These phases occur at densities n 10 atoms cm (which is typical of 

—3 

the intercloud medium) and n 10 atoms cm (which is characteristic of 

the clouds respectively; the intermediate state is obviously unstable. 

It has also been shown by Field et al (1969) that the equilibrium tempera- 

+ o 

ture critically depends upon the abundance of cooling agents like C , 0 , 

Si and Fe in the gas, and that their depletion in dense clouds as a 

result of accretion of these trace elements in the grains can Increase 

the temperature of the cold phase considerably. 

The above calculations on the thermal equilibrium of interstellar gas 

suggest that the ionization rate required is 'v 10 sec Since the rate 

of loss of energy through ionization increases as the energy of the particle 

decreases, only very low energy cosmic rays below 50 MeV/n are of importance; 

here heavy nuclei contribute substantially. It has also been suggested (Silk 

and Werner, 1970) that low energy x-rays in the region 100-250 eV, which are 

effectively absorbed in the interstellar medium, could also be an important 

source of heating and ionization. However, the recent observations by 

Copernicus satellite (Rogers on et al, 1973) indicate that low energy cosmic 

rays and x-rays cannot be the main source of heating of interstellar gas, 

because they lead to higher ionization states which are not observed, and 

-18 -1 

set a limit of '^/ 10 sec for the rate of ionization by these processes. 

It may perhaps be relevant to consider quantitatively the mechanism suggested 
by Wentzel (1971), in which heating of gas results from the dissipation of 
hydromagnetic waves created by cosmic rays during their propagation in 
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Interstellar space (Section 3.2.2); however, since bulk of the cosmic rays 
have energies 1 GeV, the heating by this process is associated with very 
little ionization. 
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